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PLUME  MECHANICS  AND  AEROSOL  GROWTH  PROCESSES 


I.  OBJECTIVES  OF  RESEARCH 

Work  under  this  contract  was  intended  to  add  to  the 
technological  base  of  the  U.  S.  Army's  programs  in  aerosol  technology, 
smoke  and  protection,  as  well  as  to  contribute  to  fundamental  knowledge 
of  aerosol  physics  and  chemistry.  In  this  study,  the  following  areas 
have  been  investigated  : 

(1)  Formation  of  liquid  and  solid  particles  by  nucleation  and 
growth . 

(2)  Influence  of  smoke  plume/puff  dispersion  on  aerosol 
growth . 

(3)  Nonlinear  interactions  of  laser  beams  with  aerosols. 

(4)  Aerosol  filtration. 

Areas  (1)  ana  (2)  above  deal  with  key  problems  in  smoke 
obscuration.  The  efficiency  of  an  aerosol  obscurant  for  electromagnetic 
radiation  depends  on  many  factors  including  particle  size,  shape,  com¬ 
position  and  concentration.  For  an  obscurant  propagating  in  the  atmos¬ 
phere,  particle  size,  shape,  composition,  and  concentration  are  deter¬ 
mined  by  the  processes  of  particle  generation  and  growth  as  well  as  by 
advection,  dispersion  and  other  processes  of  atmospheric  motion.  While 
some  qualitative  features  of  these  various  processes  are  recognized,  our 
knowledge  is  incomplete.  Improvements  in  aerosol  obscurant  technology 
must  depend  in  part  on  development  of  the  basic  knowledge  of  the  various 


processes  cited  above. 


Most  aerosol  obscurants  are  generated  by  the  condensation  of 


atomic  or  molecular  species  initiated  by  a  nucleation  process.  The 
description  of  this  process  remains  today  as  one  of  the  unsolved 
classical  problems  of  physics.  Subsequent  to  their  formation  by 
nucleation,  the  stable  particles  grow  by  coagulation  and 
condensation/ evapora t ion  processes.  For  liquid  particles  of  a  single 
chemical  species  which  coalesce  on  collision,  these  processes  for 
restricted  homogeneous  systems  can  be  described  quantitatively.  However, 
even  in  this  limited  case,  this  has  only  recently  been  done.  For 
particles  that  do  not  coalesce  on  collision,  no  general  theory  is 
available.  Such  particles  will  form  larger  particles  with  complex 
morphology  -  chains,  branched  structures,  random  aggregates,  etc.  -  whose 
description  in  the  context  of  a  dynamic  model  has  not  yet  been  achieved. 
When  the  aerosol  is  formed  from  a  mixture  of  chemical  species,  the 
problems  cited  above  for  single  component  aerosols  are  compounded.  The 
study  of  the  evolution  of  aerosol  mixtures  is  only  in  its  initial  stages. 

When  aerosol  obscurants  are  formed  in  the  atmosphere,  the 
formation  and  growth  processes  discussed  above  will  occur  in  a  turbulent 
environment.  First  these  processes  take  place  in  turbulence  produced  by 
the  obscurant  generation  device.  Then,  atmospheric  turbulence  takes 
over,  and  the  subsequent  evolution  occurs  in  this  environment.  The 
description  of  turbulence  also  remains  as  one  of  the  unsolved  classical 
problems  of  physics.  While  qualitative  descriptions  of  aerosol  evolution 
in  a  turbulent  environment  are  available  in  special  cases,  no  general 
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quantitative  description  is  currently  available.  Only  empirical 
descriptions  of  limited  validity  can  be  found  in  the  literature.  The 
evolution  of  an  aerosol  in  an  atmospheric  plume  is  a  complex  process 
whose  average  characteristics  have  only  recently  been  studied 
quantitatively. 

Practical  application  of  work  under  area  (3)  is  found  in  the 
problems  of  transmission  of  intense  laser  beams  through  the  atmosphere 
and  in  possible  use  of  aerosols  as  protective  screens.  If  one  wishes  to 
propagate  an  intense  laser  beam  in  the  atmosphere,  the  presence  of  water 
aerosols,  fog  and  cloud  droplets  poses  a  difficulty.  In  this  case  an 
understanding  of  rapid  aerosol  evaporation  by  laser  beams  is  essential. 
Laser  beams  are  totally  absorbed  by  plasmas  and  scattered  by  particle 
explosion.  These  phenomena  present  possibilities  for  the  use  of  aerosol 
as  protective  screens  against  intense  laser  beams. 

Area  (4)  is  concerned  with  the  problems  of  protecting  personnel 
by  use  of  mask  air  filters  against  chemical  and  biological  aerosol  threat 
agents.  In  order  to  ensure  such  protection,  it  is  essential  to 
understand  filtration  mechanisms  for  fibrous  and  membrane  filters  and 
factors  which  could  contribute  to  enhanced  penetration  of  filters  by 
aerosol  particles. 

Work  under  this  contract  has  addressed  important  aspects  of 
these  areas  of  application.  In  the  following  report,  the  research 
performed  is  outlined.  Full  details  are  available  for  the  most  part  in 
the  various  publications  resulting  from  this  research. 
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II.  SUMMARY  OF  RESULTS 


Highlights  of  major  accomplishments  from  the  research  are 
presented  here  along  with  a  listing  of  publications  arising  from  this 
work.  This  section  is  followed  by  a  more  detailed  exposition.  The  major 
accomplishments  are  listed  in  approximate  chronological  order. 

A.  Major  Accomplishments. 

Under  contract  from  the  Chemical  Research  Development  and 
Engineering  Center (CRDEC)  some  major  accomplishments  of  our  work,  in 
approximate  chronological  order,  are  as  follows: 


*  A  theoretical  basis  and  numerical  schemes  have  been  developed 
for  simulating  the  atmospheric  dispersion  of  plumes  of  volatile  aerosol 
such  as  diesel  smoke  plumes.  Computer  code  for  this  simulation  has  been 
transferred  to  CRDEC  for  their  use  in  planning  and  evaluation  as  a 
replacement  for  inaccurate  Gaussian  plume  models.  See  publications  CRDEC 
1 ,  2  and  3 • 

*  A  theoretical  basis  and  computer  code  have  been  developed  for 
dispersing  atmospheric  plumes  emitting  infrared  or  other  radiation  as  a 
result  of  chemical  or  thermochemical  reactions.  See  publication  CRDEC  4. 


been  carried  out  showing  the  approach  to  asymptotic  limit  distributions 
in  the  continuum  and  free  molecule  regimes.  See  publication  CRDEC  5. 

*  The  first  simulations  of  the  dynamics  of  the  Ostwald  ripening 
process  with  simultaneous  coalescence  have  been  carried  out  showing  that 
scaling  arguments  are  insufficient  to  establish  the  dominant  asymptotic 
growth  mechanism  in  the  continuum  regime.  See  publication  CRDEC  6. 

*  Limitations  in  the  use  of  Fokker-Planck  equations  for 
describing  growth  of  clusters,  as  used  in  development  of  classical 
Becker-Doring-Zeldovich  homogeneous  nucleation  theory,  have  been 
demonstrated  .  See  publication  CRDEC  7. 

*  Moment  methods  have  been  extended  to  include  aerosol 
evaporation  for  the  first  time.  This  permits  easy  and  rapid  numerical 
simulations  of  evaporation  processes.  See  publication  CRDEC  8. 

*  The  production  of  ultrafine  metallic  particles  by  evaporation 
of  metallic  samples  in  the  gas  phase  has  been  employed  in  research  and 
commercial  applications  for  at  least  50  years.  In  spite  of  this  long 
history,  the  processes  determining  particle  size  concentration  and 
morphology  have  not  been  understood.  Through  an  experimental  and 
theoretical  program,  we  have  developed  a  iv-ie.  which  explains  all  the 
important  features  of  particle  formation  and  growth  in  these  applica¬ 
tions.  See  publication  CRDEC  9. 


*  The  first  thorough  experimental  and  theoretical 
investigations  of  ferromagnetic  particle  formation  by  borohydride 
reduction  in  aqueous  phase  have  been  carried  out.  New  regimes  of 
particle  growth  have  been  discovered.  Characteristics  of  primary 
particles  have  been  explained  through  a  model  based  on  a  chemical 
mechanism  for  reduction  of  cations  by  borohydride  ions.  See  publication 
CRDEC  10. 

*  The  first  quantitative  experimental  data  have  been  obtained 
for  binary  aerosols  formed  by  homogeneous  nucleation  of  vapor  mixtures. 
The  experimental  results  are  explained  by  a  theoretical  model  for  binary 
homogeneous  nucleation  and  growth.  See  publications  CRDEC  11  and  12. 
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*  We  have  carried  out  the  first  quantitative  simulation  of  the 
explosion  of  water  droplets  caused  by  absorption  of  intense  laser 
radiation.  Substantial  departure  from  classical  blast  wave  theory  is 
shown  to  arise  from  homogeneous  vapor  nucleation  early  in  the  shock 
expansion.  See  publication  CRDEC  13- 
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*  We  have  carried  out  the  first  quantitative  simulation  of 
rapid  evaporation  of  small  droplets  exposed  to  intense  laser  beams. 
Computer  code  IRDE  (radial  isobaric  radial  droplet  evaporation)  has  been 


i 

I 


developed  for  this  simulation.  See  publication  CRDEC  14. 


*  The  first  simulation  of  the  evolution  of  the  aerosol  size 
distribution  in  the  presence  of  a  laser  beam  has  been  carried  out.  See 
publication  CRDEC  15. 
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*  The  first  complete  theoretical  investigation  of  fibrous  and 
membrane  filter  filtration  mechanisms  has  been  carried  out  numerically. 

See  publication  CRDEC  14. 
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III.  DETAILED  SUMMARY  OF  CRDEC  CONTRACT  WORK 


Complete  details  of  work  performed  under  CRDEC  Contract  DAAK- 
83-K-0006  are  given  in  the  publications  listed  in  Section  II.  In  this 
section,  a  detailed  summary  is  given  of  the  general  results  from  this 
research  in  the  four  areas  listed  in  Section  I.  First,  work  done  on 
particle  formation  by  homogeneous  nucleation  and  growth  is  discussed. 
Extension  of  this  work  to  atmospheric  dispersion  of  evolving  plumes  is 
presented  next.  Then,  work  done  on  nonlinear  interactions  of  intense 
laser  beams  with  aerosol  particles  is  summarized.  Finally,  first  results 
from  investigations  of  aerosol  filtration  from  the  standpoint  of  protec¬ 
tion  is  outlined. 

A.  Aerosol  Growth  Processes  and  Plume  Mechanics 

The  most  important  and  usual  means  for  production  of  particles 
less  than  the  order  of  1  m  involves  homogeneous  nucleation  in  a  vapor  or 
liquid  and  subsequent  growth  by  condensation/ evaporation  and  coagulation. 
Hiis  part  of  our  research  has  been  directed  toward  investigation  of  these 
processes  and  their  incorporation  in  models  for  atmospheric  dispersion  of 
aerosols.  In  this  section  work  on  homogeneous  nucleation  and  particle 
growth  is  discussed  first.  This  begins  with  work  on  nurleat.ion  and 
growth  of  liquid  aerosol  followed  by  work  on  metalli'  aero so’ s  and  ryq ro- 
sols.  Then  simulation  work  on  atmospheric  1  is  per  si  a  . s  presen  1 ed . 
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A-1.  Aerosol  formation  and  growth  in  a  laminar  coaxial  jet. 

In  nature  and  many  technical  applications,  aerosol  is  formed  by 
homogeneous  nucleation  of  a  vapor.  A  few  examples  include  condensation 
in  plumes  from  industrial  sources,  in  rocket  nozzles,  in  explosions  and 
in  various  natural  processes.  Subsequent  to  formation  by  nucleation,  the 
aerosol  may  grow  by  condensation/ evaporation  and  coagulation.  In  spite 
of  the  importance  of  these  processes,  there  are  few  available 
experimental  data  suitable  for  tests  of  theories  of  parcicle  formation 
and  growth.  A  number  of  studies  have  been  carried  out  in  turbulent  flows, 
e.g.  Amelin1,  Higuchi  and  O'Konski^,  Hidy  and  Fried  landed,  Delattre  and 
Friedlander1*,  Ismagilov  et  al.5,  Anisomov  et  al.^  and  Sutugin  et  al.?; 
however,  measurement  and  interpretation  of  particle  formation  and  growth 
in  turbulent  flows  pose  difficulties  which,  until  now,  have  not  been 
overcome  completely. 

The  principal  purpose  of  our  work  has  been  to  investigate  the 
dependence  of  the  disperse  properties  of  binary  aerosols  formed  by 
homogeneous  nucleation  on  the  physicochemical  properties  of  the  mixtures. 
We  have  carried  out  studies  by  following  the  condensation  of  oil  vapors 
in  a  laminar  coaxial  jet.  A  laminar  coaxial  jet  provides  a  fluid 
dynamical  state  which  can  be  modelled  numerically  with  considerable 
accuracy  and  also  eliminates  possible  influences  of  wall  effects  in  the 
initial  stages  of  aerosol  growth.  The  principal  theoretical 
uncertainties  in  the  analysis  of  such  a  flow  system  are  found  in  the 


description  of  homogeneous  nucleation  of  the  vapor  and,  less  importantly, 
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in  expressions  for  the  dynamical  coefficients  of  coagulation  and 
condensa  tion. 


a.  Experimental  procedure. 


A  schematic  diagram  of  the  experimental  system  is  shown  in 
Figure  A  —1  — 1 .  The  principal  elements  are  the  liquid  oil  and  carrier  gas 
metering  systems,  the  oil  vaporizer,  coaxial  jet  system,  and  the  sampling 
and  aerosol  measurement  systems. 
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Figure  A— 1  —1  .  Schematic  diagram  of  experimental  laminar  coaxial 
jet  system. 


Liquid  oil  was  delivered  at  accurate  volumetric  flow  rates  in 
the  range  0.001  -0.02  cm3  min-1  by  a  calibrated  syringe  pump.  Pure 
nitrogen  carrier  gas  was  metered  through  calibrated  rotameters.  The 
liquid  oil  was  vaporized  in  the  nitrogen  carrier  gas  in  a  heated  packed 
column.  Proper  design  of  this  vaporizer  was  essential  to  obtain  a 
steady,  laminar  coaxial  jet  from  the  nozzle.  Also  essential  was  accurate 
temperature  control  of  the  nozzle  heater  section;  this  was  achieved 
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through  use  of  a  DC  power  supply  controllable  to  +  0.01  W  at 


approximately  19  W  total  power  level  coupled  to  a  small  thermocouple 
sensor  attached  just  below  the  outer  exit  rim  of  the  nozzle.  In  the 
experimental  aerosol  growth  chamber,  superheated  oil  vapor,  of  a 
concentration  known  accurately  from  the  liquid  oil  metering  rate 
and  nitrogen  carrier  gas  flow  rate,  issued  as  a  laminar  coaxial  jet  from 
a  0.2  cm  diameter  circular  nozzle  (N  in  Figure  A —1  — 1 )  along  the  axis  of  a 
7.6  cm  diameter  glass  tube  (C  in  Figure  A  —  1  —  1  ) .  A  coaxial  flow  of 
filtered  air  at  25°C  was  maintained  at  determined  flow  rates.  With  this 
arrangement,  the  laminar  coaxial  aerosol  jet  can  be  caused  to  travel 
along  the  axis  of  the  larger  tube  for  its  entire  length  of  1.25  m. 
Experimental  and  theoretical  analysis  show  that  jet  temperature  and 
velocity  are  equilibrated  with  the  coaxial  air  flow  at  distances  greater 
than  several  hundred  nozzle  diameters. 

Aerosol  sampling  was  carried  out  with  a  cone  sampler  (S  in 
Figure  A  - 1  — 1 )  which  collected  the  total  coaxial  flow  and  produced 
(depending  on  the  coaxial  air  flow  rate)  fhom  10  to  25-fold  dilutions  of 
the  aerosol  in  times  of  the  order  of  a  few  milliseconds.  A  larger 
aerosol  jet  diameter  would  have  permitted  point  sampling  to  investigate 
the  radial  variation  of  aerosol  dispersity.  However,  it  was  found  that 
sufficiently  large  nozzles  posed  the  practical  problem  of  unwanted 
thermal  convective  motion,  and  we  were  not  able  to  overcome  this 
difficulty.  Hence,  only  the  average  aerosol  dispersity  over  the  jet 
cross-section  was  investigated.  Within  times  of  the  order  of  0.01  s  the 
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diluted  stream  from  the  cone  sampler  was  subsequently  diluted  by  filtered 


air  in  a  two-stage  isokinetic  dilution  system  by  dilution  factors  in  the 
range  1000-10,  000.  The  reason  for  this  large  dilution  is  given  below. 
The  exact  dilution  was  determined  by  metering  in  CO  as  a  tracer  gas,  as 
shown  in  Figure  A-1-1,  and  determining  the  CO  dilution  by  means  of  a 
calibrated  non-dispersive  CO  sensor.  The  dilutions  used  were  too  large 
to  permit  accurate  measurement  of  dilution  by  methods  described  in  the 
1 i terature . 

A  major  interest  in  this  study  was  the  disperse  properties  of 
aerosols  after  nucleation  and  condensation  growth  was  substantially 
complete;  the  sampling  technique  used  here  was  satisfactory  for  this 
purpose.  Order  of  magnitude  estimates  for  the  volatility  of  oils  used 
show  no  measurable  evaporation  of  the  aerosols  during  or  after  sampling, 
as  substantiated  by  the  experimental  measurements. 

The  principal  means  employed  for  measurement  of  particle  size 
distribution  was  a  Cliraet  208  optical  particle  counter  (0PC  in  Figure  A- 
1-1)  whose  output  signals  were  processed  by  a  Nuclear  Data  pulse  height 
analyzer  and  minicomputer  system  giving  512  channels  over  the  useful 
range  of  the  Climet  (0.3-5  m  particle  diameter).  It  was  found  to  be 
essential  that  calibration  curves  were  obtained  for  the  particular  oil 
aerosol  to  be  measured.  Determination  of  particle  size  distribution  from 
Climet  signals  was  found  to  depend  strongly  on  particle  refractive  :  nd  ex . 
Calibration  of  the  Climet  signals  was  obtained  by  generating  monod  ispe rse 
aero  sols  of  the  oils  by  means  of  a  parallel  plate  electrical  mobility 
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analyzer  whose  performance  was  further  checked  by  using  monodisperse 
polystyrene  latex  spheres  in  the  range  0.2-5  ;.m  with  diameters  determined 
prior  to  use  by  electron  microscopy.  In  using  the  Climet  208,  it  was 
also  found  to  be  essential  that  the  aerosol  sampled  was  diluted  to 
particle  number  concentrations  less  than  100  cm~3.  a t  higher  sampled 
number  concentrations,  spurious  particle  counts  are  obtained  which  give 
particle  concentrations  incorrect  by  factors  of  10  to  100.  A  TS I  3030 
electrical  aerosol  analyzer  was  also  used  in  studies  of  early  stages  of 
aerosol  growth  in  the  jet,  although  these  results  are  only  of  qualitative 
interest  and  are  not  presented  here.  An  electronic  cascade  impactor  (ECI 
in  Figure  A— 1  —1 )  was  used  in  a  limited  series  of  studies  to  measure  com¬ 
position  as  a  function  of  particle  size  in  binary  aerosol. 

The  description  of  the  laminar  jet  is  a  classical  problem  in 
fluid  dynamics  and  a  vast  literature  exists  on  the  subject  (e.g.  Pai  1  3, 
Birkhoff  and  Zarantonello14) .  The  laminar  coaxial  jet  employed  here  can 
be  completely  described  by  the  following  conservation  equations  of  mass, 
momentum,  and  energy: 

*  ■  v  =  0  (  A  — 1  -1 ) 

•  .  v  \  +.  p  +  -  *  :  g  =  0  (  A  -1  -2 ) 

•  °v  -  ’  'T  =  • '  k  T  :  • v  (  A-1  -3) 

where  is  the  fluid  mass  density,  v  the  mass  average  velocity,  p  the 


pressure,  the  shear  stress  tensor,  £  the  gravitational  acceleration,  cv 
the  heat  capacity,  T  the  fluid  temperature,  and  k  the  thermal 


conductivity.  The  vapor  and  aerosol  mass  concentrations  in  the  jet  are 


approximately  three  orders  of  magnitude  less  than  the  mass  concentration 
of  the  carrier  gas;  consequently,  the  vapor  and  aerosol  may  be  regarded 
as  passive  contaminants  and  do  not  influence  the  dynamics  of  the  coaxial 
jet  described  by  equations  (A— 1  -1)— (A— 1  —3). 

A  computer  program  has  been  developed  which  permits  solution  of 
equations  (  A-1  -1  )-( A-1  -3)  by  modeling  an  ax  isymmetr  ic  ,  two-dimensional, 
steady  laminar  flow.  The  numerical  calculations  are  carried  out  using  a 
finite  difference  method  and  a  variable  mesh  size  to  accommodate  the 
rapid  changes  taking  place  in  the  entrance  region.  Figures  A-1 -2,  A -1 -3 
and  A-1 -4  show,  respectively,  the  calculated  radial  velocity,  axial 
velocity  and  temperature  distributions  for  the  laminar  coaxial  jet  for  an 
initial  average  nozzle  velocity  of  94.4  cms*1  and  temperature  of  270°C 
with  an  average  velocity  in  the  annular  regionof  2.34  cms*1.  These 
initial  values  are  those  used  for  the  aerosol  measurements  reported  here. 

Such  calculations  permit  detailed  simulation  of  aerosol 
formation  and  growth  with  spatial  resolutions  in  the  velocity  and 
temperature  fields  not  possible  by  present  experimental  techniques.  As  a 
check  that  the  theoretical  model  approximated  the  experimental  system, 
average  temperatures  in  the  aerosol  jet  were  measured  by  means  of  small 
thermocouples.  Figure  A-1 -5  shows  the  approximate  agreement  between  the 
calculated  and  observed  average  jet  temperatures  as  a  function  of 
downstream  distance.  These  measurements  therefore  also 
computed  velocity  fields  in  the  jet. 
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confirm  the 


Calculated  radial  velocity  in  experimental  laminar  co¬ 
axial  jet  system.  VR  is  radial  velocity,  VNA  average 
velocity  at  nozzle  exit,  Z  is  downstream  distance,  R  is 
radial  distance,  and  RN  is  nozzle  radius  (0.1  cm). 


Figure  A-1-4.  Calculated  temperature  distribution  in  experimental 
laminar  coaxial  jet  system.  T  is  temperature,  TW  i3 
wall  temperature,  and  TN  nozzle  temperature.  For  other 
nomenclature  see  Figure  A  —1  —2  legend. 


Figure  A  —  1  — 5 .  Calculated  an  J  observed  average  jet  tern  per  •  •  .ire  a  ?  i 
function  of  ax  ia  1  downstream  distance  in  I  am  mar  o.;x  i  < 

ft  . 
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Result  3 


binary  aerosol  after  nueleation  and  condensation  growth  of  the  aerosol 
was  completed,  lit  the  particle  number  concentrations  produced  in  our 
experiments,  coagulation  was  relatively  unimportant  over  the  relevant 
time  scales;  our  measurements  were  concerned  with  growth  times  less  than 
1  s  whereas  characteristic  coagulation  times  for  the  aerosol  produced 
were  >  10  s,  even  during  homogeneous  nueleation,  according  to  our  model 
calculations.  As  pointed  out  above,  the  sampling  and  dilution  system 
used  was  adequate  for  accurate  sampling  of  these  aerosols. 

The  use  of  CO  tracer  permitted  accurate  determination  of  total 
dilution  of  the  aerosol  measured  by  the  OPC.  Particle  number 
distributions  were  then  calculated  using  the  measured  dilution  factor.  A 
confirmation  of  these  measurements  was  provided  by  calculating  the  mass 
distributions  from  which  the  total  mass  concentration  was  in  turn 
calculated.  For  those  oils  undergoing  total  condensation  to  aerosol, 
these  miss  concentrations  agreed  to  within  less  than  5$  with  the 
total  mass  concentrations  found  from  the  accurate  liquid  oil  and  nitrogen 
feed  rates. 

A  series  of  phthalate  esters  was  used  -  dibutyl  (DBP),  dihexyl 
1  CMP ' ,  and  ( i octyl p hthal ate  (POP)  -  along  with  squalene  (SQ)  (spinacene) 
and  oi*i:  acid  (.A),  and  their  mixtures.  For  all  but  dibutyl  phthalate, 
the  vapor  pressures  were  sufficiently  low  in  the  temperature  field  of  the 
jet  tn-i  t  all  condensation  growth  was  completed  at  axial  distances  greater 
than  approximately  ,?„  nozzle  diameters.  Measurements  at  downstream 
position:-  stowed  :hange  in  partial**  s  i  •*-  distributions,  except  for 


Dibutyl  phthalate  has  a  sufficiently  high  vapor  pressure  that 


condensation  growth  was  not  complete  at  the  usual  sampling  point  at  an 
axial  distance  of  20  cm  from  the  nozzle.  Figure  A  —1  —  6  shows  the  evolu¬ 
tion  of  the  particle  number  and  mass  distributions  at  axial  downstream 
distances  of  47.5,  70  and  100  nozzle  diameters  for  a  liquid  oil  feed  rate 
of  D3?  of  0,001  36  cm^min-1  and  178  cm^min-1  (298k,  1  atm.,  nitrogen 
carrier  gas).  The  "mixing"  ratio  defined  here  is  the  liquid  oil  feed 
rate  divided  by  the  nitrogen  carrier  gas  feed  rate.  The  mixing  ratio  in 
this  example  is  7.64  x  10’^.  Figure  A-1-7  shows  the  percentage  of  total 
mass  concentration  of  DBP  vapor  condensed  as  aerosol  as  a  function  of 
downstream  distance.  Total  number  concentration  is  constant  for  the 
three  positions  at  4.6  x  10^cm"^,  as  is  consistent  with  growth  of  an 
aerosol  solely  by  condensation. 

Dihexyl  (DHP)  and  dioctyl  (D0P)  phthalate  have  sufficiently  low 
vapor  pressures  in  the  temperature  field  of  the  jet  that  condensation  is 
complete  at  axial  distances  greater  than  approximately  20  nozzle 
diameters.  An  important  aspect  of  particle  formation  and  growth  is  the 
dependence  of  particle  size  distribution  and  its  moments  on  the  initial 
vapor  "'’■ncentration.  Figures  A-1 -8  and  A -1-9  show  particle  number  and 
mass  distributions  for  DOP  and  DHP  respectively  at  100  nozzle  diameters 
(20  cm)  downstream.  In  these  measurements,  nozzle  temperature  was  270°C, 
nitrogen  feed  rate  was  178  cm3  min’1  at  293'K,  1  atm.,  an  1  liquid  ail  feed 

rate  0  00206,  0.0034  and  0.0051  ern^min’1  (mixing  ratios  of  1 . 1 6  x  10’^, 
I.--’  x  10’J  3nd  2.87  x  10”6).  in  all  cases ,  total  mass  concentration 
’tl  .jilted  from  the  measured  particle  size  d ;  str  ihut  ior.r  igr«-.»d  to  within 
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character  of  the  results  obtained  from  our  experimental  system.  .A 
complete  presentation  of  binary  aerosol  data  with  a  more  complete 
analysis  will  be  given  in  additional  papers.  In  the  next  section  partial 
analysis  of  some  of  the  results  presented  will  be  given. 


c.  Analysis  of  results. 


A  complete  analysis  of  experimental  data  obtained  from  the 
laminar  coaxial  aerosol  jet  would  involve  numerical  solution  of  the 
spatially  inhomogeneous  dynamical  equations  for  aerosol  growth  (e.g., 
Brock^)  coupled  with  the  solution  of  equations  (A— 1  — 1  )  —  CA— 1  —3  )  •  An 
analysis  of  this  stystem  considering  only  coagulation  of  an  aerosol  was 
published  earlier  (Hajuiohi  and  Brock).  Here  an  approximate  analysis 
will  be  given  by  consiiering  an  average  over  the  aerosol  jet  at  a  given 
downstream  distance.  kit'  justifiable  neglect  of  axial  and  radial 
diffusion  of  particles,  the  average  properties  of  the  aerosol  jet  at  a 
downstream  position  will  be  modeled  by  obtaining,  from  the  numerical 
solution  of  equations  •  A -1  !  -  ( A- 1  -3  ) ,  the  average  axial  transport  time 
and  temperature. 

A  dynamical  model  for  evolution  of  the  particle  s i  z e 
distribution  can  be  developed  along  these  lines,  but  such  a  program  will 


be  carried  out  elsewhere. 


Here  we  present  results  based  on  an 


approximate  moment  formulation  of  the  spatially  averaged  (at  a  downstream 


:ross-sec  tier. )  single  com  ponen  t  part  i  cl  • 
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being  the  axial  downstream  distance)  is  n(x,t)dx.  We  study  the  evolution 


in  this  time  through  the  moments  of  n(x,t): 


M0  = 

n(  x  ,t )  dx 

(  A-1  -9  ) 

0 

Ml  = 

/  xn(x , t) dx 

(A-1  -5  ) 

m2  = 

co  ^ 

f  x~n(x,t)dx 

Q 

' A-1 -6  ) 

C  +  M 

1  =  constant. 

'A-1  -  ’  : 

where  C  is  the  vapor  mass  concentration.  For  a  single  component  aerosol, 
the  evolution  equations  are  obtained  from  the  equation  for  n'x,t' 
appropriate  for  the  experimental  conditions  studied: 

■n/  <t  =  -( V  ~ix)  (vn)  +  1/2  (.2/.x2)(,n)  +  R  (A-1 -8) 


where  vis  the  particle  growth  rate,  a  is  the  diffusion  coefficient  of 
particle  growth,  and  R  the  nucleation  rate.  Neglecting  the  diffusion 
term  and  evaporation,  one  obtains  from  ( 4 )—  (6  )  and  (8): 

dM0/dt  =  R  (A-1-9) 

dM./dt  r  /  ’-ndx  +  x„R  ( A— 1  — 1 0 ) 

'  o  ^ 


Also  : 


dt^/dt  =  2  /  x Vndx  +  xcR 


(  A -1  -11  ) 


d(M1+C)/dt  =  0 


(  A-1  -12) 


where  the  nucleation  process  is  assumed  to  produce  a  single  time 
dependent  nucleus  size,  x  (t): 


R  =  l'(x-x  (t))R 

c 

where  c  is  the  delta  function  and  R  is  the  nucleation  rate.  R  is 
calculated  according  to  the  classical  theory  of  Becker-Doring-Zeldov  ich 
( e  ,g .  Abraham1 ) . 

The  particle  growth  rate  has  the  form  (Davis  et  al.?): 

-  -1 

;  =  (  rrD 2 / 4 )  t'  C  (S-eKe)n1/2Kn 

y  p  1  v 

where  is  particle  diameter,  mean  thermal  speed  of  vapor  molecules, 
C  the  concentration  corresponding  to  the  bulk  vapor  pressure,  S  the 

V 

supersaturation  ratio,  Ke  the  Kelvin  and  Kn  the  Knudsen  numbers. 
Coagulation  is  unimportant  for  the  conditions  of  this  study. 

Equations  (A  — 1—9)— (A— 1  — 12)  are  not  a  closed  set  for  the 
moments.  Closure  is  only  possible  in  the  free  molecule  regime  if  ip  does 
not  contain  the  Kelvin  evaporation  term;  this  situation  does  not  exist 
for  our  experimental  system. 

Various  lumped  parameter  and  mean  growth  models  of  nucleation 
and  growth  have  been  proposed  (e.g.  Takahashi 1 5).  However,  we  have  found 
that  such  models  cannot  explain  our  experimental  observations.  Closure 
of  equations  (A— 1  —4  )— (A— 1  —12)  is  effected  here  by  introducing  an  assumed 
particle  size  distribution.  This  procedure  is  similar  to  methods  pro¬ 
posed  previously  for  coagulation  and  condensation  (e.g.  Clark^).  We  have 
used  the  lognormal  distribution  in  these  simulations,  although  other 
distributions  may  be  superior  for  this  purpose.  For  simple  condensation 
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With  these  these  definitions,  equations  ( A -1  -9  )-( A-1  -1 2)  are 


closed  and  may  be  solved  in  integral  form.  Such  solutions  provide  at 
least  a  qualitative  explanation  of  our  experimental  results.  This  model 
has  also  been  extended  to  multicomponent  aerosols. 

Figure  A  —  1  —  1 0  shows  the  calculated  nucleation  rates  and 
particle  number  concentration  as  a  function  of  t(z)  for  single  component 
aerosols  of  DOP  and  DHP  at  a  mixing  ratio  of  2.1  x  10“^  in  the  laminar 
coaxial  jet  from  equations  ( A-1  -9 )-( A-1  -1 5).  Figure  A-1-11  displays  the 
correspond ing  geometric  mass  mean  particle  diameters  and  standard  devia¬ 
tions.  It  can  be  seen  that  the  nucleation  process  is  completed  in  times 
of  the  order  of  a  few  milliseconds. 

The  characteristics  of  the  two  aerosols  produced  in  the 
com  pari  son  of  Figures  A  —  1  —10  and  A-1-11  are  the  result  of  competitive 
rate  processes:  nucleation  and  condensation.  Upon  initiation  of 

nucleation,  all  of  the  available  vapor  goes  to  the  creation  of  particles. 
As  the  particle  concentration  increases,  the  vapor  begins  to  condense  on 
existing  particles  as  well  as  to  contribute  to  creation  of  new  particles. 
The  condensation  rate  increases  and  eventually  condensation  depletes  the 
vapor  bringing  about  a  rapid  decrease  in  new  particle  formation  by 
nucleation. 

The  cooling  rate  of  vapor  in  the  jet  plays  a  key  role  in 
determining  the  maximum  nucleation  rate  and  hence  the  number  of  particles 
produced.  Maximum  cooling  rates  in  our  experiments  were  of  the  order  of 
1014  Ks*1.  Large  cooling  rates  permit  very  rapid  increases  in  nucleation 


Figure  A-1 -1 1 .  Calculated  (a)  mass  average  particle  diameters  and  (b) 
standard  deviations  for  dioctyl  phthalate  (DOP)  and 
dihexyl  phthalate  (DHP)  condensing  in  experimental 
laminar  coaxial  jet  for  conditions  described  in  Figures 
A-1 -3  through  A-1 -5.  Mixing  ratio  =  21  E-5.  Key:  DOP 
( _ .),  DHP  ( - ). 

rates  while  at  the  same  time  the  total  condensation  rates  increase  much 
less  rapidly.  In  general,  larger  cooling  rates  produce  higher  particle 
number  concentrations  and,  for  a  fixed  initial  vapor  concentration, 
smaller  ultimate  mean  particle  sizes.  Turbulent  jets  create  much  larger 
cooling  rates  than  those  found  in  laminar  jets;  consequently,  for  the 
same  initial  vapor  concentrations  and  other  conditions,  turbulent  jets 
produce  much  higher  particle  number  concentrations  and  smaller  mean 
particle  sizes  than  are  found  here  for  a  laminar  jet. 

In  the  laminar  coaxial  jet,  particle  growth  continues  long 
after  cessation  of  nucleation.  The  plots  of  mean  particle  diameter  and 
standard  deviation  shown  in  Figure  A-1 -1 1  indicate  the  apparent 
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complexity  of  the  growth  process,  although  the  theory  used  here  Is 


admittedly  approximate.  The  shoulders  in  the  curves  for  particle 
diameter  may  be  due  to  the  existing  large  nucleation  rates  and  consequent 
rapid  introduction  of  smaller  particles,  which,  however,  rapidly  increase 
in  size  at  the  high  supe r sa tur a t ion  prevailing  at  that  time.  At  the 
times  of  these  shoulders,  the  standard  deviations  are  seen  to  exhibit 
similar  changes.  The  standard  deviations  are  seen  to  decrease  during 
late  stages  of  particle  growth.  This  agrees  with  the  well-known  fact 
that  condensation  growth  of  an  aerosol  produces  a  narrowing  of  its 
distribution . 

In  this  study,  as  well  as  in  previous  studies  of  nucleation  and 
particle  growth,  it  is  found  that,  other  conditions  being  the  same,  less 
volatile  species  yield  aerosols  with  higher  number  concentrations  and 
smaller  mean  particle  sizes.  Figure  A  —1  — 1 2  shows  a  comparison  of 


calculated  and  observed  total  particle  number  concentrations  for  DOP  and 
DHP  as  a  function  of  mixing  ratio.  As  is  evident,  the  qualitative  trend 
is  reproduced  by  the  model,  although  there  is  a  quantitative  difference. 
This  difference  may  arise  from  several  sources.  Classical  nucleation 
theory  is  used  in  the  calculation,  although  the  classical  theory  is  known 
to  diverge  widely  from  experiment  at  the  high  nucleation  rates  believed 


to  exist  in  our  experiments.  A  decrease  of  10^  in  the  nucleation  rate 
calculated  by  the  classical  theory  for  DHP  would  bring  about  a  much 
better  quantitative  comparison.  Also,  the  theory  applied  here  is 
approximate  and  conditions  are  assumed  constant  over  a  cross-section  of 
the  vapor-aerosol  jet. 
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Figure  A  —1  —12.  Comparison  of  calculated  and  observed  particle 
concentrations,  with  all  vapor  condensed  to  a ercs< 
dioctyl  phthalate  and  dihexyl  phthalate  acres 
experimental  laminar  coaxial  jet  as  function  of 
ratio.  Experimental  conditions  same  as  desrrit 
Figures  A— 1  —3  through  A— 1  —5.  Key:  DCP  oalcu;  at-*  : 
DOP  observed  ( - ),  DHP  calculated  ( - Drip 


The  extension  to  binary  aerosols  follows  the  procedure  ■ 
above.  The  binary  density  function  n^(x^,  x~>,  t.1  is  define.: 

x2>  dx-|  dxo  is  the  number  of  particles  having  masses , 
species  1  and  x^  of  species  2  in  the  range  dxl7  dx,  about  x,  and  x 
Moments  may  be  defined: 


Mij  =  ■ x i  +  xj5  "V^T-Xp.M  dx,  !x- 

The  evolution  of  t  ne>  moments  :r.  t.m>-  :  .<■  •  ; :  : 

the  evolution  equati  on  fo  r  r. .  ^  .  , 

en  tal  oond i  tions  : 


Hi))  -  —  --  ?  rin)  +  R  i  - 


‘ .  i  n  i 


The  moment  equations  studied  have  the  form: 


dM  02 


■  1 n  1 2 ( X1  •  x2 »  t)  dxi  Ox.,  +  x  K 


.2n12(x1'  x2 ’  tJdXTdx,  +  x.,  K 


;  ;  x1vin]2(x1,  x2,  t)  dXjdx,  +  x,  R12 


2  r  r  xi^2n12(x1»  x2 *  t) dx i  dx2  +  x2  R]2 


■  1  and  .p  are  the  particle  growth  rates  for  species  1  and  2.  Details 
may  be  found  in  the  references  cited  in  Section  II.  R1  ^  is  the  classical 
binary  nucleation  rate,  a  st  r  a  ight  fo  r  ward ,  but,  unfortunately,  limited 
extension  of  the  single  component  classical  nucleation  theory  cited 


above . 


Closure  of  the  moment  equation  above  is  obtained  by  assuming  a 


lognormal  product  for 


t< 


)  • 


n  1 2 ^  x  1  »  x5*  ^  =  (Mq/2  In In  ,  ?  Xy  x2 ) 

exp  C - ( lnx  1  _  lnxy)2/2lnr' 
exp  (-(lnx?  -  lnx2  )^/21n^ .  ,  ) 

a  consistent  density  function  for  binary  internally  mixed  aerosols  as  in 
the  present  study. 

Detailed  presentation  of  results  is  too  lengthy  to  repeat  here. 
Some  limited  comparisons  are  given  her j  of  predicted  and  experimental 
values  of  the  total  number  concentration,  Mqq  as  shown  in  Figures  A —1  —1  3 
and  A-1-14  for  DHP  in  DOP  and  DHP  in  SQ  (squalane).  In  view  of  the  well- 
known  inadequacies  of  current  binary  classical  nucleation  theory,  the 
agreement  exnibited  in  these  figures  is  surprisingly  good.  It  should  be 
noted  that  these  are  independent  comparisons  with  no  adjustable 
parame  ters . 


Figure  A —1  —1  3 .  Asymptotic  number  concentration  as  a  function  of  volume  % 
dihexyl  phthalate  (DHP)  in  dioctyl  phthalate  (DOP). 
.  theory;  points  with  error  bars,  experimental 


Z-3L  X  (E..UJ3/#) 


Ac(x)  =  2 "Dp PCV [ 1  +  Dp  /i/2  M0<Dp>  1, 


where  Dp  is  the  particle  diameter:  x  =  (1/6)"Dp^d;  V  ,  the  diffusion 
coefficient  of  monomer  vapor  in  host  gas;  Mq,  the  particle  number  concen¬ 
tration;  <D  >,  the  mean  particle  diameter;  and  Cv,  the  equilibrium  mass 
concentration  of  monomer  vapor  for  a  plane  interface.  The  terra  m 
brackets  in  this  expression  accounts  for  competitive  effects.  For 
particle  sizes  restricted  to  those  appropriate  for  aerosols  and  for 
sufficiently  small  number  concentrations,  e  D  / 1  /2'Mq<D  >  <<  1  (and 
apparently  competitive  effects  can  be  neglected).  Alternatively,  the 

equilibirium  volume  fraction  (Marqusee  and  Ross2?),  a  (C-  •  .  ,_c  )d 

in  i  l  13  i  v 

fulfills  the  condition  •  <<  1.  For  situations  where  e  is  not  small  for 
aerosols,  competitive  effects  become  important;  however,  this  corresponds 
to  larger  values  of  Mq,  for  which  Equation  ( A  —2  — 1 )  is  no  longer  the 
correct  evolution  equation,  and  the  coagulation  process  must  be  included. 
The  importance  of  competitive  growth  for  aerosols  when  coagulation  is 
present  has  not  yet  been  examined  completely,  although  this  is  discussed 
next  in  Section  A-3. 

In  any  event,  our  simulations  are  restricted  to  the  regime 
where  competitive  effects  may  be  regarded  as  negligible.  For  continuum 
growth,  therefore: 


Ac(x)  =  2”DpCCv 

For  kinetic  (free  molecule)  growth: 

H  -  i 

V*)  =  1/14  3P  '•  — 

mi 


1  /2 


,?/? 


•k" 


(A-2-3) 


;  a  -2-i 


Here,  mi  is  the  monomer  molecular  mass  ami  •  T  the 


A  .V  A  .* 


thermal  energy. 


Equation  (A-2-1)  is  coupled  to  the  conservation  equation  for 


monomer : 


-/  v<  x,t )  n(  x,t )  dx  -  x*  T  ( x»  ,t  )n(  x*  ,t )  , 
x* 


(  A -2  -5 ) 


x*  is  the  mass  of  the  smallest  particle  that  obeys  the  growth  law, 
Equation  (A-2-2). 

Two  important  moments  of  n(x,t)  are  the  total  particle  number 
concentration,  Mq,  and  mass  concentration,  M  i  .  From  (1)  and  (5): 

dM0 

-  =  v  (x*  ,t)n(  x*  ,t)  .  (A-2-6) 

dt 

d 

—  (Cvs  +  M, )  =  0.  (A-2-7) 

dt 


Equation  (A-2-7)  reflects  the  conservation  of  total  mass  concentration. 

Equations  (A-2-1)  and  (A-2-5)  constitute  a  nonlinear  integro- 
differential  equation  for  n(x,t).  The  asymptotic  properties  of  (A-2-1) 
and  (A-2-5)  with  the  linearization: 

exp  (K/x1/3)  1  +  K/x1/3  (A-2-8) 

and  x*  =  0,  were  apparently  first  investigated  correctly  by  Lipschitz  and 
Slyozov^  (LS);  subsequent  investigations  have  not  modified  tne  LS 
results.  Asymptotic  solutions  are  obtained  in  terras  of  the  similarity 
var  lab  i.  es  : 

'  =  3  in  (Dp/Dpc(t))  , 

=  Dp/Dpp » 

is  implicitly  related  to  t  through  the  dependence,  DpC(t).  D^c  is  the 
critical  particle  diameter  separating  the  regions  of  evaporation  and 


F rom  '  2 ) : 


Dpc  =  4  ./-Tins  (A -2 -9) 

or,  with  the  linear  approximation.  Equation  (A -2 -8): 

Dpc  '  •:Y/<T(s-1  ) . 

The  distribution  is  transformed  with  the  assumption  that  asymptotically 
(t-->  ■): 

r.(x,t)  =  (2/  d)Mg(  t)p(p  )/DpcDp.  (A-2-10) 

For  the  continuum  growth  law:  A(x)  3 ,  and  the  LS  result 


is 

33e 

p;.  ’  =  -  exp  [-1  /(  1-2  .73)] /(.■ +3  ) 7/3 ( 3/2- .  ) 1 1  73  (A-2-11) 

■P  i  3 


and  Mg  ?/t,  valid  for  the  linearization.  Equation  ( A  —2  —8 ) . 

For  kinetic  {  free  molecule)  growth:  A(x)  x^/,3>  an(j  it  is 
easy  tc  show  that 

?(.  )  =  24,exp[-3/(2-.  )]/(2-.  )5  (A-2-12) 

and  MV't"  t  “  '■  ,  valid  for  the  linearization,  Equation  ( A— 2  — 8 ) . 

.vs  rrt  here  results  of  numerical  investigations  of  Equations 

(A-  -’.-'A-'-  w  1  *  h  a;.:;  without  the  linearization.  Equation  (A-2-8),  for 

botn  continuum  f.  :  kinetic  'free  molecule)  growth  laws.  In  addition,  the 
important  effect  an  approach  to  the  asymptotic  limit,  Equation  (A-2-10), 

of  the  u :  : .  ;n  af  t.ne  initial  size  distribution  is  noted.  The 

.nume-'  1  -  a’  ‘  -  :  -.r. ;  ley*  :  for  solution  of  (A-2-1)  and  (A-2-5)  is  a 

com:.:'  r  r.,:  i  «n  method  that  we  have  described  elsewhere 


reported  a  numerical  solution  of  (A-2-1) 


~ qua t ion  A-2-8),  for  the  continuum 


m 


growth  law  with  ~ ;  —  o  i  a  1  restrict  i-ans  on  ini*  : 


introduce  mrv*  >f 


in  our  .numerical  simul  at  Lons,  the  in  l  t  i  .  :  -.l :  s  t.  r  : 


defined  in  *  era’s 


' L  st r  ibu t ; on  w  i  t r 


r.d  hav-;  the  far:;,: 


:’fPp.t  =  0 


P  .  ,  in  .  r.  1 

-  1 J'-p»  p*  & ' 1 


where  Dp  and  •;.c  are,  respectively,  the  mean  an  1  gn-ometr 
deviation  of  a  lognormal  distribution.  The  lognormal  list r 
been  suggested  (Jranquist  and  Bun r man  1  '}  to  .J  as  crib*  pc 

distributions  from  various  homogeneous  nucleating  processes. 

Table  A-2-1  summarizes  the  various  cases  studied  in 
tions.  Cases  1-6  cover  the  continuum  growth  law,  A'x,.  x1 
indicated  variations  in  the  Kelvin  term  and  initial  particle 
bution.  The  kinetic  growth  law,  aJx)  x-  '  3,  is  applied  i 


where  principally  the  effect  of  linearization  of  the  Kelvin 
existence  of  the  asymptotic  limit,  equation  (A were 
all  cases,  the  value  of  Cys  -  M,  is  the  same. 


The  questions  we  wish  to  examine  in  this  study  a 
are  the  dynamics  of  the  approach  to  the  asymptotic  1  i  m  i t :  (.I 
role  of  initial  particle  size  dispersion  on  approach  to  tr. 
limit  distribution?  ■})  What  is  the  effect  of  usinr  t  -■  .’err¬ 


or  rr.  for 


ipprox  it  •.  1 1  o r. ,  u  ; :  •  *  L 


TABLE  A -2-1  :  Cases  Studied  in  Simulation  of  Qstwald  Ripening 


Case 

S 

Mi  (  gm/ cm3) a 

Mq  w//cm3)a 

'  g 

n 

(um) 

Ke  lv  in 
effect 

g  ro  vt  h 
law3 

1 

1.1 

1.0  x  10-7 

2.25594x107 

1.4 

0.2 

1  inear 

Eq.  (6) 

cor.f  inuum 

2 

1.1 

1.0  x  10*7 

2.255  94x107 

1.4 

0.2 

nonlinear 
Eq.  (2) 

continuum 

3 

1.1 

1.0  x  10-7 

1.30565X107 

3-0 

0.2 

linear 

Eq.  (6) 

cont  inuum 

14 

1.1 

1.0  x  IQ-7 

1 .305  65  xl  07 

3-0 

0.2 

nonl inear 
Eq.  (2) 

con  t  inuu  m 

5 

1.10 

1.0  x  10-7 

1.87755x107 

2.0 

0.2 

nonl inear 
Eq.  (2) 

continuum 

6 

1.12 

2.31  37x10-7 

1 .65999x1 07 

5.0 

0.09 

nonl inear 
Eq.  (2) 

cont inuu  m 

7 

1  . 1 

1.0  x  IQ-7 

2.2559X107 

1.4 

0.2 

nonl inear 
Eq.  (2) 

kinetic 

8 

i .  i 

1.0  x  IQ-7 

2.2559x107 

1.4 

0.2 

linear 

Eq  .  ( 6  ) 

kinetic 

a 


The  numerical  values  used  for  parameters  in  the  study: 

Cv  =  3»3UE-6g/cn)3;  t  =  U2  3K ;  =20  dyn/cm;  d  =  1.0  g/cm^;  m-,  =6.43E-22g. 


b.  Results. 


Figures  A-2-1  -  A-2-3  show  the  approach  to  the  asymptotic  limit 
distribution,  Equation  (A-2-1)  with  continuum  growth,  A(x)  x1//3j  for 
cases  1-4  of  Table  A-2-1.  Figures  A-2-1  and  A-2-2  differ  only  in  that 
the  linearization.  Equation  (A-2-8),  is  employed  in  obtaining  the  results 
of  Figure  A-2-1  ,  while  the  full  nonlinear  term  is  used  in  Figure  A-2-2. 
In  both  cases,  the  asymptotic  limit  distribution.  Equation  (A-2-11),  is 
found  for  nondimensional  times,  t  -  15  within  the  limits  of  accuracy  of 
our  numerical  simulations.  In  Figure  A-2-3,  cases  3  and  4  are  presented 
for  relatively  large  initial  particle  size  dispersion,  -  =  3.  As  is 

evident,  the  asymptotic  limit,  Equation  (A-2-12),  is  not  achieved  at 
I  -  15,  as  was  found  for  cases  1  and  2  where  =  1.4.  Our  estimates 

indicate  that  much  larger  values  of  (  ~  100)  would  be  necessary  to 
achieve  the  asymptotic  limit  distribution,  Equation  (A-2-12). 

The  positive  moments  Mq,  H,,  M2,  M3,  M ^ ,  M^,  and  as  a 
function  of  simulation  time  are  shown  in  Figures  A-2-4  through  A-2-7, 
where 


Mj  =  xin(  x,t ) dx  . 

x  x* 

The  lower  bound  corresponds  to  a  particle  diameter  of  0.001  ,.ra. 

In  Figure  A-2-4,  Mq  and  M-|  are  shown  for  the  continuum  and 
kinetic  growth  laws  for  all  cases  in  Table  A-2-1.  It  is  easily  seen  that 
there  is  an  initial  "adjustment  phase"  during  which  Mq  remains  constant. 
This  reflects  the  finite  time  for  evaporating  particles  to  reach  x*  so 

.hat  Mq  p3n  begin  to  decrease  according  to  Equation  (A -7-*'  .  For  cur 

assumed  initial  particle  size  distributions,  n(  x#  ,t  =0  )  ;  •  ■  ■  :■  •  0  ‘  i  a  1 1  y 


-  _  -  t=0 .01  sec  (;  =1  .828) 

_  -  _  t=0.05  sec  (: =3*768) 

o-o-o-o  t=0.4  sec  (~  =4.913) 

_ _ Asymptotic  similarity  solution,  Equation  ( A  —2  —  11) 
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Figure  A-2-1.  Approach  to  asymptotic  limit  distribution  for  continuum 
growth  law,  case  1,  Table  A-2-1. 


- t  =  0  (1  =0) 

—  - _ -  t=0.2  sec  (t=2.35) 

-  - - t=0.05  sec  (t  =3-65) 

—  o-o-o-o  t=0.40  sec  ("=4.72) 

-  Asymptotic  similarity  solution,  Equation  (A-2-11) 
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*  a.-.v. met  otic  limit  distribution  for  continuum 


Case  3: 


Figure  A-2-3.  Approach  to  asymptotic  limit  distribution  for  continuum 
growth  law,  cases  3  and  4,  Table  A-2-1. 


oooooooooo  case  1  (continuum  growth) 

-  case  2  (continuum  growth) 

-  case  4  (continuum  growth) 

case  5  (continuum  growth) 

!  !  !  !  I  i  case  6  (continuum  growth) 

—1 _  __  cases  7,  8  (kinetic  growth) 
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case  2,  g=1,4 

case  4,  =3.0 
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Figure  A-2-5. 


5  10- 16 
\  10' 17 


icr2  to'*  io°  10’  io2  103  104  105 

Time  (sec ) 

Second  moment,  as  a  function  of  time  for  continuum 

growth. 


zero  and  remains  so  until  the  evaporation  "wave  front"  reaches  x*.  Cases 
1-6  and  7-8  all  have  the  same  value  of  (Cvs+Mi)*  Equation  (A-2-3),  so 
that  the  LS  theory  predicts  a  universal  behavior  as  t  -->  00  and  for 
continuum  growth,  Mq  ~  1/t.  For  kinetic  growth,  Mq  ~  asymptotical¬ 

ly.  This  is  confirmed  in  Figure  A-2-4. 

f"01"  g—  1  •  4 ,  for  continuum  growth,  Mq  1/t  for  simulation  times 

t>1  sec  (->15).  Much  longer  times  are  required  for  >:,  =  3.0  and  £_  =  5.0,  as 

6  © 

evident  in  Figure  A-2-4.  As  is  also  evident  from  Figure  A-2-4,  M1 
achieves  a  universal  quasi sta tionary  state  for  cases  1-6,  and  does  this 
more  rapidly  than  M^  assumes  its  asymptotic  state. 

According  to  the  LS  theory  for  the  linearized  Kelvin  term, 
equation  7  A  -7  -8  . ,  the  moments  M|,  are  easily  shown  to  have  the  following 
time  dependence: 


(  A -2-1  4) 


-  t*(  i_1 )  ,  i  =  0,  +1,  +2,  .... 

where 

=  1  continuum  grwoth; 
a  -  3/2  kinetic  (free  molecule)  growth. 

Some  of  these  positive  moments  are  shown  in  Figures  A -2 -5  through  A-2-7 
for  continuum  growth,  and  are  found  numerically  to  follow  the  relation, 
Equation  (A-2-14),  for  long  times.  It  should  be  noted  that  M.  obeys  (A- 
2-14)  strictly  when  but  a  single  particle  remains  in  the  system,  a  state 
that  can  never  be  obtained  in  an  infinite  system.  As  expected,  the 
larger  the  initial  dispersion  in  particle  size,  the  longer  the  time 
necessary  for  the  moments  to  achieve  the  asymptotic  state. 

- case  2,  >  -  =  1  •  4 

- case  4,  *;*-  3.0 
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Figure  b-2-6.  Third  and  fourth  moments,  M3(t)  and  Mu(t),  as  a  function 
of  time  for  continuum  growth. 
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- case  2,  eg=1,1< 

- case  4,  £g=3.0 


I0‘*  10°  10*  to* 
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Figure  A-2-7.  Fifth  and  sixth  moments,  M^(t)  and  Mg(t),  as  a  function 
of  time  for  continuum  growth. 

For  continuum  growth,  the  effect  of  the  linearization.  Equation 
(A-2-8),  on  the  approach  to  the  limit.  Equation  (A-2-11),  is  seen  in 
Figures  A-2-8  and  A-2-9  where  the  distribution  function  3xn(x,t)  is 
plotted  as  a  function  of  particle  diameter.  In  these  figures,  for  t  “  1 
sec,  Dpc  ~  <Dp> ,  the  mean  particle  diameter.  Therefore  the  region  to  the 
left  of  the  mode  in  each  distribution  represents  evaporating  particles. 
The  nonlinear  Kelvin  term  in  Equation  (A -2-2)  at  x*  is  much  greater  than 
the  linearized  term.  Equation  (A-2-8):  e1®  >>  11.  This  causes  the 

distribution  for  the  nonlinear  term  for  long  times  to  be  smaller  than  the 
corresponding  distribution  with  the  linear  approximation.  Equation  (A-2- 
8).  As  can  be  seen  from  Figure  A-2-8  for  continuum  growth,  this  is  not 
true  for  relatively  small  times. 
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Clearly,  for  particle  diameters  less  than  the  order  of  the  Kelvin 
diameter,  the  linearization,  Equation  (A -2-3),  is  not  at  all  an  accept¬ 
able  approximation.  Therefore,  the  basis  for  the  development  of  the  LS 
theory,  Equation  ( A  —2  —8 ) ,  is  physically  unrealistic,  particularly  in  view 
of  the  fact  that  particles  disappear  at  some  very  small  particle  size, 
where  Equation  (A-2-6)  is  incorrect.  However,  we  find  the  similarity 
solutions,  based  on  the  linearization,  Equation  (A-2-8),  to  yield  good 
approximations  to  the  correct  nonlinear  development  according  to  Equation 
(A-2-2).  The  explanation  for  this  is  as  follows. 

We  find  for  the  nonlinear  and  linear  cases  for  both  kinetic  and 
continuum  growth  laws  for  long  times  that 


dM0  dM0 

-  nonlinear  —  linear. 

dt  dt 


case  1,  Linearized  Kelvin  term,  Equation  (A-2-8) 
case  2,  nonlinear  Kelvin  term,  Equation  (A-2-2) 

iO’  | - 1 - 1 - 1 - 1 - 1 


(A-2-15) 


Figure  A-2-8. 
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Effect  of  linearization  of  Kelvin  term  on  evolut  ;  or 
number  distribution  function  3xn(  x  ,t)  for  co.ot  in 
growth  law  for  cases  1  and  .?,  Table  a -2-1 . 
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case  3,  linearized  Kelvin  term,  Equation  (A-2-8) 
case  4,  nonlinear  Kelvin  term,  Eauation  ( A-2-2 ) 


Figure  A-2-9- 
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Effect  of  linearization  of  Kelvin  term  on  evolution  of 
number  distribution  function  3xn(x,t)  for  continuum 
growth  law  for  cases  3  and  4,  Table  A-2-1. 


That  is,  the  rate  of  disappearance  of  particles  is  approximately  the  same 
for  nonlinear  and  linear  Kelvin  terms.  From  Equation  (A-2-6)  it  follows 


[  „  (  x*  ,t)  n(  x*  ,t)  ]  non^  inear  ~  [’+  (  x*  ,t)  n(  x*  ,t)  ]  nnear  *  (A-2-16) 

The  re  fore,  near  to  x*  ,  the  LS  distribution,  Equation  (A -2 -11),  and  the 
asymptotic  result,  Equation  (A-2-12),  cannot  be  correct.  However,  for  t 
—  >  ■,  the  large  differences  between  n(  x»  ,t)  nonl inear  and  n(  x»  ,t)  x  inear 

become  confined  to  a  "boundary  layer"  near  x* ,  and  the  asymptotic  limits, 
Equations  (A -2 -11}  and  (A-2-12),  hold  to  a  good  approximation  outside  of 
this  "boundary  layer."  This  is  indicated  implicitly  by  Figures  A-2-8  and 
A-2-9  for  the  continuum  growth  law.  Cases  5  and  6  represent  two  very 
dissimilar  initial  lognormal  distributions  studied  with  the  continuum 
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Figure  A-2-11.  Approach  to  asymptotic  limit  distribution  for  kinetic 
(free  molecule)  growth,  case  7,  Table  A-2-1. 


Linear  Kelvin  terms  for  kinetic  growth  (cases  7  and  8).  For  small  times, 


there  are  significant  differences  in  the  tails  of  the  distributions  for 


the  linear  and  nonlinear  cases. 


It  is  possible  to  form  estimates  for  the  characteristic  terms 


for  the  evolution  of  n(x,t)  by  Equations  ( A-2 -1  )-( A-2-5 ).  From  examina- 


ion  of  the  moments,  Figures  A-2 -4  through  A-2 -8,  it  is  clear  that  for 


the  initial  conditions  studied,  there  are  two  regimes  in  the  evolution  of 


n  '  x  ,  t  ■  . 


In  the  first,  the  number  concentration,  Mn(t),  remains  constant 


i  j  "  i  r.g  the  time  period  necessary  for  particles  to  evaporate  down  to  x* 


lost  from  the  distribution.  The  characteristic  time,  t .  (  for  thi: 


m  oar,  be  estimated  from  the  evaporation  rate,  .  ,  Equation  (A-2-2). 


r  continuum  growth: 


For  kinetic  growth: 


(m2(o)/m0(o)) 


lk2 


'vd(^,T/  m-,  ) 1  /2Dp(?  -■/.  T) 


where  M?(0)  is  the  initial  second  moment. 

In  any  event,  it  is  clear  that  increasing  the  aerosol  vapor 
pressure  over  that  used  in  our  simulations  will  lead  to  a  decrease  in  the 
time  necessary  to  approach  closely  the  asymptotic  regime.  Also,  a 
decrease  in  initial  particle  dispersion  decreases  the  time  necessary  to 
achieve  the  asymptotic  regime. 

From  our  numerical  studies  of  the  cond en sa ti on/ ev apo ra t ion 
problem  with  vapor  oonserva tion ,  we  conclude  the  following: 

1.  Starting  with  some  arbitrary  initial  distribution,  for  long  times,  an 
asymptotic  limit  distribution  is  approached  which  agrees,  for  the 
linearized  Kelvin  term,  with  the  analytical  similarity  solutions  for 
continuum  and  kinetic  (free  molecule)  growth  laws.  At  this  asymp¬ 
totic  limit,  details  of  the  initial  distributions  are  completely 
"  forgotten 

2.  The  time  evolution,  starting  with  some  arbitrary  initial  distribu¬ 
tion,  involves  two  epochs.  In  the  first,  the  total  number  concentra¬ 
tion  remains  sensibly  constant.  In  the  second,  particles  begin  to  be 
lost  from  the  distribution  and  the  rate  of  decrease  of  total  particle 
concentration  approaches  that  predicted  by  the  asymptotic  similarity 
theory . 

3.  The  time  necessary  to  achieve  the  asymptotic  similarity  solutions  for 
continuum  and  kinetic  (free  molecule)  growth  laws  increases  with 
increasing  dispersion  in  initial  particle  size  distribution. 

U .  For  long  times,  the  differences  in  the  distributions  resulting  from 
use  of  nonlinear  and  linear  Kelvin  terms  (Equations  (A-2-2)  and.  (A-2- 
3))  in  the  growth  law  are  confined  to  a  boundary  region  near  x* , 
outside  of  which  the  distributions  are  very  nearly  identical. 


5.  As  is  implicit  in  the  similarity  theory  development,  no  unique 
asymptotic  limit  distribution  occurs  in  the  transition  region  of 
Knudsen  numbers.  This  is  borne  out  by  our  numerical  simulations 
using  transition  region  growth  laws. 

6.  For  conditions  where  coagulation  may  be  neglected,  it  does  not  appear 
to  be  necessary  to  include  competitive  effects  in  describing  aerosol 
growth  by  condensa tion/evaporation. 

A -3 .  Aerosol  growth:  Ostwald  ripening  and  coalescence 

Numerous  studies33-^7  have  been  concerned  with  the  asymptotic 
stages  of  growth  by  a  homogeneous  nucleation  process.  Comparatively  less 
attention  has  been  given  the  associated  problem  of  growth  when  both 
condensation/evaporation  and  coalescence  occur.  For  aerosols  we  have 
presented  the  first  quantitative  study  of  the  Ostwald  ripening  problem 
with  a  realistic  coalescence  process. 

a.  Theory. 

We  have  investigated  the  isothermal  spatially  homogeneous 
growth  of  particles  suspended  in  supersaturated  monomer  vapora  nd  inert 
host  gas.  The  process  is  described  by  the  evolution  equation  for  the 
singlet  density  function  n(x,t): 

■n(  x,t)  '■  x/2 

-  +  —  [  Xx ,s)n(x  ,t )  ]  =  b( x-x' ,x' )n( x-x'  ,t) n( x'  ,t) dx' 

•t  X  x* 

-  n(x,t)  '  b(x’  ,x)n(  x’  ,t)dx’ 

x*  (A-3-1) 

where  n(x,t)dx  is  the  number  of  particles  having  masses  in  the  range  x,dx 
at  time,  t.  .(x,s)  is  the  growth  law  for  a  particle  by  the  condensa¬ 
tion/evaporation  process  and  b(x’,x)  is  the  coalescence  rate  coefficient 

hi 


for  two  particles  of  masses  x'  and  x.  x*  is  the  mass  of  the  smallest 


•SI 


^article  in  the  population  and  s  is  the  super  saturation  ratio  of  monomer 


vapor.  We  do  not  consider  here  competitive  growth  effects, ^7*^8  as  the 


theory  for  this  is  lacking  when  coalescence  is  significant. 


Equation  ( A  —3  — 1 )  is  coupled  to  the  mass  conservation  equation 


monomer  vapor: 


-v  —  =  -  /x  :+(x#3)n(  x,t)dx  -  x*  x*  ,s )  n(  x*  ,t) 
dt  x* 


( A-3-2 ) 


where  Cv  is  the  equilibrium  bulk  vapor  concentration . 


Equations  ( A  —3  —  1 )  and  (A-3-2)  obey  the  mass  conservation  law: 


—  (Cvs  +  M-j )  =  0 
dt 


wr.ere  M1  =  :  xn(x,t)dx  is  the  mass  concentration  of  particles.  An 

x* 


Important  additional  moment  of  n(x,t)  is  the  total  number  concentration, 


n(x,t)dx.  It  follows  from  Equation  (1)  that: 


—  =  .( x*  ,s  )  n(  x*  ,t )  -  /  /  b(  x’  ,x)n(x*  ,t)n(  x,t)dxdx' 

dt  x*  x* 


We  have  studied  the  evolution  of  n(x,t)  by  numerical  simulation 


using  techniques  described  in  detail  elsewhere.  These  simulations  use 


■  "a  continuum  diffusive  growth  laws: 


(x,s)  =  4"RDCv(s  -  exp(K/P)! 


where  x  =  (4r/3)R3d,  D  is  the  binary  diffusion  coefficient  of  monomer  in 
host  gas  and  K  is  the  Kelvin  coefficient  for  the  dependence  of  vapor 
pressure  on  droplet  curvature,  d  is  the  mass  density  of  a  particle  with 
radius  R.  The  rate  coefficient  for  coalescence  is  that  appropriate  to 
the  continuum  Brownian  process: 

1  1 

b(x',x)  =  (2-  T/3::)  ( x1  73  +  x'1^)  ( -  +  - ) 

x1/3  x.1/3 

where  >.T  is  the  thermal  energy  and  is  the  viscosity  coefficient  of  the 
host  gas. 

Equation  (A-3-1)  may  be  put  into  nondimensional  form  through 
the  substitutions:  p  =  xfx~0,  x  =  H  -  DCv(3/4  "dx^)1  /3t,  i  =  “xQn(x ,t) /N0- 
x0  is  the  initial  mean  particle  mass  and  N0  the  initial  number  concentra¬ 
tion  of  particles.  In  these  variables,  Equation  (a-3-1)  becomes: 


xfD/2  b^c-o’ 
0* 


L-]/3(s  -  ex p( K  '  / ,  1  x3  ) )  >■(  p,; )  ] 


(A -3-?  ) 


t  )  d. '  (p,t)  f  ,-)d. 


where 

=  kTN0x02/3/3-UDCv(3A^d)1/3 

The  coefficient  a  is  a  measure  of  the  relative  rates  of  coalescence  to 
condensa  tion/ evaporation.  In  Equation  (A-3-3),  fora  — >  ,  the  coales¬ 

cence  process  becomes  dominant  and  for  t — >  0,  only  the  conde-u  i- 
tion/ evaporation  process  remains  with  results  previous  described .35 , ^ 


A  question  addressed  by  others4 J»40  using  scaling  arguments  is 


which  mechanism  -  condense tion/evaporation  or  coalescence  -  will  be 
dominant  for  intermediate  values  of  The  conclusion  from  these 

analyses  is  that  asymptotically  coalescence  will  dominate  at  higher 
particle  mass  concentrations  -  that  is  for  larger  values  of  t .  We  show 
here  that  this  is  not  the  case  and  that  apparently  coalescence  can  only 
be  dominant  asymptotically  for  the  growth  laws  studied  for  vanishingly 
small  rates  of  the  condensation/ evaporation  process.  This  is  illustrated 
by  presenting  a  result  from  numerical  solution  of  Equations  (A-3-1)  and 
(A-3-2)  for  the  stated  growth  laws.  Parameters  used  yield:  x  =  0.31  ,  - 

5.3  x  1 0 3 1 ,  s(t=0)  =  10,  K/R0  =  1-0. 

Figure  A-3-1  shows  the  evolution  in  time  of  the  moments  N,  M1 

of  the  distribution  together  with  variation  with  time  of  the  interfacial 

diameter,  Dp«;  Dp*  =  2K/lnS.  After  an  initial  induction  per iod14^ ' 1,11  it 
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Figure  A-3-1  . 


can  be  noted  that  N ( t )  ~  t-^.  This  dependence  is  the  same  predicted  for 
asymptotic  continuum  diffusive  growth  by  the  condensation/evaporation 
process  according  to  the  Li  fshitz-Slyozov  (LS)  theory. 35  it  is  ai30  the 
same  dependence  found  both  numerically  and  by  similarity  theory^  for  the 
continuum  Brownian  coalescence  process.  The  increase  of  Dp*  with  time 
follows  the  law:  Dp*  ~  t^-3,  also  the  same  predicted  by  the  LS  theory 
for  only  the  condensation/ evapora tion  process.  Interestingly,  even  sub¬ 
stantial  coalescence  rates  do  not  alter  significantly  this  result  from 
condensa  tion/ evapora  tion  similarity  theory.  Since,  separately,  condensa¬ 
tion/evaporation  and  coalescence  both  give  the  result,  N(t)~t-1,  for  the 
continuum  diffusive  cases,  scaling  arguments  are  insufficient  to  estab¬ 
lish  which  of  the  two  mechanisms  is  dominant  asymptotically. 

Figure  A-3-2  shows  the  contribution  to  the  total  rate  of  change 
of  N,  dN/dt,  by  each  of  the  two  processes  -  condensation/evaporation  and 
coalescence.  The  continuous  curve  is  that  for  the  ratio  of  coalescence 
rate  of  change  of  N  to  the  total  rate  and  the  dashed  curve  that  for  the 
condensation/ evapora  tion  process  according  to  Equation  (A— 3—1 ).  Clearly, 
the  contribution  of  coalescence  to  the  change  of  N  decreases  with  time 
while  that  for  condensa  tion/  evaporation  increases.  Therefore  previous 
conjectures  on  the  dominance  of  either  of  the  two  processes  appear  to  be 
incorrect  for  these  continuum  diffusive  growth  processes. 


these  binary  aerosols  is  considered  briefly  in  terms  of  suitable 
mations  of  the  growth  process. 


a .  Theory  . 


There  have  been  a  number  of  discussions  of  the  growth  of  r  * 
aerosols  (e.g.  Fr  iedlander^ t  Gelbard  and  Seinfeld^).  jt  is  suf :  .  .. 
to  note  here  that  the  exact  description  of  condensation  growth  of  t 
aerosols  is  giver,  by  the  differential-difference  equation: 


<  1  ) 


dni,j/dt  =  -(  ‘i,j  + 


(2) 


U) 


j)  ni » j  +  1  i— 1 »jni— 1 »j  +  ai, j-1ni,j-1 


. 

{  i,  J 


i » j }  ni,j 


(1) 


'  i  +  1  » J  i  +  1  *  j 


a(2)  n. 

►i,j+l  nifj+1» 


(  A- 


where  n^j  is  the  concentration  of  particles  having  i  molecule 

species  1  and  j  molecules  of  species  2.  Here  'if!  and  ft  are  the  : 

1  t  J  1  f  J 

of  addition  of  a  molecule  of  species  1  and  2,  respectively  to  a  .  1  j 
containing  i,j  molecules.  off  and  rff  are  the  rates  of  los 

>  J  1  9  J 

evaporation  of  a  molecule  of  species  1  and  2,  respectively,  by  a  cli 
containing  i,j  molecules. 

If  it  is  supposed  that  species  2  is  involatile,  then 

(2)  _ 

.  .-0  by  definition.  Therefore  equation  ( A  —  4  — 1 )  becomes  on  dm- 

i  >  i 

the  remaining  supe ’script: 


d  n  ,  .  /  d  t  -  -  ( 


i  .  J 


1- 


j  n  i  —  1  ,  j  + 


;  V.  , 


In  the  nontext  of  thin  note  equation  ; 


ultrafine  aerosol  phase  tne  coefficients 


_  ■ ,  .  .  •  are  evaluated  for  the 

free  molecule  regime  of  Knud sen  number  assuming  the  size  of  a  particle  is 
large  compared  to  that  of  a  condensing  molecule: 

iifj  =  (8~-  T/m1)1/2(3v1/4-:)2/3(i+jv?/v1)-/3n1,  '  A-s-3  ) 

where  kT  is  the  thermal  energy,  and  m-,  the  molecular  mass  of  species  i; 
v.  and  v,  are,  respectively,  the  volume  per  molecule  of  species  1  an  i  2, 
and  for  an  ideal  mixture  would  be  independent  of  i  and  j,  the  numbers  of 
molecules  of  species  1  and  2  making  up  a  particle,  n,  is  the  number 
concentration  of  species  1  molecules  in  the  vapor  at  large  distances  from 
a  cluster. 

Similarly,  for  the  evaporation  coefficient: 

I(1  =  (8  kT/m,  )1/2(3v1/4~)2/3  (  i+  jv2/Vl  )2/3n* ,  (A-4-4) 

where  r.^  is  the  concentration  of  species  1  molecules  leaving  the  particle 
surface.  From  Defay  et  al.^1: 

L n '  r. ,  / r. i  )  =  2  v1  /kT(  3 v  1  /4 -) 1  /3  (  i+  jv2/v1  ) 1  /3  +  In-., 

-  ln(1  +  j/i)# 

where  is  the  surface  tension  of  a  particle  of  the  indicated  composition 
and  ^  is  the  activity  coefficient  of  species  1  in  the  particle  phase; 
r.*"1  is  the  surface  concentration  for  a  plane  surface  given  in  terms  of 
the  vapor  pressure  of  pure  species  1.  On  the  assumption  of  ideality  and 
large  particles,  it  is  easy  to  show  that  (5)  reduces  to  the  familiar  form 
used  in  cloud  physics  to  derive  the  so-called  Kohler  curves  (Rogers^): 

=  (3  .  T/m1 J1 /-( 3v1  /4  )2/3n+Ot 


<  A-4-5) 


=  (  i-t-A  )2/3exp(K/(  i+A  )1 /3+in  ■ , ) /( 1  +  j/ i) 

K  =  2  ’v  1  /.  t>-3v1  /4  Ti)1  73 
A  =  jv2/vr 

Owing  to  the  lack  of  a  suitable  theory  of  such  differential- 
difference  equations  and  in  the  interest  of  computational  simplicity,  it 
has  for  some  time  ( Abraham3® )  been  customary  to  study  continuous  approxi¬ 
mations  of  equation  (A -4 -2).  Introducing  the  displacement  operator 

exp(  +  a  --)  n(i)  =  n(iaa) 


and  regarding  i  as  a  continuous  variable,  (3)  becomes: 


n(  i  ,  j ) 


1  2 

—  ('—•■)  n(  i  ,  j)  +  -  — —  )n(  i  ,j) 

■i  2  i2 


1  3 

- -  )n(i,j) 

6  i3 


24  1 


(  ,  +■  )  n(  i ,  j 


(  A  -4  -6  ) 


Tne  convergence  of  equation  (A-4-5)  is  usually  assumed  and  it  has  been 

customary  (Abraham3®)  to  truncate  the  series  at  the  second-order  term  as 

a  suitable  approximation  to  such  differential-difference  equations; 

n(i,j)  1  2 

-  =  -  —  (  )n( i, j)  +  -  — -  !  +  ) n( i ,  -  .  A-<-" 

i  2  12 

It  will  be  shown  below  that  this  truncation  is  no  t  g«n  or  a '  '.7  •>  •  . 


7h-  presence  of  nn  invol  utile  component  makes  possible  the 
exist-:.  '*  •  ••  :  1  itr  ium  1 :  st ribu tions  in  certain  regions  of  vapor  con¬ 

cent  rat  :  or..  .  r  » *  :  :  t :  o  n  to  their  fundamental  interest  these  distribu¬ 
tions  l  .  *  e«:  if  the  validity  of  the  expansion,  equation  ( A  — —6 ) . 

in-  ■  ; :  :  r  nun  distributions  discussed  here  are  found  from  the 
solution  '  '-qiation  'A  -d-?:  for  dn{  ,/dt  =  0.  The  first-order  approxi- 
m  at  1  on  if  equ  i  *  i  or.  .A—  J '  is: 

hi  i,  ’’/  =  )n(  i ,  j '  /  -i. 

The  equ  id  ifcriu~  sol  ution  of  this  equation  for  a  single  value  of  j  is  the 
delta  fir.  oti^n: 

r.;  i ,  j  '  :  N  ■  ;  i-  i#  )  ,  (  A-d  -8  ) 

where  i*  is  the  real  root  of  •-  =0  and  N  is  the  total  number  concentra¬ 
tion.  For  variable  j  and  A  one  has  a  sum  of  such  delta  functions. 

Equation  A  —  -'  of  course  does  not  agree  with  results  from  equation  (A- 

4 -J  .  .  ~h-  r  ;  in  reality  dispersion  in  the  distribution  owing  to  the 

d  i  r  .  r  •  •  t e  n  i*.  in*  f  growth  process. 

in  order  to  include  these  dispersive  effects,  one  uses  equation 
A - ’  ,  w n i  •' r  it  -'qcilibrium  is  a  steady- state  Fokker-Planek  equation: 

d  *  :n'1i-  -  ’1  —  'n/di  r  Q.  A-d  -d  1 

,  r  :v  solution  of  inis  equation  for  :  h---  binary 


where  is  me  normali  zation  factor: 


+  exp 


d  i '  )  c  i 


In  the  iisc  where  the  condensing  vapor  is  conserved,  the  equation  :  . 
saturation  ratio  s  and  equation  (  A  -4  —9  )  are  coupled.  In  this  :  n  v. 
equation  is  only  a  formal  solution  as  it  is  now  a  n -..ij; 

integral  equation  for  n'i,j).  There  are,  however,  no  new  features  of 
equilibrium  distributions  introduced  by  considering  vapor  conserve 
and  the  case  of  fixed  vapor  concentration  is  usually  the  most  impor 


Figure  A-'t-l  3hows  a  typical  example  of  the  reasonable  ag 


ment  possible  between  equation  (A -4 -10)  and  the  exact  equilibrium  s 
tion  of  equation  '  A— *4  -1  i  However,  this  agreement  is  not  universal. 


s  1  or  K 


►.•qua  tion  A-U-10)  no  longer  provides  a  use:  ;1 


proxni  a  1 1  o  n  to  the  exam  result  from  equation  (  A  -U  -?  ) .  As  sc 


Figures  A-s-s  and  A— *  -  :  where  K  -  C.  1,  the  distribution  giver;  by  - 
;A-h-9)  is  o'  1  ••v-r  abl  e.  Fven  the  next  level  of  approx :  •:  .* 


G 


Figure  A-4-2.  Comparison  of  equilibrium  distributions  calculated  from 
exact  difference  equation  (equation  (A-4-2),  squares) 
Fokker-Planck  approximation  (equation  (A-4-9),  circles) 
fourth-order  approximation  (equation  (A-4-12),  tri¬ 
angles).  Parameters:  s=0.99;  K=0.1;  A  =  1 0;  j  =  5. 
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Figure  A -4 -4. 


Comparison  of  equilibrium  distributions  calculated  from: 
exact  difference  equation  (equation  (A-4-2),  squares); 
Fokker-Planck  approximation  (equation  (A-4-9),  circles); 
fourth-order  approximation  (equation  (A-4-12),  tri¬ 
angles).  Parameters:  s=0.9999;  K  =  1.0;  A  =  1 00;  j  =  50. 


For  particular  values  of  the  parameters  in  a  and  3,  as  s  in¬ 
creases  above  1,  equilibrium  distributions  are  not  possible  and  un¬ 
rest  rio  te  i  growth  of  an  initial  distribution  occurs.  It  is  conventional 


to  lesoript' 


;rowtn  of  an  aerosol  by  condensation  with  the 


••q  i  •  t  ion: 


n/  t  -  -(.-.)  n/  i 


( A -4-1  3) 


■,  ere  first  term  in  the  expansion  of  equation  (A-4-2) 
■a,-..--  investigators  have  neglected  the  second-order 


j n  for  this  neglect  has  been  discussed  by 


Here  the  question  of  neglect  of  the  diffusive  term  of  equatio' 

C  A  — —7 )  for  modelling  unrestricted  condensation  growth  is  examined  using 
a  lognormal  model  for  the  distribution:  n(i)  =  (N /,  Jmi  1  n-; )  exp  (-  1  n  i  i>  ' 
where  i  and  -  are,  respectively,  the  geometric  mean  number  and  stand  •  ■-  1 
deviation.  A  convenient  criterion  for  this  neglect  based  on  the  rale., 
of  the  derivatives  d- (  ■+  )n/di^,  d(.x-3  )n/di,  evaluated  at  i  would  to  : 
1/ i(  lnt  )2<<  i ,  which  will  De  satisfied  for  i  and  sufficiently  large.  A 
more  detailed  study  would  need  to  compare  terms  in  equation  (A-^-7'  for-  a 
particular  growth  law. 

In  simulations  of  particle  growth  for  other  aerosol  growth 
models  (to  be  presented  elsehwere)  it  is  found  that  equation  (a-4-7; 
gives  an  excellent  representation  of  the  exact  equation  (a-4-1)  from  the 
nucleation  step  to  growth  to  large  particles.  The  necessity  for  discrete- 
continuous  simulations  is  not  therefore  always  clear.  One  can  argue  :r 
favour  of  including  the  diffusive  term  in  simulations  of  aerosol  growth 
under  all  conditions  in  which  growth  may  occur  at  sharp  fronts  or  boun¬ 
daries  of  distributions.  This  would  be  in  the  spirit  of  inclusion  of 
viscosity  terms  in  the  hyperbolic  conservation  equati'...s  for  numerical 
treatment  of  shock  propagation  (von  Neumann  and  Richtmeyer^). 

An  interesting,  if  intuitively  obvious  aspect  of  unrestricted 
growth  of  binary  aerosols  is  the  appearance  of  multimodal  distributions 
from  an  initially  urii.T''dal  distribution  when  the  aerosol  is  art  ex  t  -  -  r: 
mixture  of  a  single  volatile  compojnent  and  two  or  more  invo, 
sp-=‘c  ies.  An  external  mixtures  means  that  the  particles  k  1  v  •• 
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Figure  A-4-5.  Development  of  bimodal  distribution  during  condensation 
growth  from  initial  unimodal  distribution  of  mixture  of 
two  aerosols,  I  and  II.  Parameters:  Aerosol  I:  s  =  2; 
R  =0 . 1 ;  A=25;j=20.  Aerosol  II:  s=2;  R=0.4;  A=25;  j=20. 
In  both  cases  y,  taken  from  data  for  Hg-Bi  mixtures  at 
320  °C . 

aerosol  contain  different  masses  of  the  involatile  species.  Figure  A-4-5 
shows  the  development  of  a  bimodal  distribution  for  an  initial  aerosol 
consisting  of  two  different  involatile  species,  although  all  the 
parameters  in  -t  and  p  are  the  same  except  for  K:  A  =  25,  j  =  20  and  Y1 
is  taken  from  experimental  data  for  Hg-3i  mixtures  at  320°C.  In  the 
simulation,  K  =0.1  and  0.4.  In  general  for  differences  in  A,  K  and  >i 
even  for  the  same  values  of  j  multimodal  functions  can  arise  for  long 
periods  of  unrestricted  growth. 

From  the  results  of  this  study  it  is  concluded  that  for  binary 
ultrafine  particles  with  one  involatile  component: 

(1)  Equilibrium  distributions  are  approximated  satisfactorily 
by  the  steady-state  Fokker-Planck  equation,  equation  t'A-!l- 
9),  for  sufficiently  large  K  and  s  sufficiently  less  than 
1.  In  practice,  K  >  1  and  s  <  0.4. 


(2)  As  s  approaches  1  and  K  approaches  0  the  series  expansion 
of  the  difference  equation  does  not  provide  a  useful 
approximation  to  the  equilibrium  distributions. 

(3)  Neglect  of  the  second-order  diffusive  term  in  describing 
unrestricted  growth  of  binary  ultrafine  particles  is 
justified  only  for  sufficiently  large  clusters  with 
relatively  broad  distributions  in  the  absence  of  sharp 
fronts  or  boundaries  in  the  growth  process. 

(4)  Under  restricted  condensation  growth,  initial  differences 
in  properties  of  ultrafine  nuclei  may  lead  for  some  time 
to  the  appearance  of  multimodal  distributions  from  initial 
unimodal  distributions. 

A -5 .  Growth  of  ferromagnetic  particles  in  gas  phase 

The  production  of  finely  divided  particles  by  evaporation  into 
a  background  gas  dates  back  to  the  early  1  930's  when  Pfund^  made  the 
famous  and  well=known  metallic  blacks  of  a  series  of  metals  by 
evaporating  into  air  at  a  pressure  of  1  Torr.  Since  that  time  many  studios 
have  been  carried  out  in  this  process. -^“^5  jn  Spite  of  the  large  amount 
of  previous  work  and  the  frequent  use  of  the  evaporation  method  in  ex¬ 
perimental  procedures,  the  detailed  processes  leading  to  the  formation  of 
small  primary  particles  in  the  gas  surrounding  a  heated  sample  have  not 
been  well  understood.  Primary  particles  are  defined  as  the  spherical 


particles  observed,  either  as  single  isolated  particles  or  as  their  chain 


Tr.  work  reporte 1  here,  experimental  studies  were  undertaken 
to  confirm  observations  of  previous  investigators  on  the  formation  of 
primary  particles  and  their  chain  agglomerates  when  metal  samples  are 
neated  in  ir.-'rt  gas  at  reduced  pressures.  These  experimental  studies  are 
reported  first,  together  with  a  summary  of  the  observations  and  discus¬ 
sion  of  the  results.  Next,  previous  theoretical  analyses  on  the  forma¬ 
tion  of  primary  particles  are  summarized  followed  by  the  development  of  a 
new  theoretical  analysis  which  explains  the  present  observations  as  well 
as  those  of  earlier  investigators  on  the  formation  of  primary  particles 
in  evaporation  chambers. 

a.  Experiment. 

The  formation  of  primary  particles  and  their  chain  aggregates 
was  studied  experimentally  using  an  evaporation  chamber  for  particle 
generation  ana  electron  microscopy  for  determination  of  particle  mor¬ 
phology  and  size  distribution.  The  experimental  systems  and  procedures 
used  in  the  study  are  outlined  here. 

The  evaporation  chamber  is  shown  schematical lyin  Figure  A  —5—1 . 
The  chamber  was  similar  to  that  used  in  conventional  vacuum  evaporation. 
Here,  the  chamber  was  a  Pyrex  "T"  with  a  diameter  of  75  mm  and  a  length 
of  755  mm  for  the  horizontal  portion.  Feedthroughs  were  provided  in  the 
bas-*  pi  ate  for  host  gas  and  electrical  power  and  a  vacuum  line  was 
attached  to  one  of  the  side  plates.  The  heater  for  the  evaporation 
source  was  tungsten  w :  -  e  25  mm  in  length  and  0.5  mm  in  diameter.  The 


detail 

Figure  A-5-1  .  Schematic  diagram  of  experimental  metal  evaporation 
chamber.  Points  1,  2,  and  3  are  particle  sampling  posi¬ 
tions  . 

wire  was  bent  to  form  a  "V"  in  its  center  where  the  metal  sample  was 
placed.  This  sample  typically  had  a  mass  of  2-3  mg.  Usually,  the 
chamber  was  evacuated  down  to  IE  -  5  Torr  initially,  then  filled  with 
host  gas  to  300  Torr,  reevacuated,  and  then  refilled  with  host  gas  to  the 
desired  pressure.  The  heating  current  was  applied  for  about  ten  seconds 
to  evaporate  the  metal  sample.  The  pressure  inside  the  chamber  was 
measured  with  ionization  and  thermocouple  gauges.  The  temperature  of  the 
evaporation  source  was  determined  with  an  optical  pyrometer.  Particles 
resulting  from  evaporation  of  the  metal  sample  were  collected  on  carbon- 
coated  transmission  electron  microscope  grids  at  various  positions  in  the 
chamber  relative  to  the  source.  A  shutter  was  used  to  control  the  am  met 
of  deposit  on  a  grid.  The  entire  region  could  be  subjected  to  magnetic 
fields  up  to  1200  G.  Grids  with  deposited  particles  were  stored  in 
nonoxidizing  conditions  and  observed  using  a  transmission  electron 
microscope  (JEDL  JEM -200  CX)  as  well  as  a  scanning  electron  m  icron  ; 
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(JEOL  JSM-35C).  For  determination  of  primary  particle  size 
distributions,  a  sem  iautomated  procedure  developed  by  Hazlett  et  al.^ 
was  used.  In  this  de term  ination,  images  on  negatives  From  electron 
micrographs  of  particle  samples  were  digitized  with  a  video  camera. 
Analysis  proceeded  with  edge  enhancement  and  edge  detection  algorithms 
which  ask  For  visual  user  conFirmation  oF  particle  contour  and  acceptance 
oF  particle  size. 

b.  Results  . 

The  Following  qualitative  observations  correspond  to  those  oF 
previous  workers.  Upon  heating  oF  the  source,  a  well-deFined  spherical 
smoke  shell  was  observed  to  Form  concentrically  with  the  source. 
Particulate  samples  were  readily  obtained  when  the  specimen  grids  were 
placed  on  or  near  this  smoke  shell.  When  a  grid  was  placed  inside  the 
smoke  shell,  the  deposits  resembled  those  seen  in  thin- film  micrographs. 
This  strongly  implies  that  particles  are  not  Formed  in  the  gas  phase 
inside  the  smoke  shell  and  that  nucleation  and  particle  growth  occur 
directly  on  the  grids  in  this  region.  As  the  host-gas  pressure  was 
lowered  below  1  Torr,  the  visual  smoke  shell  disappeared.  At  host-gas 
pressure  below  0.01  Torr,  it  was  inferred  that  no  particles  were  Formed 
in  the  gas  phase  as  the  deposits  on  the  grids  had  the  characteristic 
appearance  of  thin-filra  deposits.  At  higher  pressures  in  the  vicinity  of 
the  smoke  shell,  the  particles  formed  were  more  or  less  linear  chains 
even  in  the  absence  of  an  external  magnetic  field;  inside  the  smoke 


shell,  particle  chains  were  rarely  observe!.  T  r;  an!  ►».<  t  *;r  ior  to  the 
smoke  shell,  application  of  an  external  magnetic  fie  1 1  produced  linear 
chains.  However,  it  was  found  that  small  scale  linearity  of  the  chains 
was  not  greatly  improved  even  at  1200  G. 

The  size  of  the  primary  particles  making  up  the  linear  particle 
chains  depended  upon  the  host-gas  pressure,  the  nature  of  the  host  gas, 
and  the  source  temperature.  The  mean  size  of  the  primary  particles 
increased  with  host-gas  molecular  weight.  Mean  primary  particle  size 
increased  as  the  host  gas  was  changed  from  helium  to  argon  to  xenon. 
Primary  particle  size  also  increased  as  the  temperature  was  increased 
above  1800°C.  Below  approximately  1800°C  very  few  particles  could  De 
observed  on  the  grids.  The  primary  particles  have  quite  narrow  size 
distributions.  Figure  A-5-2  shows  primary  particle  size  distributions 
for  some  representative  chains.  The  three  distributions  represent 
particles  produced  in  the  chamber  at  pressures  of  2,  8,  and  20  Torr  and  a 
temperature  of  1 950°C  (in  order  of  increasing  mean  particle  diameter). 

c .  Theory . 

There  have  been  a  number  of  previous  studies  or,  the  formation 
of  particles  in  evaporation  chambers.  The  previous  theoretical  analyses 
have  not  been  successful  in  explaining  all  of  the  qualitative  observa¬ 
tions  cited  above.  For  example,  Kawamura^?  posed  a  simple  model  in  which 
the  vapor  super  saturation  in  the  evaporation  chamber  determined  the  size 
of  the  metallic  fine  particles.  This  model  explained  the  dependence  of 
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DIMENSIONLESS  DIAMETER, D/DAVG 

A-5-2.  Experimental  size  distributions  of  primary  particles 
formed  in  evaporation  chamber.  D  is  particle  diameter, 
DAVG  is  the  number  mean  diameter.  Curve  A  (line): 
DAVG  =  12i4  A,  standard  dev iation=0.22,  host-gas  pressure=20 
Torr;  Curve  B  (dash);  DAVG=99  8,  standard  dev ia tion  =  0.25, 
host-gas  pressure=6  Torr  Curve  C  (dot):  DAVG=6  %, 

standard  dev  ia  tion=0.265,  host-gas  pressure=2  Torr. 


ary  particle  size  on  the  chamber  pressure,  but  failed  to  explain  its 
science  on  source  temperature.  The  model  also  did  not  contain  a 
ription  of  particle  formation  and  growth.  Ozawa^®  obtained  particle 
i  J i str  ibutions  by  solving  the  growth  equations  for  the  primary 
ides  using  steady-state  solutions  for  the  spatial  distributions  of 
•■.nature  and  vapor  concentration.  This  model  predicted  that  primary 
ide  formation  occurred  in  a  narrow  region  close  tothe  cooled  wall  of 
■camber,  a  finding  contradicted  by  experimental  observations.  In 
dor.,  the  particle  growth  equations  were  not  coupled  to  the  energy 
v »por  conservation  equations  and  did  not  include  any  particle 
■r  art  mechanisms.  Yatsuya  et  alA-  observed  a  convective  flow  around 


■  j  jrce. 


a  heate-i  tungsten  boat,  and  concluded  that  the  smoke  shell 


observed  occurred  m  a  reg icr.  wr.ere  tru 


lively  coci  host  gas  net  -  r  a  :  •.  I  flow  ■:  m  •.  t  : .  .  -  :•  v 

out  from  the  source.  Howe.s-r,  when  the  turigsten  tout  was  re,  - 
tungsten  wire  as  heater,  n  convective  f  '  w  was  observed  •«  1  ♦ 


smoke  she 2 1  was 


that  primary  particles  grow  mostly  by 


in  the  sir  rt 


conclusion  based  not  on  simulation  of  particle  dynamics,  but  p.n  . 

son  of  residence  time  in  the  smoke  shell  with  characteristic  con 
times.  Similar  conclusions  were  reached  by  Granqvist  et  al.,’Jh  ag  a:  •.  t 
the  absence  of  simulation  of  particle  dynamics. 

In  this  section,  a  model  is  developed  for  the  formation  o 
primary  particles  in  evaporation  chambers  which  contains  all  f  r  t 
relevant  dynamical  processes  of  particle  formation,  growth,  as 
transport.  The  physical  system  to  be  modeled  consists  of  a  spr.e":  a 
chamber  of  radius  r0  filled  with  Inert  host  gas  at  some  pressure  P. 
spherical  source  is  located  at  the  center  of  the  chamber  with  radius  r 
and  is  maintained  at  temperature  Tg  while  the  chamber  walls  are  kept  . 
temperature  Tg.  At  the  surface  of  the  spherical  source,  the  partia 
pressure  of  the  metal  vapor  is  assumed  to  be  the  equilibrium  v  ..  • 
pressure  of  the  metal  at  temperature  T  The  contribution  of  the  m  i  : 

vapor  to  the  chamber  pressure  is  assumed  to  be  negligible.  Fro:-;  ‘  h 
source,  both  he  a  t  a  n  d  m  e  ‘  i  1  vapor  d  i  f  f 1 :  s  <?  radially  o  u  t  wa  r  1 . 
gradients  of  temperature  and  vapor  result  in  spatial  gradients 
saturation  which  give  ris--  t  homogeneous  nuc  1  eat  Lon  and  r  • 
primary  particles  whose  r  ite;-  vary  with  5 ;  d  ance  t  rom  the  sc  ,>-  •>  . 


•P 


gradients  of  particles  and  temperature  give  rise  to  radial  particle- 
fluxes  due  to  Brownian  diffusion  and  thermophoresis.  Based  nr. 
observations  and  simulation  results  to  be  presented,  th-5  formation  of 
particle  chains  only  occurs  at  larger  radial  distances  after  tn-  jr^wti. 
of  the  primary  particles  is  completed. 

T.oe  particle  formation,  growth,  and  transport  processes 
outlined  are  described  by  the  following  evolution  equation  for  the 
primary  particle  size  distribution  function  n(x,r,t),  where  n(x,r,t)dx  is 
the  number  of  particles  having  masses  in  the  range  x,dx  at  radial 
position  r  at  time  t : 


r. x  ,  r  ,  t } 


t 


1  -(r2N) 

(__)  .  +  o.5  b(  x  ,x' )  n(  x-x' ,r  ,t )  n(  x' ,r  ,t )  dx  ' 

r2  r  x* 


-  n  (  x  ,  r  ,  t )  b(  x’  ,x)n(x'  ,r  ,t)dx' 

x* 


[r(x,cv,t)n(x,r,t)] 
•  x 


( A-5-1  ) 
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Here  N  is  the  total  diffusional  flux  of  particles.  The  second  and  third 
terms  on  the  right-hand  side  account  for  particle  coalescence  and  the 
fourth  term  for  growth  by  condensation  of  vapor.  R  is  the  rate  of 
production  of  particles  by  homogeneous  nucleation.  For  all  experimental 
runs,  the  particle  Knudsen  number  Kn  is  such  that  Kn>>10.  Consequently, 
ail  particle  growth  and  transport  processes  are  described  by  free 
molecular  theory.22  The  diffusional  flux  of  particles  is  given  by  the 


expression 


wher--  is  the  free  molecular  Brownian  •.  •■••••  .and  % 

f 

the  tree  molecular  thermophoreetic  velocity.’’  '  im  i’.  arly  ir.  ;  •.  *  •  r. 
5-1),  the  coalescence  coeffi"  ient  blx,x'  a  r.  j  *  h--  cop  Jens.-,  t :  r.n  '  ~v  apo 
t  ion  coefficient  :  are  given  by  the  fret  ",  ;  <•  ■ul  ar  ex  pres.-.; '  r.s. 
homogeneous  nucleation  rate  P  is  assumed  ’  '  r.  y  *  .he  ->xpr>*ss  lor 

P  =  P*  [x  -  x*'  r,‘  )  ,  A-c 

where  is  the  delta  function,  x*  is  the  critical  nucleus  si z- ,  and  P1 
the  ra*-1  if  appearance  of  critical  nuclei  ty  to-  ■'lass:  r»’.  t. r*- 

of  y  - 1  m  er-?ec  ker-for  i  r.g-Zeldov  icn. '  ^  Equat :  -  A-c-i'  ir  v>up  l^c  to 

conservation  equation  for  metallic  vapor,  wh .  n-.n  tr.e  f-rr 

■  cv  )  n  ,  t )  1  (r^Ny) 

-  1  =  —  ( —  )  ■  —  {  x ,  c  y  t  r  N  r.  x  ,  r  ,  *  d  x  -  x  ■*  t  x  .  A  — 

t  r'  r  x*  x* 

The  vapor  diffusion  flux  is 


Ny  = 


where  £  is  the  diffusion  coefficient  for  vapor,  the  molar  Conner, 
t  i  or.  of  the  host  gas,  p_  and  M,  the  density  »r.  i  m  ol  ecu!  ar  w-  .sot  -  f 
host  gas,  Mv  the  molecular  weight  of  vapor,  an:  k-r  th-=*  thermal  d  1  ffu 
ratio  of  vapor  in  the  host  gas.  In  all  calculations  t-mp-ra* 

dependence  of  the  transport  coefficients  is  accoanted  for. 

Finally,  the  energy  equation  is  coupled  to  Equations  ’  A-‘ 


and  ■  A  -  5  — h  '■ 


t  has  the  form 


-  i  . 


M,  r,a)  -  xr'o  x,r,t)dx 

x* 

wr.-re  *  :  -=  .tr:.-  itrt"  y  >v.>  1  :;um  t  er ,  M  ;_•?  one  t:< 

of  particles,  .*•*.  <  tr.*.’  :r.  <.s^  concertina  tic::  of  parti  cl  *‘3,  and  V 
m  oz  or.  t  related  to  the  variance.  It  ;s  clear  from  Equ  lticr. 
straightforward  application  of  the  moment  net  hoi  to  Equation 
not  result  in  a  closed  set  of  equations  for  tnese  tnre< 
fractional  moments  arise  from  the  various  particle  lynarr.  ioal  > 
which  involve  fractional  powers  of  the  particle  mass.  Mor 
solutions  of  this  problem  are  possible;  however,  here  it  is  su 
employ  a  very  approximate  theory  since  only  the  qualitative 
the  present  and  previous  experimental  results  are  sought.  T. 
the  simulations  presented  here  are  for  times  less  than  the  o 
ms  so  that  particles  io  not  grow  larger  than  approximately  2-n 
In  this  case,  the  particle  dynamical  coefficients  are  taker, 
stants  and  the  average  particle  diameter  can  be  expressed  in  t; 
three  moments  with  the  assumption  of  a  lognormal  particle  si  z 
t  i  o  n  for  n '  x  ,  t : : 

-  1  >2 

nlx.tl  =  M  -  __  (1/xln  )  exp  [ -  ( lnx  -  Inxav  )^/.?  In'  1, 

where  xav  is  t.ne  geometric  mean  particle  mass  an  i  in.  a..- 
standard  :e / i  : f  i  •■r. .  From  Equations  and  'A  ;  t  • 


(A-5-10  , 


M  ,  ;  vj^xav  exp  O  In  -./?), 

M-,  -  Mg  xav  -  exp  (18  lr.2..).  f  A-5-1  "! 

xi.  ana  ear.  be  eliminated  from  F.quations  (A-5-10)  and  ( A  —  5  —  1  1  1  to 
obtain  the  geometric  mean  particle  diameter  d  • 


Jn  =  (6M 


-A  /2 
0  " 


A-8- 1 


Particle  diameters  are  calculated  only  in  the  region  of  the  smoke  shell 
since  the  numerical  error  in  the  calculation  will  be  small  there.  At 
other  positions  the  particle  concentrations  are  found  to  be  very  low  and 
the  numerical  error  due  to  the  very  stiff  nature  of  this  problem  is 
large  . 


Discussion  . 


The  moment  equations  for  Mg,  m^,  and  M2,  the  mass  conservation 
equation,  and  the  coupling  equation  between  local  mean  particle  tempera¬ 
ture  and  host -gas  temperature  are  solved  numerically.  Spatial  distribu¬ 
tions  of  vapor  concentration,  temperature,  supersaturation,  and  par-  i  -1 
moments  are  reported  as  a  function  of  the  logarithm  of  the  dimer.?:  : 
distance  from  the  source,  In  ( r/rs)/ln(  rg/rg)» 

As  shown  in  Figure  A -5-3,  the  distribution  'f  • 

has  a  maximum  between  the  source  and  chamber  wall’.  :  * 
the  super sa  burst  ion  value  at  the  maximum  r  :n  * 
temperature,  tut  is  relatively  inruns:  -  iv- 
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Figure  A-5-3*  Distribution  of  supersaturation  ratio  in  model  evaporation 
chamber  with  time  and  host-gas  pressure  as  parameter. 
Ts=2000  K,  nitrogen  host  gas.  Curve  A:  P=500  dyn/cm2, 
t=0.05  ms;  Curve  B:  P=10000  dyn/cm2,  t=0.05  ms;  Curve  C: 
P  =500  dyn/ cm2,  t=0.1  ms;  Curve  D;  P=10000  dyn/cm2,  t  =  0.1 
ms. 


properties.  It  is  in  this  super saturation  maximum  that  particle 
formation  and  growth  is  dominant.  As  a  consequence  of  this  maximum 
point,  the  particle  concentration  is  vanishingly  small  everywhere  but  in 
a  narrow  spatial  region,  termed  above  the  so-called  "smoke  shell." 
Figure  A-5-4  shows  that  the  vapor  concentration  curves  exhibit  a  cross¬ 
over  as  the  host-gas  pressure  is  changed.  As  host-gas  pressure  is  in¬ 
creased,  with  reference  to  Figure  A-5-3  the  supersaturation  maximum  moves 


toward  the  source  and,  as  shown  in  Figure  A-5-5,  the  smoke  shell  becomes 
thinner.  The  other  particle  moments  display  similar  trends  to  that 
indicated  for  the  particle  number  concentration  in  Figure  A-5-5. 

The  super sa turation  ratio  and  particle  moments  are  very 
sensitive  to  the  source  temperature,  as  would  be  expected  from  the 
exponential  dependence  of  metal  vapor  pressure  on  temperature  and  the 
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Figure  A-5-^.  Distribution  of  supersa tura tion  metal  vapor  concentration 
in  model  evaporation  chamber.  Crossover  effect  with 
variation  of  host  gas  pressure  is  shown.  T3=2000  K, 
nitrogen  host  gas.  Curve  A:  P=500  dyn/cm2,  t=0.05  ms; 
Curve  B:  P=10000  dyn/cm2,  t  =  0.05  ms;  Curve  C:  P=500 
dyn/cm2,  t=0.1  ms;  Curve  D:  P  =  10000  dyn/cm2,  t=0.1  ms. 


LN  (REDUCED  DISTANCE  FROM  SOURCE) 

Figure  A-5-5.  Distribution  of  total  number  concentration  of  primary 
particles  in  model  evaporation  chamber  with  pressure  and 
time  as  parameters.  Narrow  concentration  distributions 
l  correspond  to  experimentally  observed  "smoke  shell"  sur¬ 

rounding  the  vapor  source.  Tg=2000  K,  nitrogen  host  gas. 
Curve  A:  P=500  dyn/cm2,  t=0.05  ms;  Curve  B:  P=10000 
dyn/cm2,  t  =  0.05  ms;  Curve  C:  P=500  dyn/cm2,  t  =  0 . 1  ms; 
Curve  D:  P  =  10000  dyn/cm2,  t=0.1  ms. 
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large  variation  of  particle  nucleation  rate  with  supersaturation  ratio. 


At  lower  source  temperatures  supersaturation  ratios  are  not  large  enough 


to  produce  high  concentrations  of  particles  and  a  smoke  shell  about  the 


source  will  not  be  observed.  However,  as  shown  by  Figure  A-5-6,  if 


sufficiently  high  particle  concentrations  exist  so  that  the  smoke  shell 


could  be  observed  at  all,  the  width  of  the  smoke  shell  will  not  vary  with 


source  temperature;  this  phenomenon  has  been  observed  experimentally  by 


Yatsuya  et  al.,  but  not  explained  before  now. 


Mean  particle  diameters  have  been  calculated  from  the  particle 


moments  only  in  the  region  of  high  particle  concentration,  the  smoke 


shell.  The  calculations  show  that  increases  of  host-gas  pressure,  source 


temperature,  or  molecular  weight  of  the  host  gas  results  in  increases  of 


the  mean  particle  diameter.  This  is  illustrated  for  increasing  host  gas 


pressure  in  Figure  A-5-7.  These  trends  have  been  observed  experimentally 


here  as  well  as  by  previous  investigators,  but  up  to  now  have  not  been 


satisfactorily  explained 


The  published  literature  on  gas  evaporation  systems^-^  has 


emphasized  the  role  of  merging  convective  flow  of  host  gas  and  diffu- 


sional  flow  of  the  evaporating  vapor  in  explaining  the  formation  of  the 


so-called  "smoke  shell."  The  present  numerical  model  shows  clearly  that 


these  previous  explanations  are  not  correct  and  that  thermophores;  s  oh 


the  partic.es  is  responsible  for  formation  of  the  smoke  shell;  t h 


illustrated  by  Figure  A  -‘l  where  the  spatial  distribution  of 
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Figure  A-5-6.  Distribution  of  total  number  concentration  of  primary  par¬ 
ticles  in  model  evaporation  chamber  with  source  tempera¬ 
ture  and  time  as  parameters.  Width  of  "smoke  shell"  does 
not  vary  with  source  temperature.  P  =  500  dyn/cm^,  nitro¬ 
gen  host  gas.  Curve  A:  Ts=2000  K,  t  =  0.05  ms;  Curve  B: 
Tg= 1 900  K,  t=0.05  ms;  Curve  C:  Ts=2000  K,  t=0.1  ms;  Curve 
D:  Ts  =  1900  K,  t=0.1  ms. 
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Figure  A-5-7.  The  effect  of  host-gas  pressure  on  average  diameter  of 
primary  particles.  Curve  A:  ?  =  10000  dyn/cm^;  Curve  P: 
P=500  dyn/cm?. 
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Figure  a-5-8.  Distribution  of  total  number  concentration  of  primary  par¬ 
ticles  in  model  evaporation  chamber  with  and  without  par¬ 
ticle  thermophoresis.  Illustration  of  the  principal  role 
of  particle  thermophoresis  in  forming  the  "smoke  shell." 
P  =  500  dyn/em2,  T  =  200  K;  Nitrogen  host  gas.  Curve  A: 
Without  thermophoresis,  t=0.05  ms;  Curve  B:  With  thermo¬ 
phoresis,  t=0.05  ms;  Curve  C:  Without  thermophoresis, 
t=0.1  ms;  Curve  D:  With  thermophoresis,  t=0.1  ms. 
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concentration  is  shown  with  and  without  particle  thermophoresis.  In  the 
absence  of  thermophoresis,  particles  have  a  broad  spatial  distribution 
and  the  particle  free  zone  around  the  source  observed  experimentally  is 
not  found. 

Over  the  time  scale  represented  by  the  position  of  the  smoke 
shell  (up  to  approximately  0.1  ms)  coalescence  of  clusters  to  produce 
primary  particles  is  not  important,  as  verified  by  calculations  with  and 
without  the  coalescence  process.  Primary  particle  growth  clearly  occurs 
by  condensation  of  metal  vapor  atoms;  condensation  growth  is  well  known 
to  produce  narrow  particle  size  distributions.  Experimentally,  it  is 
always  observed  that  the  size  distributions  of  primary  particles  in 


99 


evaporation  chambers  are  narrow  and  that  the  primary  particles  are 
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spherical.  For  longer  times,  as  particles  move  further  away  from  the 
source,  coagulation  (not  coalescence)  becomes  important  as  evidenced  by 
the  formation  of  particle  chains. 

The  present  model  explains  present  and  previous  experimental 
observations  on  primary  particle  formation  in  metal  evaporation  systems. 
The  dependence  of  primary  particle  size  on  host-gas  pressure,  host-gas 
molecular  weight  and  source  temperature  has  been  explained.  Particle 
thermophoresis  has  been  shown  to  cause  the  formation  of  the  so-called 
"smoke  shell"  observed  surrounding  the  source.  Particle  coalescence  nas 
been  shown  not  to  be  an  important  growth  process  for  the  primary  par¬ 
ticles.  Particle  coagulation,  of  course,  for  longer  times  is  an  im¬ 
portant  growth  process  in  the  formation  of  the  chains  of  primary  par¬ 
ticles  observed  experimentally.  Particle  chain  formation  and  its 
dependence  on  the  system  parameters  will  be  the  subject  of  an  additional 
study . 


A -6.  Growth  of  ferromagnetic  particles  in  liquid  phase 


The  addition  of  a  reducing  agent  to  a  solution  containing  metal 
:ations  produces  metal  atoms  which  condense  to  form  particles.  The 
'e-action  between  metallic  salts  and  sodium  borohydride  in  the  liquid 
.nas--  to  produce  ultraf ine  me* a!  lie  particles  has  been  known  since  the 


f  . '  e  h  1  e  s  i  r.g  e  r  e  t.  a  1 . 


investigators  demonstrated  that 


.  "i  tor  .  h  •/ 1  r  i  1  r--  :  4  -  i ;  i  1  v  wit’1,  '-r  ta  in  first-row  transition 


metal  salts  in  aqueous  phase  to  yield  finely  divided  black  precipitates 
thought  to  be  metal  borides.  It  is  now  known  that  through  variationof 
ionic  concentrations,  temperature,  and  reducing  agent  ultrafine  particles 
of  a  given  mean  primary  particle  size  can  be  produced.  The  simplicity 
and  rapidity  of  the  ohydride  reduction  reaction  have  led  to  its  wide¬ 
spread  use  in  many  applications,  such  as  electroless  plating,  preparation 
of  catalysts,  and  production  of  magnetic  recording  media. 

Most  studies  on  application  of  borohydride  reduction  in 
catalyst  preparation  have  been  directed  toward  procedures  for  preparation 
of  very  highly  dispersed  mate r ia 1 .73 -7 6  Research  on  use  of  metallic 
cation  reduction  for  preparation  of  magnetic  recording  coatings  has  been 
focused  mainly  on  the  magnetic  properties  and  oxidation  resistance  of  the 
particles  produced.  ??  However,  in  previous  research  the  dependence  of 

particle  properties  on  the  various  experimental  parameters,  such  as 
cation  and  anion  nature  and  concentration,  solvent,  stabilizers,  addi¬ 
tives,  and  temperature  have  not  been  investigated  systematically. 

Work  described  here  was  designed  to  develop  fundamental  in¬ 
formation  on  particle  formation  and  growth  in  the  reduction  of  metal  It' 
cations  by  borohydride.  Of  particular  interest  were  the  growth 
mechanisms  of  the  primary  particles,  the  effect  of  the  experimental 
parameters  noted  above  on  the  properties  of  the  primary  particles,  ,ir.d 
the  role  of  primary  particle  properties  on  subsequent  growth  ..  f 
aggregates.  First,  the  experimental  procedures  used  in  this  invest  v  »- 
tier,  are  described.  This  is  followed  by  a  summary  and  analysis  v  ‘  : 


! 
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experimental  results.  Finally,  a  theoretical  model  is  developed  and  used 


to  interpret  qualitatively  these  results. 


a .  Experiment  . 


The  reaction  between  the  borohydride  ion  and  metallic  aniens  is 
initially  rapid.  Therefore,  in  order  to  obtain  reproducible  results 
independent  of  the  fluid  dynamical  state,  a  mixing  junction  was  designed 
similar  to  that  used  for  stopped-flow  reactors  for  the  study  of  fact 
reactions.®1  This  junction  was  in  the  form  of  a  "tee"  formed  from  brass 
tubing  with  an  inside  diameter  of  1  mm.  Mixing  of  the  two  reactant 
streams  could  be  effected  on  time  scales  down  to  a  few  microseconds.  The 
flow  from  the  mixing  junction  was  directed  to  an  agitated  ultrasonic  bath 
where  the  slower  final  phase  of  particle  growth  took  place.  External 
magnet5ic  fields  could  be  applied  bot  to  the  "tee"  reactor  and  to  the 
ultrasonic  bath.  The  particle  suspension  in  the  bath  was  washed  with 
deionized  water  to  remove  the  soluble  reaction  by-products  and  the  solid 
particles  were  than  finally  dispersed  in  acetone;  all  these  steps  were 
carried  out  in  the  absence  of  oxygen.  This  method  of  particle  formation 
was  the  principal  one  used  in  these  studies.  However,  an  alternative 
method  for  bringing  about  reaction  with  rapid  mixing  was  also  employed. 
In  this  method,  reactants  were  aerosolized  in  two  separate  nebulizers. 


The  aerosol  jets  from  these  nebulizers  were  directed  towards  each  ot ner 
As  a  result,  droplet  collision  and  coalescence  occurred  so  that  to 
reduction  reaction  was  confined  to  the  individual  droplets.  ■ 
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droplets  were  of  the  order  of  1  p  in  diameter  so  that  the  diffusional 
time  scale  was  less  than  the  order  of  1  msec. 

Particles  produced  in  the  "tee"  reactor  system  were  studied  by 
a  number  of  techniques  including  transmission  electron  microscopy  (TEM). 
The  first  step  here  was  to  sample  the  particles  collected  from  the 
reactor  as  a  dispersion  in  acetone.  The  most  obvious  method  was  to  place 
a  drop  of  the  acetone  dispersion  onto  a  grid;  however,  as  the  acetone 
evaporates  the  particles  in  the  droplet  tend  to  agglomerate  on  the  grid 
and  this  makes  it  difficult  to  obtain  satisfactory  micrographs  for 
analysis.  Another  method  involved  placing  a  drop  of  the  acetone  disper¬ 
sion  onto  a  glass  microscope  slide,  spreading  the  drop  over  the  3lide  by 
shearing  and  then  sampling  the  sheared  deposit;  however,  it  was  found 
that  the  shearing  action  tended  to  break  the  original  chains  formed  in 
the  reaction  system.  The  most  successful  method,  and  the  one  used, 
involved  nebulizing  the  acetone  dispersion  as  approx ima tely  10  m 
droplets  which  were  then  collected  on  the  TEM  grids  by  direct  impaction; 
this  method  was  found  to  minimize  alteration  of  particle  characteristics 
noted  with  the  other  methods  and  could  be  applied  successfully  to  any 
concentration  of  the  acetone  dispersion.  Grids  with  deposited  particles 
were  stored  in  nonoxidizing  conditions  and  observed  using  a  JEOL  JEM- 
200CX  transmission  electron  microscope. 

In  addition  to  the  study  by  TEM,  the  saturation  magnetiziation 
(SM)  of  the  particles  was  also  measured.  In  preparation  for  the  SM 
determination,  the  acetone  dispersion  was  dried  under  a  nitrogen  atmos¬ 
phere  at  room  temperature,  although  in  some  experiments  drying  was 
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carried  out  at  elevated  temperatures  to  investigate  the  effects  of  ther¬ 


mal  annealing  on  SM.  In  this  manner,  particle  samples  of  approximately 
0.02  and  0.05  grams  were  placed  in  a  sample  holder  and  used  in  measure¬ 
ment  of  SM  in  a  vibrating  sample  mangetometer  (Model  FM-1  Princeton 
Applied  Research).  This  instrument  was  initially  calibrated  with  a 
standard  nickel  sample.  Reproducibility  of  SM  data  with  these  techniques 
was  very  good.  A  portion  of  the  sample  used  for  the  SM  determination  was 
at  the  same  time  used  to  get  information  on  chemical  composition  of  the 
particles.  The  bulk  and  surface  chemical  compositions  were  determined  by 
Inductively  Coupled  Plasma  Atomic  Emission  Spectroscopy  (ICP-AES)  and  X- 
ray  Photoelectron  Spectroscopy  (XPS). 

The  reactants  were  all  reagent  grade  chemicals  and  were  used  in 
the  experiments  without  additional  purification.  Solutions  of  ferrous 
chloride,  ferrous  sulfate,  cobaltous  chloride,  and  nickelous  chloride 
we'”-'!  used  as  the  sources  of  the  metallic  cations.  The  reducing  agent  was 
sodium  borohydride;  the  standard  reducing  solution  used  in  the  experi¬ 
ments  was  0.1  M  to  which  was  added  0.1  wtH  sodium  hydroxide  to  maintain  a 
basic  solution  and  prevent  hydrolysis  of  the  sodium  borohydride.  In  the 
experiments,  cation  and  anion  concentrations  were  varied.  In  order  to 
study  the  effects  of  additives,  chromium  chloride  was  used  as  an  anti- 
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xiiant  and  gelatin  and  sodium  citrate  were  used  as  stabilizers.  Neo- 


iun  chloride  was  also  used  in  separate  studies  of  supermagnetic 


.als.  In  most  experiments,  deionized  water  was  the  solvent;  how- 


me  runs  were  carried  out  in  which  mixtures  of  deionized  water  and 


vrade  ethyl  alcohol  were  the  solvents. 
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As  noted  above,  two  methods  were  used  for  carrying  out  the 
borohydride  reduction  reaction  to  produce  particles.  One  of  these  in¬ 
volved  nebulization  of  the  reactants  in  opposed  jets  so  as  to  bring  about 
subsequent  collision  and  coalescence  of  the  reactant  droplets.  This  mode 
of  reaction  produced  very  small  primary  particles  (approximately  50  A  in 
diameter)  which  were  confined  on  the  TEM  grids  to  the  areas  of  the 
evaporated  acetone  droplets  which  were  around  0.2  to  0.5  pm  in  diameter. 
The  second  method  used  the  "tee"  reactor  system.  With  this  system,  the 
primary  particle  sizes  varied  from  approximately  100  to  1500  A  in 
diameter,  depending  on  the  experimental  conditions.  In  most  of  the 
experimental  runs  using  the  "tee"  reactor  system,  the  primary  particles 
formed  essentially  linear  chains  with  some  branches.  External  magnetic 
fields  were  necessary  to  produce  straight  chains  without  branching.  The 
primary  particles  in  a  given  chain  were  found  to  have  narrow  size  distri¬ 
butions  . 

Table  A -6-1  shows  the  bulk  weight  percent  of  the  elements  in 
particles  produced  by  the  "tee"  reactor  system  as  determined  gy  ICP-AES 
analysis;  for  this  analysis  all  the  samples  were  obtained  from  the 
"standard"  reducing  solution  and  0.1  M  metallic  cation  solutions.  The 
alloy  samples  were  prepared  from  equimolar  mixtures  of  the  cations  an  J 
used  the  "standard"  borohydride  solution  with  0.1  M  metallic  cation 
solutions.  Nitric  acid  was  used  as  solvent  for  the  particle  sample  for 
ICP-AES  analysis;  the  ratio  of  the  mass  from  the  analysis  to  the  t„‘  i'. 
mas?  ~,r  *  r;e  original  sample  shows  that  some  minor  components  of  t  r.e 


Table  A  —6  — 1 :  Bulk  Coraposi tion  of  Particles 


Sample 

Number 

Sample 

I.D. 

Na 

Fe 

Co 

wt .% 

Ni 

B 

Total 

% 

Me  tal/B 
( Atomic) 

I 

( Fe , Co) ?b 

<0.02 

34.71 

53.59 

0. 1 1 

9.37 

97.78 

1.82 

II 

Co2S 

<0.03 

0.02 

83. 16 

0.18 

8.84 

92.20 

1.72 

III 

Ni2B 

<0.02 

0.04 

<0.01 

73.65 

3.14 

81  .83 

1.66 

IV 

(Fe,Ni)2B 

<0.02 

32.35 

<0.01 

53.04 

9.03 

94.  42 

1.83 

V 

Fe2B 

<0.02 

89.09 

<0.01 

0.02* 

7.96 

97.07 

2. 16 

*  Close 

to  the  detection 

limit 

partic , 

_es  were 

insoluble.  Atomic 

ratios 

of  metal  to 

boron  are 

approximately  2 

for  all 

sam  pies. 

It  can  al so 

be  seen 

that  in 

the  alloy 

samples,  the  cobalt  or  nickel  content  is  always  higher  than  that  of  iron 
even  though  the  reactant  cation  concentrations  were  equimolar.  Analysis 
of  the  samples  by  XPS  for  surface  composition  showed  that  the  ratio  of 
metal  to  boron  on  the  surface  was  one  or  smaller.  This  indicates  that 
the  boron  was  deposited  mainly  on  the  surface  of  the  metal  particles  and 
that  the  particles  are  probably  not  metal  borides  but  are  an  admixture  of 
boron  and  metal.  In  the  absence  of  thermal  treatment,  all  the  particles 
except  those  prepared  at  very  low  metallic  cation  concentration  did  not 
show  any  crystalline  diffraction  pattern  with  TEM.  All  the  dried  parti¬ 
cles  were  pyrophoric  in  air.  Also  all  the  particles  were  ferromagnetic 
except  particles  formed  from  nickel  cations;  these  particles  became 
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ferromagnetic  when  heated  for  20  min.  in  nitrogen  at  around  100°C. 


In  all  of  the  studies  of  the  production  of  particles  by 


borohydride  reduction,  it  was  found  that  variation  of  the  anions  (sulfate 


and  chloride)  associated  with  the  metallic  cations  did  not  affect  the 


particle  characteristics.  However,  variation  of  the  cations  (ferrous, 


cobaltous  and  nickelous)  led  to  large  changes  in  primary  particle  size, 


chain  forming  ability  of  the  primary  particles,  and  saturation  magnetiza¬ 


tion  of  the  particles.  In  the  study  of  the  effect  of  variation  of  the 


metallic  cation,  various  amounts  of  cobaltous  chloride  or  nickelous 


chloride  were  added  to  ferrous  chloride  while  maintaining  the  same  molar 


concentration  of  total  metallic  cation  (0.1  M).  As  shown  in  Figure  A-6- 


1,  significant  changes  occur  in  the  mean  primary  particle  diameter  as  the 


compositions  are  varied.  Chain  length  was  found  to  diminish  gradually 


with  increasing  content  of  cobaltous  or  nickelous  ion.  Particles  formed 


from  mixtures  of  cobaltous  and  nickelous  ions  in  the  absence  of  iron  did 


not  form  chains  even  with  magnetic  fields  up  to  1200  Gauss.  For  parti¬ 


cles  formed  from  single  component  metallic  ion  solutions,  the  saturation 


magnetization  (5M)  values  decreased  in  the  sequence  Fe  (130  emu/g)  >  Co 


(30  emu/g)  >  Ni  (0.4  emu/g).  The  variation  in  SM  values  of  the  Fe  +  Co 


alloy  particles  ranged  from  1 30  to  30  emu/g,  decreasing  as  the  fraction 


of  Fe  decreased. 


Figures  A -6 -2  and  A-6-3  show  the  effect  of  variation  of  ferrous 


cation  concentration  on  the  mean  primary  particle  size  and  saturation 


magnetization  (SM)  respectively.  Without  any  of  the  two  additives, 
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Effect  of  ferrous  ion  concentration  on  mean  primary 
particle  size. 
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Figure  A-6-3.  Effect  of  ferrous  ion  concentration  on  saturation  magneti¬ 
zation  of  particles. 
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sodium  hydroxide  and  chromium  (III)  chloride,  the  mean  primary  particle 
size  is  not  affected  by  the  variation  of  ferrous  cation  concentration. 
With  the  "standard"  reducing  solution,  the  mean  primary  particle  size 
shows  great  variation  with  ferrous  ion  concentraiton.  In  this  case,  as 
shown  in  Figure  A-6-2,  above  a  critical  cation  concentration,  the  mean 
primary  particle  size  remains  nearly  constant,  while  below  that  concen¬ 
tration,  the  particle  size  increases  rapidly.  These  large  particles  tend 
not  to  form  chains  in  the  absence  of  an  external  magnetic  field,  indi¬ 
cating  that  these  particles  are  of  multidomain  size.  The  transition 
noted  takes  place  at  higher  ferrous  ion  concentrations  and  is  steeper 


when  chromium  (III)  chloride  (0.002M)  i3  added  to  the  ferrous  chloride 
reactant  solution.  In  the  absence  of  the  two  additives,  the  saturation 
magnetization  (SM)  did  not  vary  with  ferrous  ion  concentration.  Upon 
addition  of  NaOH,  the  SM  drops  rapidly  at  approximately  the  same  critical 
ferrous  cation  concentration  observed  for  variation  of  mean  primary 
particle  size.  Below  the  critical  ferrous  ion  concentration,  where 
larger  particles  are  produced,  it  was  found  that  the  mean  size  and  SM 
values  were  highly  dependent  on  residence  time  in  the  reactor  system. 
Below  the  critical  ferrous  cation  concentration,  a  gel-like  green  non¬ 
magnetic  material  first  appeared  and  then  slowly  changed  into  black  or 
dark  brown  ferromagnetic  particles.  The  lifetime  for  the  green  gel 
increased  as  the  ferrous  cation  concentration  was  lowered.  In  this  same 
regime,  at  a  constant  ferrous  ion  concentration  in  the  reactant  solution, 
a  higher  initial  NaOH  concentration  led  to  longer  lifetimes  of  the  green 
gel.  The  existence  of  this  phenomenon  was  also  corroborated  by  electron 
micrographs  and  SM  measurements.  Precipitate  obtained  for  small 
residence  times  in  the  reactor  system  (approximately  1  min.)  appeared  in 
electron  micrographs  as  sheet-like  nets  with  very  small  spherical  parti¬ 
cles  in  some  of  the  nets;  this  is  shown  in  Figure  A-6-4.  Subsequently, 
for  larger  residence  times,  these  small  particles  appeared  to  grow  at  the 
expense  of  the  net  material  while  the  SM  values  increased  and  finally 
reached  the  values  (about  130  emu/g)  characteristic  of  particles  formed 
above  the  critical  value  of  ferrous  ion  concentration.  All  the  data 
shown  in  Figures  A-6-2  and  A-6-3  were  obtained  from  samples  whose 
residence  time  in  the  reactor  system  was  5  minutes. 
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Figure  A-6-4.  Primary  particle  formation  below  critical  ferrous  ion 
concentration.  Standard  reducing  solution  and  0.01  MFe++. 
Ferromagnetic  particles  are  shown  in  growth  sites  in 
hydroxide  get  sheet-like  nets.  Residence  time  in  reactor 
system  was  5  min. 

Since  particles  are  produced  as  a  result  of  chemical  reaction, 
it  might  be  expected  that  the  reaction  temperature  would  be  an  important 
parameter  in  the  formation  of  primary  particles.  In  the  study  of  the 
effect  of  temperature,  the  entire  reactor  system  was  maintained  at  a 
fixed  temperature,  which  could  be  varied  from  5  to50°C.  Figure  A-6-5 
shows  the  variation  of  mean  primary  particle  size  with  temperature  for 
0.1  M  ferrous  ion  initial  reactant  concentration  and  the  standard  boro- 
hydride  solution  (0.1  M  with  0.1  wt.  %  NaOH).  As  the  reaction  tempera¬ 
ture  increased  the  mean  primary  particle  size  decreased,  being  at  50°C 
less  than  half  that  at  5°C.  At  the  same  time  the  chaining  tendency  and 
SM  values  were  not  affected  by  the  temperature  of  the  reactor  system. 
Since  the  reaction  was  exothermic  and  the  reactor  system  did  not  permit 
maintenance  of  the  system  free  of  temperature  gradients,  the  actual 
reaction  temperatures  were  slightly  greater  than  the  nominal  values  shown 
for  Figure  A-6-5. 


TEMPERATURE  (  K) 

Figure  A-6-5.  Effect  of  reactor  temperature  on  mean  primary  particle 
size.  Ferrous  ion  (0.1  M)  reactant  solution  in  standard 
reducing  solution. 

Additives  in  the  borohydride  reduction  reaction  have  been  used 
by  others  to  modify  particle  characteristics.  A  limited  investigation 
was  carried  out  on  the  effects  of  various  additives  on  particle  charac¬ 
teristics.  When  a  small  amount  of  gelatin  (0.5  wt.  %)  is  added  to  0.1  M 
ferrous  ion  reactant  with  the  standard  reducing  solution,  random  clusters 
are  formed  instead  of  chains.  Similar  behavior  was  found  with  addition 
of  sodium  citrate  to  the  ferrous  ion  reactant.  The  addition  of  Cr+^  at 
0.002  M  in  0.1  M  ferrous  ion  reactant  solution  led  to  increases  in  the  SM 
values  of  the  resultant  particles.  However,  when  the  chromium  ion  con¬ 
centration  was  increased  to  0.005  M  in  the  0.1  M  ferrous  ion  reactant 
solution,  the  chaining  tendency  of  the  resultant  particles  was  sig¬ 
nificantly  reduced  as  was  the  SM  value.  Addition  of  chromium  was  also 
found  to  prevent  oxidation  of  the  particles.  Without  chromium,  from 


electron  micrographs  of  the  particles  a  noticeable  film  was  observed  to 
form  on  particle  samples  after  an  exposure  time  to  air  of  one  day.  This 


film  had  an  amorphous  structure  from  the  electron  diffraction  pattern, 


and  based  on  Olsen's  findings®^,  was  believed  to  be  iron  oxide.  When 
chromium  was  added  to  the  ferrous  ion  reactant  solution,  no  film  was 
observed  to  form  on  the  particles  produced  even  after  exposures  of 
samples  to  air  for  as  long  as  six  weeks. 

If  borohydride  reduction  of  ferromagnetic  cations  was  carried 
out  in  solvent  mixtures  of  water  and  ethyl  alcohol,  it  was  found  that 
mean  diameters  of  the  primary  particles  decreased  as  the  fraction  of 
ethyl  alcohol  in  the  solvent  increased.  There  is  currently  no  explana¬ 
tion  for  this  solvent  effect.  However,  the  SM  values  of  particles  pro¬ 
duced  with  from  0  to  50  volume  %  ethyl  alcohol  as  solvent  remained 
constant . 

External  magnetic  fields  greater  than  approximately  500  Gauss 
applied  to  the  reactor  system  produced  highly  linear  chains  with  little 
or  not  branching.  The  external  field  was  varied  over  the  range  0  to 
12500  Gauss,  but  no  appreciable  changes  in  primary  particles  or  their 
chain  formation  were  observed  after  the  field  exceeded  500  Gauss. 

The  saturation  magnetization  (SM)  values  for  particles  produced 
from  0.1  M  ferrous  and  0.1  M  cobaltous  reactant  solutions  were  130  and  30 
erau/g  respectively.  These  values  are  small  compared  to  the  values  for 
pure  bulk  iron  and  cobalt.  The  lowered  SM  values  are  attributed  here  to 
the  admixture  of  boron  in  the  particles  and  to  the  small  particle  size. 
The  small  size  of  the  cobalt  and  nickel  particles  imply  a  larger  Brownian 
diffusivity  than  that  for  the  larger  iron  particles.  This  higher 
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diffusivity,  together  with  lower  SH  values,  are  believed  to  be 
responsible  for  the  reduction  in  chain  formation  for  the  nickel  and 
cobalt  particles.  Reduction  in  chain  formation  was  also  observed  for  the 
larger  iron  particles  formed  from  reactant  solutions  with  ferrous  ion 
concentrations  below  the  critical  value  shown  in  Figure  A  -6~2.  These 
larger  particles  did  not  form  chains  even  though  their  3d  values  were  the 
same  as  the  iron  particles  produced  from  higher  concentration  ferrous  ion 
reactant  solution  ( 1 30  f>nu/g).  This  failure  to  form  chains  is  attributed 
to  the  multidomain  nature  of  these  large  particles. 

The  explanation  for  the  effect  of  temperature  on  mean  primary 
particle  size  is  straightforward.  Higher  reaction  temperatures  are 
believed  to  result  in  higher  reaction  rates  which  in  turn  result  in 
higher  rates  of  production  of  metal  atoms  available  for  nucleation. 
Therefore,  more  nuclei  are  formed  at  higher  temperatures  and,  as  a 
result,  fewer  metal  atoms  are  available  for  condensation  on  the  nuclei. 
This  picture  is  supported  also  by  the  variation  of  mean  primary  particle 
size  with  the  cation  and  the  chemical  analysis  of  t  associated  par~ 
tides.  As  Table  A-t  -1  shows,  the  atomic  ratios  of  Co/.e  li/Fe  are 

always  greater  than  one,  even  though  the  reactions  occurred  reactant 

solutions  that  were  equimolar  in  the  two  cations.  The  higher  proportion 
of  nickel  and  cobalt  in  the  particles  suggests  that  iron  cations  are 
reduced  more  slowly  than  the  cobalt  and  nickel  cations.  In  addition,  it 
may  be  noted  that  pure  nickel  and  cobalt  cation  reactant  solutions 
yielded  smaller  m-an  primary  particle  sizes  than  *h~se  found  for  the 


ferrous  ion;  this  confirms  the  view  that  higher  react: or.  ntes  result  in 
smaller  mean  primary  particle  sizes. 

From  the  experiments  on  variation  of  concentration  of  ferrous 
ion  in  reactant  solution,  summarized  in  Figure  A -6-?,  it  appears  that 
there  are  two  quite  separate  mechanisms  for  formation  of  primary 
particles  by  borohydride  reduction.  Below  the  critical  ferrous  ion 
concentration  indicated  in  Figure  A-6-2,  the  first  stc  r  in  formation  of 
primary  particles  is  the  formation  of  a  green  gel,  probably  an  hydroxide, 
which  is  nonmagnetic  and  has  the  appearance  of  a  sheet-like  net  in 

Q  o 

transmission  electron  micrographs  (see  Figure  A -6-4'.  Paltz  et  al.  3 
observed  this  sheet-like  net  in  some  of  their  sampler  and  identified  it 

Oil 

as  FeOOH  by  electron  diffraction.  However,  Detourney  et.  al  as  well  as 
Bernal  et  alw5  observed  this  green  gel  in  an  electrolytic  bath  for  pro¬ 
duction  of  iron  particles  and  found  it  to  be  a  more  complex  material  than 
FeOOH  which  slowly  changes  in  FeOOH.  On  the  basis  of  studies  of  this 
process  at  various  reaction  times,  it  is  suggested  here  that  below  the 
critical  cation  concentration  the  hydroxide  complex  forms  first,  and  that 
possibly  active  sites  in  the  hydroxide  gel  become  the  -enters  for 
subsequent  reduction  to  iron  to  form  primary  particles.  Above  *he  criti¬ 
cal  ferrous  ion  concentration,  the  sheet-like  nets  of  ^reer.  hydroxide  gel 
are  not  observed;  instead  primary  particles  are  form"-,  iireotiy  in  the 
reaction.  Therefore,  at  the  higher  ferrous  cation  concentrations,  it  is 
inferred  that  the  rate  of  the  reduction  reaction  is  h’.zner  than  that  of 


the  hydroxide  complex  reaction. 


As  noted  previously,  in  all  the  reaction  systems  the  primary 
particles  collided  to  form  aggregates.  The  details  of  the  various  aggre¬ 
gation  processes  observed  in  these  experiments  are  not  well  understood. 
Therefore,  a  measure  of  these  aggregate  structures  may  be  of  aid  in 
studying  the  aggregation  mechanisms.  Such  a  measure  is  the  Hausdorf  or 
fractal  dimension®®  of  an  aggregate  which  has  been  applied  in  investiga¬ 
tion  of  a  variety  of  aggregation  processes,  such  as  formation  of  branched 
polymers,  the  sol-gel  transition  formation  of  soot  particles, etc.  To 
date,  most  of  the  investigations  of  this  concept®?'®®  have  involved 
computer  simulations  of  model  growth  processes  which  give  scale  invariant 
aggregate  structures  that  can  be  described  as  fractals.  In  the  regime  of 
scale  invariance,  the  mass  of  a  cluster  or  aggregate  is  predicted  to 
scale  as  where  L  is  a  measure  of  the  cluster  size  and  D  is  the 

Hausdorf  or  fractal  dimension.  From  the  present  experimental  electron 
micrographs,  estimates  of  the  Hausdorf  dimension  of  many  aggregates  were 
obtained  by  correlating  the  number  of  particles  in  any  aggregate  to  a 
characteristic  length  of  the  aggregate  —  in  this  case  the  geometric  mean 
of  the  longest  linear  dimension  of  the  aggregate  and  the  length  perpen¬ 
dicular  to  that  axis.  The  correlations  were  carried  out  for  aggregates 
produced  by  the  standard  solution  borohydride  reduction  in  the  reactor 
system  of  pure  0.1  M  ferrous  ion  and  pure  0.1  M  cobaltous  ion  in  the 
absence  of  any  external  magnetic  field.  These  correlations  are  shown  in 
Figures  A-6-6  and  A-6-7,  and  the  data  are  well  described  by  a  power  law 
as  would  be  expected  for  fractal  objects.  The  fractal  dimension  for  the 
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6-6.  Fractal  geometry  of  iron  particle  aggregates  prepared 
from  0.1  M  Fe4+  reactant  solution  in  standard  reducing 
solution.  Fractal  dimension  is  1.5*1. 


Figure  A -6 -7. 


Fractal  geometry  of  cobalt  particle  aggregates  prepared 
from  0.1  M  Co'M’  reactant  solution  in  standard  reducing 
solution.  Fractal  dimension  is  1.72. 


aggregates  of  the  iron  particles  was  1.5*4  while  that  for  cobalt  particles 
was  1.72.  It  should  be  noted  that  the  Hausdorf  dimension  of  linear  chain 
particles  would  be  1.0,  and  this  is,  of  course,  observed  for  the  iron 
particles  under  external  magnetic  fields  greater  than  500  Gauss.  It  is 
clear  from  Figures  A-6-6  and  A-6-7  that  the  iron  particle  aggregates, 
with  their  larger  saturation  magnetization  values,  exhibit  a  greater 
degree  of  "linearity"  than  do  the  cobalt  particle  aggregates.  The 
Hausdorf  dimension  in  this  case  may  be  a  useful  measure  of  the  trend 
toward  linearity  of  ferromagnetic  particle  aggregates  and  the  magnetic 
properties  of  the  primary  particles. 

b .  Theory 

In  this  paper,  only  the  processes  leading  to  the  formation  of 
the  primary  particles  are  considered  in  the  context  of  a  theoretical 
model;  this  model  is  intended  to  provide  only  a  qualitative 
interpretation  of  the  experimental  results  presented  above  for  particle 
formation  above  the  critical  cation  concentration.  Additional  work  will 
be  reported  concerning  models  for  chaining  and  aggregation  of  these 
primary  particles. 

In  the  borohydride  reduction  of  metallic  cations,  the  cations 
are  red -iced  to  metal  atoms.  In  the  liquid  phase,  these  atoms  have  such 
low  solubility  that  nucleation  occurs  with  effectively  zero  activation 
energy  and  the  dimer  is  the  smallest  stable  cluster.  As  reaction 
proceeds,  clusters  may  grow  by  accretion  of  atoms  and  other  clusters.  A 
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central  question  m  consideration  of  the  formation  of  primary  particles 
is  the  relative  contribution  to  the  growth  process  of  condensation  of 
atoms  and  collision  with  other  clusters  (i.e.  coagulation) .  In  aqueous 
phase  in  the  absence  of  external  magnetic  fields,  collisions  between 
clusters  are  subject  to  three  principal  interaction  energies  between  a 
pair  of  clusters.  These  are  van  der  Wallas,  hydrodynamic  and  double 
layer  energies.  The  van  der  Waals  forces  arise  from  the  synchronized 
dipoles  created  by  fluctuations  in  the  electron  clouds  of  the  constituent 
atoms  composing  the  clusters.  Hydrodynamic  forces  arise  from  the 
interaction  between  a  pair  of  particles  with  relativemotion  in  a  viscous 
medium.  Double  layer  forces  are  due  to  the  interactions  between  the 
electrical  double  layers  surrounding  each  particle.  Of  these  three 
forces,  only  the  van  der  Waals  force  i3  attractive.  Analysis  of  these 
forces  is  the  starting  point  for  consideration  of  the  stability  of  col¬ 
loidal  dispersions.^, 90  generai  picture  is  that  van  der  Waals 

attraction  predominates  at  small  and  large  interparticle  distances;  at 
intermediate  distances,  the  electrical  double  layer  forces  may  be 
dominant,  depending  on  the  relative  magnitudes  of  these  two  forces.  The 
hydrodynamic  force  is  normally  incorporated  by  a  redefinition  of  the 
Brownian  diffusion  coefficients  of  a  pair  of  particles.  If  the  electri¬ 
cal  double  layer  forces  are  dominant,  the  total  energy  as  a  function  of 

interparticle  distance  displays  a  positive  maximum  and  the  disperse 
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system  is  stable;  the  system  exhibits  what  is  termed  slow  coagulation. 

If  the  double  layer  forces  are  not  dominant  at  intermediate  distances, 
the  disperse  system  is  unstable  and  exhibits  what  has  been  termed  rapid 


coagulation.  These  considerations  are  of  obvious  importance  in  the  model 
development  presented  below.  Intercluster  magnetic  forces  are  not  con¬ 
sidered  since  the  clusters  are  below  the  critical  domain  size. 

The  physical  system  to  be  described  consists  initially  of  a 
spatially  homogeneous  mixture  of  borohydride  ions  and  metallic  cations  in 
basic  aqueous  phase;  the  cation  concentration  is  taken  to  be  greater  than 
the  critical  value  so  that  the  hydroxide  mechanism  is  not  included.  The 
ensuing  reduction  reaction  constitutes  the  source  of  metal  atoms  that 
undergo  a  zero  activation  energy  nucleation  process  to  form  clusters. 
These  clusters  grow  by  accretion  of  additional  atoms  and  possibly  by 
accretion  of  other  clusters.  Conservation  equations  can  now  be  applied 
which  describe  the  formation  and  growth  of  the  primary  particles. 

The  conservation  equation  for  metal  and  boron  atoms  has  the 

form: 

ac/dt  =  RXN  -  (1/ra^  J*x*  v  ( x,c,t)  n(  x,t)  dx  -  2R  (A-6-1) 

where  c  is  the  number  concentration  of  metal  and  boron  atoms.  m-|  is  the 
effective  atomic  mass.  It  is  assumed  that  metal  and  boron  are  produced 
by  the  reaction  and  are  both  incorporated  in  the  particle  formation 
process.  RXN  is  the  rate  of  production  and  metal  rate  coefficients.  The 
second  term  on  the  right  hand  side  represents  loss  of  atoms  by  attachment 
to  existing  particles,  being  the  growth  rate  of  a  particle  mass  x  and 
n(x,t.)  the  mass  distribution  function  of  particles.  R  represents  loss  of 
atoms  by  nucleation  of  clusters. 


There  have  been  no  rate  studies  which  have  established 


quantitative  values  for  RXN  for  the  borohydride  reduction  of  metallic 
cations.  Ivanov  et  al.92  have  carried  out  a  mechanistic  study  for  the 
reduction  of  nickelous  cations  by  borohydride  in  basic  aqueous  solution. 
Their  mechanism  is  derived  from  mass  spectroscopic  studies  of  the  re¬ 
action  in  deuterium  oxide.  They  proposed  the  following  mechanism  which 
is  used  in  the  present  analysis: 

BH^  +  4H0H  +  2NI2  +  =  2Ni  +  2H2  +  4H+  +  B02  +  2H20 
BH4  +  H+  =  B  +  2.5H2 

BH4  +  4H0H  =  B02  +  4H2  +  2H20  (A-6-2) 

from  this  mechanism  the  following  rate  expression  is  used  for  RXN: 

RXN  =  6k i  [BH4]  [Ni2+]2  ( A-6-3) 

The  rate  coefficient  k -|  is  an  unknown  parameter  which  is  varied  in  this 
analysis . 

In  liquid  phase,  particles  grow  in  the  continuum  regime.  The 
growth  rate  for  this  regime  has  the  form: 

.  =  BkTm^7^1  73c/3u  :  x  1  / 3 c  (A-6-4) 

The  nucleation  term  has  the  form: 

ri  =  R*  '  (  x  -  x* )  (  A-6-5) 


where  x*  is  the  mass  of  the  smallest  cluster,  a  dimer  in  this  case  and  R# 


is  the  coagulation  rate.  ■'  is  the  delta  function.  The  quantities  in 
equation  (A-6-5)  are  defined: 

x*  =  2m  1 

R«  =  bnc2,  bn  =  8k T/ 3 

where  kT  is  the  thermal  energy  and  p  the  viscosity  coefficient. 

The  particle  evolution  equation  is  given  by: 

5n(x,t)  1  x 

■  -  -  I  E  b(  x-x’  ,n)  n(  x-x’  ,t)  n(  x'  ,t)  dx' 

5t  2  x* 


-  n(x,t)  /  E  b(  x' ,x) n( x* ,t) dx' 


—  [^(x.c.t)  n(  x,t)  ]  +  R 
3x 


( A-6-6) 


The  first  two  terms  on  the  right  hand  side  account  for  coagulation  of 
particles  and  the  rate  coefficient  is  given  by  the  continuum  regime 
expression. 93  g  i3  the  collision  efficiency,  (0  <  E  <  1).  b  is  the 
Smoluchowski  collision  rate.  The  third  term  represents  particle  growth 
by  condensation. 

For  a  qualitative  picture  of  the  formation  process  of  primary 
particles  it  is  sufficient  to  obtain  an  approximate  solution  of  these 
equations  through  a  moment  technique  for  the  particle  phase.  The  k-th 
moment  of  the  particle  size  distribution  function  is: 


Mu(t) 


xkn( x,t) dx 


( A -6-7) 


The  set  of  moment  equations  is  closed  through  the  assumption  that  the 
particles  evolve  through  time  with  a  lognormal  size  distribution.  The  k- 
th  moment  of  the  size  distribution  is  therefore  given  in  terms  of  the 
lognormal  parameters,  Xg  the  geometric  mass  and  c  the  geometric  standard 
deviation : 

Mk(t)  =  M0xgkexp(  9k2ln^o/2)  (A-6-8) 

Any  fractional  moment  of  the  distribution  can  be  expressed  in  terms  of 
the  three  moments  Mq,  m-|  and  M2,  and  the  moment  equations  are  therefore 
closed.  The  mean  particle  diameter,  Dm,  can  be  expressed  in  terms  of 
these  moments: 

Dm3  =  6M0-3/2Mi2M2-1/2/  p7T 

The  atom  conservation  equation  and  the  moment  equations 
analyzed  here  have  the  form: 


dc/dt  = 

RXN  -  Mj/^icmi  -  biic2 

( A-6-9) 

dM0/dt  = 

bi(M02  +  ^ T  / 3M _ -j  / ^ )  +  &i 1 °2/2 

(A-6-10) 

dM^dt  = 

M1/3^1c  +  m1b1 1°2 

( A-6-1 1) 

dM2/dt  = 

2M4/3>^ic  +  2b](Mi2  +  M4/3M2/3)  +  2mi2b-|iC2 

(A-6-12) 

The  set  of  equations  ( A  —6  —9) — ( A-6-12)  can  now  be  solved  numeri¬ 
cally  to  give  a  qualitative  picture  of  the  growth  processes  of  primary 
particles.  Solutions  have  been  obtained  for  two  cases  which  correspond 
to  slow  coagulation,  E=0  (designated  as  Case  1)  and  rapid  coagulation, 
E=1  (designated  as  Case  2).  The  initial  conditions  used  were: 


c(  0)  r-  Mq  r  O'  ^  M  i  ( 0)  =  MofO)  =  0;  Fe2  +  =  BH^  =  0.1M. 


Quantitative  information  on  the  rate  coefficients  for  borohydride 
reduction  of  metallic  cations  was  not  available.  Certain  qualitative 
features  of  primary  particle  growth  can,  however,  be  examined.  First, 
there  is  the  observation  that  higher  reduction  rates,  whether  brought 
about  by  higher  temperatures  or  change  of  cation  from  iron  to  cobalt  or 
nickel,  lead  to  smaller  mean  diameters  of  primary  particles.  Second,  it 
is  uniformly  observed  that  primary  particles  produced  under  given  re¬ 
action  conditions  have  a  narrow  particle  size  distribution.  An  accept¬ 
able  mechanism  for  primary  particle  formation  and  growth  must  account  for 
these  observations. 

The  moment  equations,  ( A— 6— 9)  —  ( A-6-1 2) ,  were  solved  for  Cases  1 
and  2  for  two  orders  of  magnitude  variation  in  the  rate  coefficient,  k-j. 
Since  quantitative  rate  data  were  not  available,  the  values  are  only 
esitmates  chosen  to  give  a  reasonable  range  of  particle  production  rates. 
For  Case  2,  E  =  1,  solution  of  equations  ( A  —6  —9)  —  ( A  —6  —  1  2)  result  in  vari¬ 
ations  of  primary  geometric  mean  particle  diameter  with  reaction  rate 
that  are  exactly  opposite  to  experimentally  observed  variation;  in  the 
simulations,  mean  diameters  increase  with  increasing  reaction  rates. 
Also,  for  Case  2  the  geometric  standard  deviations  are  much  larger  than 
those  found  for  Case  1,  E=0,  and  for  the  smallest  reaction  rate  studied 
continuously  increase  with  time.  The  solution  of  equations  ( A  —6—9) -( A-6- 
12)  for  Case  1,  E=0,  give  variations  of  geometric  mean  primary  particle 


diameters  and  standard  deviations  with  time  and  reaction  rate  shown  in 


Figures  A-6-3  and  A-6-9.  These  results  agree  qualitatively  with  tr.e 
experimental  observations  that  mean  particle  diameters  decrease  with 
increasing  reaction  rates  and  the  particles  are  narrowly  distributed  in 
si  ze . 

These  numerical  and  experimental  results  are  consistent  with 
the  hypothesis  that  primary  particle  growth  takes  place  primarily  by  the 
process  of  condensation  of  atoms  onto  nuclei  and  that  collisions  between 
clusters  is  suppressed.  This  suppression  is  attributed  to  the  electrical 
double  layer  forces  between  clusters. 

c .  Discussion  . 

This  study  has  shown  that  primary  ferromagnetic  particles  with 
diameters  from  5  nm  to  1  ,m  and  with  saturation  magnetizations  from  0  to 
1 30  emu/g  can  be  prepared  by  control  of  various  experimental  parameters 
including  the  metallic  cation(s),  the  metallic  cation  concentration,  the 
reaction  temperature,  and  the  use  of  additives.  The  tendency  of  the 
primary  particles  to  form  chains  is  strongly  affected  by  the  saturation 
magnetization,  the  size  of  the  primary  particles,  and  the  presence  of 
external  magnetic  fields.  The  chaining  mechanism  was  studied  indirectly 
through  determination  of  the  fractal  dimensions  of  iron  and  cobalt  aggre¬ 
gates  of  known  saturation  magnetization  formed  in  the  absence  of  an 
external  magnetic  field;  iron  aggregates  with  higher  saturation  mag  - 
netization  ( 1  j0  emi.’g;  had  a  lower  fractal  dimension  .'■</  than  cobalt 
-i  g  g  regales  ;  f  r  acta'  1  1  m  -  r.  s  :  ■.  t.  <.  t  1 .  'a  w  1 *  n  !  ■  r  ■  ■  i  i  r  ■ t  ’  n 
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Figure  A-6-8. 


Variation  of  geometric  mean  primary  particle  diameter 
with  time  for  three  values  of  chemical  rate  coefficient 
for  Case  2,  collision  efficiency  between  clusters,  E=0. 
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Figure  A-6-9. 
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Variation  of  geometric  standard  deviation  with  time  for 
three  values  of  chemical  rate  coefficient  for  Case  2, 
collision  efficiency  between  clusters,  E=0. 
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magnetization  (30  emu/g).  Volume  and  surface  compositions  of  primary 
particles  were  determined  by  Inductively  Coupled  Plasma  Emission  Spectro¬ 
scopy  and  X-ray  Photoelectron  Spectroscopy.  It  was  found  that  atomic 
ratios  of  metal  to  boron  was  around  2  and  that  the  boron  was  concentrated 
on  the  surfaces  of  the  particles.  Alloy  particles  prepared  from  equi¬ 
molar  Fe-Co  and  Fe-Ni  solutions  always  were  enriched  in  Co  or  Ni,  indi¬ 
cating  a  higher  reduction  rate  for  Co'M‘  and  Ni‘M‘  as  compared  to  Fe++. 

Two  mechanisms  for  formation  for  primary  particles  have  been 
identified.  One  is  found  at  low  ferrous  cation  concentrations  with  basic 
borohydride  solutions;  in  this  regime,  metallic  particle  formation  ap¬ 


pears  to  occur  in  sites  in  net-like  hydroxide  gel  structures.  The  par¬ 
ticles  formed  in  this  regime  do  not  form  chains  in  the  absence  of  an 
external  magnetic  field  and  are  of  multidomain  size  with  a  crystalline 
diffraction  pattern.  At  higher  metallic  cation  concentrations  in  basic 
borohydride  solutions,  formation  of  ferromagnetic  primary  particles 


occurs  directly  from  the  metal  atoms  produced  by  the  reduction  reaction; 
these  particles  are  of  single  domain  size  and  do  not  exhibit  a 
crystalline  diffraction  pattern.  In  the  absence  of  sodium  hydroxide  in 
the  reducing  solution,  the  mean  primary  particle  size  is  not  affected  by 
variation  of  cation  concentration;  in  this  case  tin-  borohydride  hydroly¬ 
sis  reaction  is  important. 

Mean  primary  particle  size  decreased  and  saturation  magnetiza¬ 
tion  decreased  as  the  cation  was  changed  from  Fe  to  Co  or  Ni ;  alloy 


particles  of  Fe+Co  and  Fv+N'i  similarly  decreased  in  mean  size  and 
saturation  magnetization  as  the  fraction  of  iron  in  the  reactant  solution 
was  decreased.  The  addition  of  chromium  chloride  to  the  ferrous  chloride 
reactant  solution  resulted  in  primary  particles  that  exhibited  no  oxida¬ 
tion  upon  exposure  to  air  for  six  weeks;  in  the  absence  of  chromium 
chloride,  collections  of  particles  were  pyrophoric  and  iron  oxide  was 
identified  in  micrographs  of  individual  particles  after  exposure  to  air 
for  one  day.  In  all  cases,  increase  of  the  reaction  temperature  resulted 
in  smaller  primary  particles.  When  the  reduction  reaction  was  carried  out 
in  solvent  mixtures  of  water  and  ethyl  alcohol,  mean  primary  particle 
size  decreased  as  the  fraction  of  ethyl  alcohol  increased. 

A  relationship  has  been  demonstrated  between  fractal  dimension 
of  aggregates,  their  saturation  magnetization,  and  the  chaining  tendency 
of  the  primary  particles.  Aggregates  of  primary  iron  particles  with 
saturation  magnetization  of  1  30  emu/g  had  a  fractal  dimension  of  1.5M 
whereas  aggregates  of  primary  cobalt  particles  with  saturation  magnetiza¬ 
tion  of  30  had  a  fractal  dimension  of  1.72.  Linear  chain  particles  have 
a  fractal  dimension  of  1.0. 

A  model  for  growth  of  primary  particles  above  the  critical 
cation  concentration  in  basic  solution  has  been  developed  which  gives 
results  in  agreement  with  ex perimenta 1  observations.  The  model  indicates 
that  primary  particle  „*  r  w t h  takes  place  in  this  regime  primarily  by 
condensation  of  atoms  or, to  nuclei  and  that  collisions  between  clusters  is 


probably  suppressed  L  y  ••  1  •  •  f  r  i o a  1  double  layer  forces. 
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process  whose  simulation  is  of  practical  importance.  ••  ,.ecus  fogs  from 
cooling  towers  ini  drier  industrial  sources  can  1  reduction  in 
visibility,  with  concomitant  ha  carls.  In  lnflustna. ,  1  -ultural,  and 

military  technology  an  understanding  of  the  dynamics  of  volatile  aerosol 
plumes  is  necessary.  Realist  ic  simulation  of  the  pro  •>  saes  involved  ;n 
these  applications  is  possible  only  by  numerical  te-'hn ;  ques .  Conse¬ 
quently,  this  paper  reports  our  investigation  of  rei  i«  .-  and  accurate 
numerical  schemes  for  such  simulation. 

.tie  problems  arising  in  numerical  simulation  •;  tne  d  ispersi  r. 
of  a  plume  of  volatile  aerosol  are  clearly  shown  by  st.::  ..y,  a  crosswind 
line  source.  Tins  choice  obviates  the  need  for  a  sutg :  :  e  mdel  to 

approximate  3  point  source;  any  error  is  strictly  i  •  n  .  n.  *■  r  :  a  . 

methods  used  in  the  simulation.  Also,  a  number  of  ;  •  •  .  1 .  s  s  •  : 

aerosol  plume  dynamics  are  a , a  i  1  a  t 1 e  for  t  »  •  .  ,  ’  .  .  n- 
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First,  the  relevant 


I'uat’.ons  for  this  process  are  discusser]. 


Next,  trie 


numerical  schemes  f^r  tneir  solution  are  outlined.  Finally,  the  numeri¬ 
cal  schemes  are  tester:  against  analytical  solutions  and  are  usel  m  the 
simulation  iispersion  under  realistic  atmospheric  conditions. 

a .  Theory  . 

For  clarity,  we  choose  to  describe  the  dispersion  of  an  aerosol 
consisting  of  a  single  chemical  species  in  inert  host  gas.  The  aerosol 
concentration  is  sufficiently  low  that  it  may  be  regarded  as  a  passive 
contaminant.  Also,  the  evaporation  process  is  sufficiently  slow  that 
that  evolution  of  the  plum*3  may  be  treated  as  an  isothermal  process.  In 
addition,  it  Is  assumed  that  the  turbulent  processes  can  be  represented 
by  first-order  closure  assumptions. 

Witn  these  restrictions  we  study  the  evolution  of  the  particle 
number  uensity  function  n(m,x,z,t),  where  n( m ,x,z ,t)  dm  is  the  number  of 
particles  having  masses  in  the  range  m,dm  at  downwind  position  x  at 
height  z  at  time  t.  The  evolution  equation  is 


n '  m,  x  ,z  ,  t  n(  m ,  x  ,z  ,  t) 


"A  z)  [  — ■ —  !■—  *  mw]  *  —  [.  (  m ,  s)  n(  m  ,  x  ,  z  ,  t )  ] 
x  m 


(  A  -  ’ . 


n(m,x  ,z,  t)  n( m ,x  ,z  ,  t) 

- J  ♦  I - 


1 AV, 


v-y-u.,. 
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This  equation  is  coupled  to  the  conservation  equation  for  the  vapor 


which  we  give  in  terms  of  the  saturation  ratio  s: 

•s(x,z,t)  s(x,z,t)  -s(x,z,t) 

-  ♦  U  ( z)  [ - ]  =  —  K  ( z)  [ - ] 

•  t  'X  ■  z  ■  z 

1 

-  —  .  ( m,s)  n(  m,x  ,z  ,  t)  dm .  (a-7-2) 

cv 

Equations  (  A  —7  —  1 )  and  (A-7-2)  are  subject  to  the  following  initial  and 
boundary  conditions: 


n(  m,x ,z ,0) 

=  0, 

s(x,z,0)  = 

0, 

(  A-7-3) 

n(m,  0,h  ,t) 

=  [Q0/U(h)l  (z-h)  f (  m)  , 

( A-7-4) 

s(0, h,t)  = 

1.0, 

(A-7-5) 

n(m,x,  H,t) 

•s(x, H , t) 

(  A-7-6) 

■  % 

■  z 

nVd(m)  =  K 

■n/  ■  z,  z  =  0, 

(  A -7 -7) 

3/  Z  =  0, 

z  =  0. 

(  A  —7  —8 ) 

U(z)  is  the  x  component  of  the  mean  fluid  velocity,  and  K(z)  is  the 
vertical  eddy  diffusivity.  G-.C  m)  is  the  gravitational  settling  speed  for 
a  particle  of  mass  m.  Qq  is  the  source  strength,  h  its  height,  and  f(  m) 
the  source  size  distribution.  H  is  the  mixing  height,  Vd(tn)  is  the 
deposition  velocity  of  a  particle  with  mass  m. 

.(m,s)  is  the  condensation/  evaporation  growth  law  for  a  par¬ 
ticle.  We  use  the  following  approximate  expression,  due  to  Fuchs  and 
qq 

-  utug m : 


"s  „  .m1 1  -c.,(  s-eKe1 
fs  1  v 


1 . 33Kn+0.7i 


(  A-7-9: 


where  tr.e  3  if  fusion  :ve f:\  :  lent  T  ,  of  vapor  j  in  host  gas  is  related  to 
the  meat',  molecular  speed  v  1  and  the  mean  free  path  Lj  by  Dg j  =  1/3  VjLj. 
Ir.  Fruitier,  A  -  [<  ,  i~  tne  particle  density;  Kn  the  Knud  sen  number, 
Kn=2Lj/Dpj  and  L'p  is  the  particle  diameter.  The  saturation  ratio  s=c/cv» 
where  c  is  the  vapor  mass  concentration  and  cv  is  the  vapor  mass  concen¬ 
tration  corresponding  to  the  bulk  liquid  vapor  pressure  at  plume  tempera¬ 
ture  T.  Ke  is  the  Kelvin  number,  Ke:H'  v/D  kT,  •  is  the  surface  tension, 
v  the  particle’s  molecular  volume,  and  kT  the  thermal  energy. 

Some  of  the  important  observables  of  n(  m)  are  the  total  par¬ 


ticle  number  concentration  N, 


Nr/  n(  m)  dm ; 


the  particle  mass  concentration  M, 


Mr/  mn(  m)  dm; 

and  the  total  extinction  coefficient  ext>  obtained  from  the  normalized 
extinction  efficiency 


ext  '  ' 


qextrc2/^n(  m)  dm. 


rnese  are  ir.  the  fallowing  analysis  of  results. 
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routine  application  of  classical  numerical  methods.  Numerical  approxima¬ 


tion  of  these  terms  by  such  methods  is  well  known  to  lead  to  suostantial 
numerical  dispersion  and  diffusion.  Our  aim  in  this  study  is  to  obtain 
numerical  results  of  high  accuracy  by  techniques  that  enable  one  to 
reduce  to  any  desired  level  errors  such  as  numerical  dispersion  and 
diffusion. 

With  respect  to  simulation  of  the  advection  and  dispersion 
terms  in  Equation  (A-7-1),  Tsang  and  Brock  concluded  that  Fromm's 
method^' for  ajvection  and  orthogonal  collocation  on  finite  elements^ 
for  diffusion  gives  reliable  numerical  schemes.  Tsang  and  Brock'^14 
successfully  used  the  method  of  fractional  steps1^?  in  the  study  of 
coagulation  in  a  dispersing  plume  from  a  crosswind  line  source.  One  of 
the  purposes  of  this  study  is  to  demonstrate  the  utility  of  these  methods 
in  plumes  undergoing  the  condensation/ evaporation  process. 


The  condensation/ evaporation  term  has  long  posed  a  major  prob¬ 
lem  in  aerosol  dynamics  simulations.  Numerical  simulation  of  this  first- 
order  hyperbolic  evolution  equation  remains  as  a  challenging  problem  in 
computational  physics  and  engineering.  For  first-order  hyperbolic  equa¬ 
tions,  most  numerical  schemes  exhibit  either  numerical  dispersion  'mani¬ 
fested  as  oscillations  over  some  region)  or  numerical  diffusion  (dissipa¬ 
tion  of  the  amplitude  of  a  distribution  or  the  smearing  of  the  front). 
Both  numerical  dispersion  and  diffusion  are  obviously  ur.  lea :  r  a*  1  ••• 
simulating  the  evaporation/condensation  of  aerosols.  Recently,  Tsang  and 


Brock  devised  a  novel  numerical  scheme  for  sn-.h  siaului 


mzmm 


ref 


Comparison  of  the  numerical  solution  and  known  analytical  solution  shows 


that  the  scheme  gives  very  accurate  results.  Case  studies  at  various 


conditions  were  carried  out.  For  the  evaporation  problem,  it  was  found 
that  the  finite-element  method  (LFEM)  using  linear  basis  functions  and 
natural  boundary  conditions'^  gives  accurate  results.  Here  we  combine 
the  numerical  algorithm  for  advection  and  diffusion  used  by  Tsang  and 
Brock 1 04  with  the  finite-element  method  for  evaporation  to  study  a  vola¬ 


tile  aerosol  plume.  Only  steady  state  solutions  are  presented  in  this 


paper . 

The  numerical  solution  of  the  coupled  integro-par tial  differen¬ 
tial  Equations  (A -7-1)  and  (A-7-2)  presents  a  challenging  problem.  The 
simultaneous  solution  by  an  implicit  numerical  scheme  using  a  matrix 
technique  for  10^-10^  unknowns  is  no  trivial  task.  In  this  work  the 
method  of  fractional  steps  is  used  in  the  simulation.  Advection,  diffu¬ 


sion,  and  evaporation  are  included  as  follows: 

1.  Solve  the  advection  and  diffusion  equation  Jmax  times, 
where  is  the  number  of  size  classes. 

CU  d  A 


2.  Solve  the  evaporation  equation  by  LFEM  at  each  grid  point. 

3.  Solve  the  advection  and  diffusion  equation  for  s. 

4.  Calculate  the  integral  term  in  Equation  (A-7-2),  which  is  a 
source  term  for  vapor  due  to  evaporation,  and  update  the 
saturation  ratio  at  that  grid  point. 

if  the  saturation  ratio  of  the  vapor  does  not  change  with 
position  -  that  is,  is  constant  -  procedures  3  and  4  are  not  necessary. 
The  method  of  fractional  steps  just  outlined  decouples  Equations  '  A  —7  —1 ' 

!  : 


•I 


and  (A-7-2).  The  explicit  nature  of  procedure  ^  does  not  pose  a  problem 


because  the  source  term  is  counterbalanced  by  the  dilution  effect  of 
advection  and  diffusion. 

c.  Results  and  Discussion. 

In  order  to  validate  the  computational  procedure  outlined  in 
the  previous  section,  we  perform  some  convergence  studies  on  the  numeri¬ 
cal  scheme.  For  certain  forms  of  U(z),  K(z),  and  ((m),  analytical 
solutions  are  available.  Here,  we  assume 

U(z)  =  2m  sec-1,  K(z)  =  1  m2sec-1, 

and  L  =  -am1^  (this  represents  continuum  evaporation,  where 
a=2x1 0-1  1g2/^sec-1) .  We  also  assume  that  gravitational  sedimentation  and 
deposition  are  absent  for  this  test  problem;  that  is,  Gz( m) =0  and 
Vd(m)=0.  For  vertical  diffusion  all  the  numerical  results  in  this  work 
were  obtained  by  using  six  interior  collocation  points  (the  quadrature 
points  of  a  shifted  Legendre  polynomial)  in  each  finite  element  of  sizes 
2.07,  1 2- 9 3 »  and  85  m.  The  source  height  is  2  m  from  the  ground.  For 
these  restrictions,  the  analytical  solution  of  Equation  ( A  —7  — 1 )  is 
readily  obtained  by  Fourier  transform  methods  and  the  method  of 
characteristics.  Figures  A -7-1  and  A-7-2  show  the  comparison  between 
analytical  and  numerical  solutions  at  the  source  height  for  the  total 
particle  mass  and  number  concentrations,  respectively.  x  is  the  grid 
size  in  the  x  direction,  ' te  is  the  time  step  used  in  procedure  2  for 
evaporation.  Normally,  '  t„  for  evaporation  is  usually  less  than  that 
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Figure  A-7-1.  Comparison  between  analytical  solution  ( _ )  and  numeri¬ 
cal  solution  ( - • - ,  Ax=2.0  m,  tc=0.02  sec; - ,  Ax =0.5 

m,  tg=0.005  sec)  at  source  height  for  total  particle  mass 
concentration . 
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Figure  A -7  -2 .  Comparison  between  analytical  and  numercal  solutions 

( - * - ,  4  x  =2.0  m,  'te  =  0.02  sec; - ,  x=0.5  m  ,  te=0.005 

sec;  at  source  height  for  total  particle  number  concentra- 
t  ion . 


for  advection  and  diffusion.  In  this  work,  te  -  1  / 5 CX  t acj ,  a nd  so  in 
practice  we  perform  the  advection/diffusion  calculation  every  ‘-t  d  step 
rather  than  every  ' te  step.  This  saves  computation  time  without  sacri¬ 
ficing  accuracy.  Figures  A— 7  —  1  and  A-7-2  show  that  as  we  decrease  x  and 

t  ,  the  numerical  solution  converges  to  the  analytical  solution  conveives 
e 

to  the  analytical  solution,  thus  demonstrating  the  validity  of  the  method 
of  fractional  steps  in  volatile  plume  calculations.  Figure  A-7-3  shows  a 
comparison  between  analytical  and  numerical  solution  at  the  source  height 
for  the  number  density  function,  dN/dlogDp,  for  different  downwind 
distances.  The  excellent  agreement  with  the  analytical  result  shows  that 
the  numerical  method  is  accurate.  One  of  the  difficulties  in  solving  the 
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Figure  A-7-3-  Comparison  between  analytical  solution  (curves;  a 
numerical  sol  u  t  ion  (  sy  mb  ol  s)  at  source  height  for  number 
density  function,  dN/dlogDp,  for  different  down  vi  n  J 
distances  (-0-,  50  m;  —  100  m;  -"  •  2C0  rr.)  . 


A 


initial  boundary  value  problem  of  first-order  hyperbolic  equations  lies 


in  the  choice  of  appropriate  boundary  conditions  and  this  is  crucial  in 
obtaining  meaningful  and  accurate  solutions.  In  Figure  A-7-3,  even  at  a 
downwind  distance  of  200  m,  it  is  obvious  that  there  is  no  boundary 
effect  between  particle  diameters  of  10“6  and  10  "5  cm.  The  minimum  D 

r 

used  in  the  simulation  is  10“^  cm.  Since  computations  must  be  performed 
on  finite  domains  of  particle  diameters,  and  evaporation  increases  the 
number  of  smaller  aerosol  particles,  this  result  shows  that  natural 
boundary  conditions  are  appropriate  in  simulating  evaporation  of  an 
aerosol . 

With  the  tests  of  the  numerical  scheme  completed  for  advection, 
diffusion,  and  evaporation,  we  use  realistic  atmospheric  wind  and  dif- 
fusivity  profiles  for  approximating  evaporating  plumes  of  simulated  light 
and  heavy  oils  (diesel  and  fog  oil).  The  vertical  eddy  diffusivity 
profile  is  taken  from  the  data  presented  by  Smith. The  wind  speed  is 
from  field  measurements  at  Cape  Canaveral  under  neutral  stability 
(Panofksy ) .  Deposition  velocity  is  taken  from  the  data  of  Sehmel.1®2 
We  assume  the  mixing  height  to  be  100  m.  The  source  height  is  2  m  from 
the  ground.  Thie  line  source  size  distribution  is  characterized  by  a  log 
normal  distribution  with  source  strength  of  0.3307  g  m-1  sec-1,  geometric 
mean  diameter  of  0.6  m,  and  geometric  standard  deviation  of  2.24.  The 
particle  density  is  assumed  to  be  1.0  g  cm-3.  For  the  heavy  oil  <  fog 
oil,  or  100  pale  oil)  the  average  molecular  weight  is  300,  the  surface 
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tension  is  35  dyn  cm-'  and  a  constant  bulk  vapor  pressure  of  9.6x10“^  mm 
Hg  is  assumed.  For  calculation  of  ext>  the  particle  refractive  indexes 
vr  with  associated  radiation  wavelengths  are  vr  =  1. 5077  -2.  94  xl  0  , 

A=3. 39  ,..m  and  v  p=  i .  44  -0. 1i  ,  '•  =0. 44  1  6  m.  For  the  light  oil  (diesel)  the 
average  molecular  weight  is  221,  the  surface  tension  is  20  dyn  cm-1,  and 
a  constant  bulk  vapor  pressure  of  1.5x10-3  mm  Hg  is  assumed.  The 
particle  refractive  indexes  with  the  associated  radiation  wavelengths  are 
v  p=  1 . 522  -0. 064  i  ,  ''=3-5  urn  and  v  r=  1 .51  -O.Oi ,  '-=0.5cm. 

In  order  to  demonstrate  the  roles  of  vapor  conservation  and  the 
Kelvin  effect,  we  carried  out  simulation  studies  of  atmospheric  plume 
dispersion  at  saturated  conditions,  s  =  1.0,  and  with  s=s(x,z)  according  to 
Equations  (  A  —7  —  1 )  and  (A-7-2).  Note  that  when  s=1.0,  only  evaporation 
occurs  owing  to  the  Kelvin  effect.  Figures  A -7 -4  and  A -7 -5  show 
isopleths  of  total  single-particle  extinction  (a  sensitive  index  of 
changes  in  the  density  function)  for  the  plume  cross  section  ats  =  1.0  for 
the  two  wavelengths.  Figures  A-7-6  and  A-7-7  are  companion  plots  for  the 
complete  realistic  simulation  with  vapor  conserved  under  otherwise 
identical  conditions.  For  the  infrared  wavelengths  with  the  initial 
particle  size  distribution  employed,  cext  very  sensitive  to  the  par¬ 
ticle  concentrations  at  the  larger  diameters;  these  particles  disappear 
much  more  rapidly  for  the  conserved  case  than  for  the  case  s  =  1.0,  with 
attendant  rapid  decrease  of  cext  for  the  infrared  ijn  the  conserved  case. 

Figure  A -7 -8  shows  for  the  conserved  case  the  number  density 
function,  dN/d  leg  Dp,  for  the  oil  aerosols  at  a  downwind  distance  of  100 
m  at  the  source  height  of  h=2  m.  Under  the  conditions  studied,  the 
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Figure  A-7-4.  Isopleths  of  total  single- particle  extinction  cext  in  the 
plume  from  a  crosswind  line  source  for  s=1.0,  ext  ^ 
cm-1)  is  equal  to  values  given  on  isopleths  times  5x10" 

( _ ,  fog  oil,  >=0.441  um; - .diesel,  X=0.50Um). 
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Figure  A -7-0. 


1.-th s  of  total  single  particle  extinction  ext  m  t 
e  from  a  cross  wind  line  source  for  s=1.0,  ~ext  (in 
1  is  equal  to  values  given  on  isopleths  times  5x10 
'or  oil,  =7. <9  m; - .diesel,  '  =3-50  m)  . 
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Figure  A-7-6.  Isopleths  of  total  single  particle  extinction  oext  in  the 
plume  from  a  crosswind  line  source  with  vapor  conserved, 
cex t  (  V  cm  '  isequal  tovaluesgivenon  isopleths  times 
5x10~5  ( - -  fog  oil,  ’'-=0.441  um; - ,  diesel,  *  =0.50 

i-m)  . 


Figure  A-7-7.  Isopleths  of  total  single  particle  extinction  '  ex+  in  the 
plume  from  a  crosswind  line  source  with  vapor  conserved, 
ext  (in  era”1)  is  equal  to  values  given  on  isopleths  t  i  t,  t  s 
5x10"  ( - -  fog  oil,  •  =3*39  ..os:  - .diesel.  •  ='•. 


Kelvin  effect  causes  a  rapid  increase  of  evaporation  rate  with  decreasing 
particle  diameter.  As  indicated  in  Figure  A-7-8  for  the  heavy  oil  (fog 


oil},  this  causes  an  overall  decrease  in  N  and  an  increase  of  dN/d  log  Dp 
at  smaller  diameters,  whereas  dN/d  log  Dp  is  hardly  altered  by  evapora¬ 
tion  at  the  larger  diameters.  The  vapor  pressure  of  the  light  oil  is 
sufficiently  large  that,  as  shown  in  Figure  A-7-8,  at  x  =  100  m  the  entire 
density  function  is  shifted  to  smaller  diameters. 


Figure  A-7-8.  Number  density  function  dN/d  log  Dp  for  light  and  heavy  oil 
aerosols  in  plume  from  a  crosswincf  line  source  at  x  =  100  m, 

h-2  m  ( _ ,  source  distribution; - ,  no  evaporation  at 

100  m;  --0--,  evaporation  of  fog  oil  at  100  ra,  -  ' 
evaporation  of  diesel  at  100  m)  . 


Conclusions . 


It  has  been  demonstrate":  that  .1  isjvrsior,  ••  :  .  ;  a r  :  :  e 

aerosol  can  be  simulated  numerically  by  the  method::  -•■.  r 

the  special  case  where  an  analytical  solution  is  ;■  -  —  .  •  •  ;  ►  e«n 

shown  that  the  numerical  solution  converges  to  the  a-.  -  solution 

with  a  reduction  in  step  size;  consequently  the  . jt.  erica  1 

simulation  can  be  r  ed  uc  ed  to  any  desired  1  ev  e  1 .  Ir.  t  r  •  .t.  ;  »r i sons  , 

the  validity  of  the  method  of  fractional  steps  has  beer.  :  -  ■  ,■  *•-.  »•;  tr.  the 
attendant  advantage  over  fully  implicit  schemes.  From  ••  it. her 

studies  it  has  been  demonstrated  that  the  finite-ele-:  •  '  m  •  t  r.o  i  using 

linear  basis  functions  and  natural  boundary  conditions  s  ice  urate 

numerical  solutions  for  the  evaporation  process. 

With  the  methods  described  in  part  here,  ••  foresting 

problems  in  aerosol  dynamics  can  now  be  studied.  Ft- include 
simulations  of  dispersion  of  plumes  in  which  both  coagul  *  .  ar.  i  various 
cond  ensa  tion/ evaporation  processes  occur.  It  would  a’.  •••  ‘  interest 

to  extend  these  procedures  to  other  problems  such  as  t.v  ri,  ;i  condensa¬ 
tion/evaporation  processes  in  nonisothermal  plumes,  ■  ‘  ’  ~t.  of  some 
95 


importance . 


'ner-'  ir  ;rrcr,tl  y  interest  in  the  nonlinear  ir.  ter  an  *  i 
i  :  .  t .  v  •  :  it.  .•  •.  r.  3  : ;  er:  3  ed  particles.  For  very  intense  bem, 

p  •»  r  t  :  1  •  j  •  Dren>.  :  «r.  an  3  propagation  of  a  pi  asm  a  1  u  x  r: . 

s  ox  e  w  r.  a  t  r  :  >-  -»x  m  T  *  ns :  t  les  ,  a  particle  may  explode  or  fraotur- 

f)ri  ir.g  at  ring  m  t  ne  viscosity  of  the  parti  •  1  e.  At  1  ^  w 

Intensities,  particles  will  evaporate  rapidly  with  only  stn  a  1  . 
J  ist  urbar.ces  :r.  the  host  gas.  Work  reported  here  is  concerned  with  the 
1  a st  *  .  •  •  "  he  interact  ion  regimes  noted  above.  First,  work  or.  particle 

exp  1  s.  m  w t  .  ".  f>e  j escribe 3.  Then,  work  on  fast  evaporation  of  jroplet.s. 

- '  .  -xpl~si  r.  of  1  water  droplet  by  pulsed  laser  heating  . 

a  .  >  *  rod  u:  ti  or. 


T.~  simulation  discussed  here  addresses  the  following  problem: 
A  st-.;.  iroplet  uniformly  absorbs  energy  from  a  laser  pulse.  The 

:  ur  i*  :  r.  *  -  :>•  ;  . .  se  is  ■>  snort  that  no  fluid  motion  occurs  during  the 

pus-.  "he  jr,;  .-:  c,rc  -  .-rough  energy  that  its  final  state  after  the 

p-l  u- ,  •  ,t  before  my  x  t:cr.,  is  w«  1 1  above  its  critical  point.  Tne 

iroplet  s,.:  ae  ;  uen  tl  y  expires.  We  solve  the  conservation  equations  for 

the  r , ;  *  ; v  j.pn-r  it  t :  .  y  symmetric  explosion.  Uniform  ,  inst  antane  •  is 
h1’  a  f  i .  .  -  .  -  t  ■  >  .  n  •  : ;  ;  e  tne  hea  ting  process.  "h  u«  ,  we 

invert  id  i  *  -  V  v  s ;  •  r-  •  ••  dd-t  with  the  fluid  mot  >  »>?:  after 


»  . 


Our  resuits  show  that  in  the  explosive  regime,  large  nucleat:or 
rates  can  develop  behind  the  shock  of  the  exploding  droplet.  This  pr.j.:- 
transition  can  strongly  influence  b  >th  the  hydrodynamics  and  the  electro¬ 
dynamics,  and  thus  needs  to  be  accounted. 

In  addition,  a  number  of  interesting  experimental  phtr. 
EMR-droplet  interactions  have  been  reported  in  the  literature, 
convective  motion  is  important  to  all  of  these  phenomena.  A  realise, 
simulation  of  explosive  vaporization  can  help  to  explain  tt.e  -o 
occurrences . 

We  beg  in  with  a  brief  review  of  previous  work  on  EM  F  -d  r«.  1  :  ' 

interactions.  Next,  we  discuss  the  theory  of  these  interactions  an 
apply  the  theory  to  the  explosive  vaporization  of  single  water  droplets 
Finally,  we  present  the  results  of  our  numerical  simulation. 

Many  researchers  are  currently  investigating  how  a  very  small 
droplet  responds  to  applied  electromagnetic  radiation,  and  how  t:.< 
droplet  subsequently  alters  the  electromagnetic  fields.  This  problem 
falls  within  the  context  of  the  larger  problem  of  how  to  model 
predict  the  effect  of  a  system  of  discrete  attenuators  on  the  propagati*' ! 
of  EMR.  Thoroughly  understanding  this  phenomenon  will  lead  to  adva--.:v; 
in  the  microtechnologies  (for  example,  the  manipulation  and  eo..T:  ; 
submicrometer  sized  droplets,  as  in  optical  levitation’1  as 
advance  our  knowledge  in  the  areas  of  '  nrer  propagation  ♦  h.rcu.. . 
or  cloudy  atmosphere,  in  situ  analysis  of  combusion  pr:dn. 
astrophysics . 
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Broadly  speakirg,  EM R  interacting  with  a  small  droplet  and  its 
environment  (in  this  case  the  atmosphere)  manifests  two  types  of 
response:  electrodynamic  and  hydrodynamic.  The  Maxwell  equations  deter¬ 

mine  the  electrodynamic  response;  the  conservation  equations  for  mass, 
momentum,  and  energy  determine  the  hydrodynamic  response.  These  partial 
differential  equations,  along  with  the  constitutive  relations  that  govern 
the  physical  properties  of  the  system,  and  an  equation  of  state  determine 
the  behavior  of  the  system.  The  two  sets  of  equations  are  strongly 
coupled  through  source  terms  in  the  primitive  equations  and  through  the 
constitutive  relations. 

We  are  currently  working  on  the  hydrodynamic  aspects  of  this 
problem.  To  a  first  approximation,  the  effect  of  EMF  on  a  small  droplet 
is  to  generate  heat  within  the  droplet  due  to  resistive  heating.  The 
rate  at  which  the  heat  is  generated  determines  the  response  of  the 
droplet.  The  heating  rate  can  vary  over  a  wide  range  as  determined  by 
the  complex  index  of  refraction,  the  intensity  of  the  radiation  source, 
and  the  optical  size  of  the  droplet  (Pendleton12^  anCj  Greene  et  al.11^). 
The  optical  size,  also  known  as  the  Mie  size  parameter,  is  defined  as  x  = 
d/  ■■  where  d  is  the  droplet  diameter  and  is  the  wavelength  of  the 
incident  light.  The  droplet  shows  an  equally  wide  range  of  response, 
from  slow  evaporation  to  explosion.  Thus,  solving  the  problem  for 
artibrary  radiation  intensity  is  very  complicated. 

A  large  body  of  literature  has  developed  considering  many 
aspects  of  this  problem.  For  example,  the  articles  by  Yudanin  and 
Lax  ,  1  ^  ?r :  sh  i  va  1  ko  ,  1  -  and  Ch  i  ta  n  v  i  s  1  1  7  illustrate  the  wide  range  of 


effects  that  are  being  considered.  Most  of  the  current  wer*  addressing  a 


realistic  simulation  of  the  hydrodynamics  has  dealt  with  the  slow 
vaporization  regime  (see  for  instance  A  rm  s  t  rong 1  ^  9,  1  '  anc  Tag  eev  and 

•»  i  ]  ' 

Sein  te  Id  ‘  “  ) . 

We  are  exploring  the  explosive  vaporization  regime  of  this 
problem.  This  regime  is  interesting  because  experimental  studies  have 
uncovered  many  unusual  effects  that  seem  be  be  based  at  least  partly  on 
strong  convection.  Among  these  are  the  lowering  of  the  breakdown 
threshold  for  air  containing  small  particles,  and  the  appearance  of  jets 
of  material  emanating  from  a  droplet  (Chang1 1^)).  Under  some  conditions, 
this  behavior  can  be  explained  completely  from  hydrodynamic 
considerations,  regarding  the  EMR  simply  as  a  source  of  heat.  For 
instance,  the  lowering  of  the  breakdown  potential  has  been  attributed  to 
shock  heating  of  the  atmosphere  (Steverding12?) .  Also,  the  appearance  of 
"cloudy"  jets  can  probably  be  explained  in  terms  of  rapid,  very  localized 
vaporization  at  a  hot  spot  near  the  surface  of  the  droplet. 

We  have  chosen  a  simple  but  very  dynamic  test  problem  as  a 
first  step  towards  developing  a  general  solution  for  the  convection 
dominated  regime.  By  choosing  a  relatively  simple  hydrodynam  ical  system, 
we  can  concentrate  on  simulating  the  physical  properties  of  the  droplet 
realistically.  This  in  turn  will  help  us  determine  what  other  physics 
(such  as  nucleation  and  particle  growth)  might  arise  due  to  strong 
convective  effects.  The  physical  system  examined  consists  of  a  pure 
water  droplet  in  air  as  a  host  gas.  Initially,  the  droplet  is  subjected 
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to  a  very  short ,  very  intense  pulse  of  loser  radiation.  The  optical  size 
of  the  droplet  is  assumed  to  be  in  the  Rayleight  limit,  thus  it  absorbs 
energy  uniformly.  Because  of  the  very  short  duration  of  the  pulse,  the 
droplet  remains  locked  in  place  during  the  entire  heating  period.  This 
implies  that  the  droplet  heats  at  constant  volume.  Additionally,  the 
intensity  of  the  pulse  is  great  enough  that  the  energy  absorbed  drives 
the  water  droplet  into  a  supercritical  fluid  state.  Because  of  uniform 
absorption,  the  motion  that  evolves  will  be  spherically  symmetric,  and  a 
one  dimensional  simulation  is  sufficient  to  describe  the  flow. 

Under  these  conditions  the  hydrodynamic  simulation  can  be 
started  after  the  pulse  but  before  the  onset  of  fluid  motion,  using  the 
very  hot,  dense  fluid  as  the  initial  conditions.  In  practice  tnis 
situation  can  be  realized  with  an  aerosol  droplet  whose  Mie  size 
parameter  is  less  than  about  0.  3,  and  with  a  pulse  length  whirl;  is  much 
shorter  than  the  time  required  for  an  acoustic  wave  to  propagate  across 
the  droplet.  The  conditions  required  for  this  type  of  "  vaporiza  tion" 
have  been  discussed  in  the  literature  (Armstrong110  or  Zuev  ana 
Zemlyanov  . 

Heating  the  droplet  at  constant  volume  causes  its  pressure  to 
increase  enormously  (as  high  as  30000  atmospheres  for  some  cases  that  we 
nave  run,'.  The  high  pressure  causes  the  droplet  to  expand  violently 
against  its  host  gas.  The  strong  dynamics  give  rise  to  shock  waves  ••no 
associated  chock  heating  in  some  regions,  and  adiabatic  expansion, 
cooling  ,  a  r. :  tie  possibility  o  f  r.  uclea  tion  in  other  regions.  '**•'  will 
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consider  these  possibilities  after  we  discuss  the  equations  which  govern 
the  solution. 

b .  Theory  . 

The  conservation  equations  for  mass,  momentum,  and  energy 
determine  hydrodynamic  flows.  Written  in  conservation  form,  these  are 
-'.Bird,  Stewart,  and  Light  foot 1  1 1 ) : 

—  +  (  V  ,  ±V)  =  -  ■  '  ji  i  =  1 ,  n  (B-1  -1 ) 

t 

—  V  +  (  V  *  -  VV)  =  -  VP  -  •  •  (B-1 -2) 

•t 

—  cE  +  (  V*  cEV)  =  -  V  -  ?v  -  7*  qp  (B-1 -3) 

•t 

-  V  •  qc  -  (  T  *  V) 

where : 

.•  is  the  mass  density  of  the  i’ th  component 

n  is  the  number  of  components 

is  the  total  mass  density 

V  is  the  mass  average  velocity 

is  the  diffusive  mass  flux  of  component  i.  with 
respect  to  the  mass  average  velocity 

is  the  shear  stress  tensor 

-  is  the  fluid  pressure 

is  the  total  energy  per  unit  of  mass,  kinetic  plus  internal 


* ^  • 


is  the  radiative  heat  flux 


3r 

qc  is  the  conductive  heat  flux 

These  equations  represent  a  pointwise  mass  balance  for  each  component,  a 
force  balance,  and  a  total  energy  balance  respectively.  We  have  neglected 
terms  containing  gravity  and  are  not  allowing  the  possibility  of  chemical 
reaction.  Both  of  these  can  be  included  when  necessary. 

Notice  that  the  momentum  equation  (B-1-2)  is  nonlinear  in 
velocity.  This  gives  rise  to  much  interesting  behavior.  In  fact,  under 
many  flow  conditions  this  nonlinear  term  dominates  the  behavior  of  the 
solution.  We  will  examine  this  further  shortly. 

In  the  explosive  vaporization  regime,  molecular  transport 
effects  such  as  mass  diffusion  and  heat  conduction  are  small  compared  to 
convective  transport.  Convective  transport  is  driven  by  gradients  in  the 
pressure  field.  In  the  d ropl et/ atm ospher e  system,  convection  becomes 
important  when  the  vapor  pressure  of  the  evaporating  droplet  becomes 
comparable  to  the  pressure  of  the  surrounding  gas.  Under  these 
conditions,  the  rate  of  evaporation  becomes  greater  than  the  rate  at 
which  diffusion  can  "conduct"  mass  away  from  the  droplet.  As  the 
temperature  and  therefore  the  vapor  pressure  and  evaporation  rate  of  the 
droplet  rises,  convection  becomes  an  effective  mechanism  for  mass 
transport  away  from  the  droplet.  At  higher  temperatures  it  becomes  the 
dominant  mechanism. 

When  a  system  is  convection  dominated,  the  conservation 
equations  can  be  simplified  by  eliminating  the  terms  associated  with 
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molecular  transport.  Specifically,  the  last  term  on  the  right  hand  side 


f 
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of  equations  (  B  —  1  —  1 )  and  (B-1-2),  and  the  last  two  terms  on  the  RHS  of 
equation  (B-1-3)  are  eliminated.  Since  there  is  no  molecular  diffusion, 
the  system  of  equations  for  mass  conservation  (1)  can  be  added  to  give  a 
single  equation  giving  the  total  mass  density  at  a  point.  The  resulting 
equations  are  called  the  "inviscid  Euler  equations."  The  driving  force 
for  these  equations  is  the  pressure  gradient  term  in  equation  (B-1-2). 

Current  work  in  the  slow  vaporization  regime  incorporates 
molecular  transport  effects,  but  assumes  a  constant  pressure  field,  thus 
eliminating  convection.  The  stipulation  of  slow  vaporization  validates 
this  assumption.  The  basis  for  the  two  extremes  of  vaporization  in  a 
single  aerosol  droplet  thus  becomes  evident.  One  extreme  occurs  under 
conditions  of  constant  pressure  and  is  dominated  by  molecular  transport. 
This  arises  under  low  radiation  intensities.  The  other  extreme  occurs 
when  high  intensities  cause  large  pressure  gradients,  and  is  dominated  by 
convective  transport.  There  is  a  large  grey  area  between  these  two 
extremes  in  which  both  transport  mechanisms  must  be  incorporated. 

As  mentioned  earlier,  the  gradient  of  the  pressure  drives  the 
fluid  motion  through  the  momentum  equation.  The  extent  to  which  the 
pressure  is  reliably  predicted  for  a  compound  determines  how  well  the 
hydrodynamic  simulation  corresponds  to  observed  behavior.  With  this  in 
mind,  we  have  incorporated  an  equation  of  state  which  reliably  predicts 
the  equilibrium  properties  of  liquid,  vapor,  and  supercritical  fluid 
water  over  a  very  wide  range  of  temperatures  and  densities  (Kestin  et 
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'-•jar  e t  a  1 1 2^)  .  This  equation  of  state  allows  us  to  account 


for  deviations  from  ideal  gas  behavior. 

The  equation  of  state  is  given  as  Helmholtz  free  energy  (A)  as 
a  function  of  temperature  and  density.  The  pressure  can  be  derived  from 
this  expression  as  a  partial  derivative: 

A  =  A(T,r)  (B-1-4) 

a  A 


P  =  c2  (  — ) 

3p  , 


( B-1  -5) 


The  equation  of  state  predicts  the  proper  fluid  pressure  even  under 
conditions  of  superheating  or  supercooling.  Thus,  we  can  determine  if 
the  vapor  becomes  supersaturated,  and  its  degree  of  supersaturation. 
This  helps  us  determine  whether  condensation  is  an  important  mechanism. 
To  do  this  we  have  computed  nucleation  rates  at  every  point  in  the  system 
as  a  function  of  time.  Classical  nucleation  theory  yields  the  following 
expression  for  nucleation  rate  (McGraw  and  Saunders  ,22)  : 


JCf,)  =  K  (  f , )  exp  [-W(g»)/kT]  (B-1 -6) 

The  prefactor  K  incorporates  information  about  the  kinetics  of  the 
collision  process  giving  rise  to  nucleation: 


K(  f  i )  =  (af^/Cq)  [2c/~m^W2  (B-1 -7) 

The  activation  energy  for  the  formation  of  a  critical  sized  cluster  is 
given  from  thermodynamic  considerations  by: 
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3  ( .  ,kT  lnS)‘ 


f ]  is  the  concentration  of  the  vapor 


g*  is  the  critical  cluster  size,  the  minimum  stable  size  for  a 
liquid  particle  at  (T,P,  ) 

k  is  the  Boltzmann  constant 

a  is  an  accommodation  coefficient  or  sticking  factor 
is  the  molecular  number  density  in  the  cluster 

0  is  the  surface  tension 

m1  is  the  molecular  mass 

S  is  the  super sa turation  ratio,  the  ratio  of  the  vapor 
concentration  to  that  which  exists  at  equilibrium  over  a 
flat  interface  at  T 

The  inviscid  Euler  equations  have  been  the  bane  of  numerical 
analysts  for  decades.  The  equations  are  coupled,  hyperbolic  and 
nonlinear.  The  hyperbolic  nature  of  the  equations  manifests  wavelike 
behavior  (or  fronts)  in  the  solution.  As  an  equation  becomes  more 
hyperbolic  in  nature,  its  numerical  solution  becomes  more  susceptible  to 
spurious  oscillations  (Carey11-*)  especially  in  regions  of  the  solution 
that  contain  high  curvature  (Roache12^).  The  nonlinear  term  in  Equation 
(B-1-2)  exacerbates  the  difficulties  by  generating  shock  waves  which  are 
discontinuities  in  the  solution.  A  discontinuity  has  essentially  an 
infinite  curvature,  giving  rise  to  severe  oscillations  and  ultimately 
causing  most  standard  numerical  methods  to  fail. 

Much  research  has  gone  into  developing  methods  that  can  provide 
stable  and  well  resolved  solutions  in  the  presence  of  shock  waves. 
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Woodward  and  Collela^28  provide  a  very  good  review  and  comparison  of 
several  state  of  the  art  numerical  schemes  for  solving  the  Euler 
equations.  We  have  chosen  the  Flux  Corrected  Transport  ( FCT)  method 
devised  originally  by  Boris  and  Book. ^3  Very  simply,  FCT  works  as 
follows:  Two  solutions  to  the  problem  are  computed.  One  is  a  low  order 

solution  which  is  guaranteed  not  to  blow  up  under  any  conditions,  but 
which  gives  low  accuracy.  The  second  solution  is  of  higher  order  and 
gives  accurate  results,  but  tends  to  develop  wiggles  in  highly  structured 
(shocked)  regions  of  flow.  The  algorithm  searches  for  regions  in  which 
the  high  order  solution  develops  extrema  which  do  not  occur  in  the  low 
order  solution,  an  indication  of  spurious  oscillations.  When  one  of 
these  regions  is  encountered,  the  low  order  solution  is  used;  in  other 
regions  the  high  order  solution  is  used. 

For  many  problems  the  original  FCT  algorithm  gives  very  good 
results.  However,  because  of  the  extremely  dynamic  nature  of  our  test 
problem,  the  original  algorithm  does  not  give  acceptably  smooth  answers. 
We  have  incorporated  a  more  stable  version  of  FCT  (Book  et  al.11^)  which 
gives  very  nice  results.  Computing  a  smooth  profile  in  smooth  regions  is 
very  important  because  of  the  extremely  strong  dependence  of  nucleation 
rate  on  temperature. 

As  mentioned  earlier,  the  equation  of  state  requires  density 
and  temperature  as  input.  Unfortunately,  transport  simulations  give 
density  and  internal  energy,  but  not  temperature  directly.  Thus,  one 
must  iterate  within  the  equation  of  state  to  find  the  temperature  to 


which  the  internal  energy  corresponds.  This  takes  a  lot  of  computer 
time.  Instead  of  iterating,  we  have  formed  a  large  table  of  pressure  as 
a  function  of  density  and  internal  energy.  Due  to  the  complexity  of  the 
equation  of  state,  interpolating  from  the  table  is  much  faster  than 
direct  evaluation,  even  if  the  temperature  were  known.  We  have  used  a 
shape  preserving  interpolation  method  developed  by  Delbourgo  and 
Gregory  ^  1  ®  to  interpolate  accurately,  and  to  reduce  the  number  of  data 
points  required  . 

c .  Results  . 

We  have  computed  the  solution  to  our  test  problem  for  a  number 
of  initial  energies  and  droplet  size3.  The  hydrodynamic  results  are 
given  for  a  droplet  initially  heated  by  a  laser  pulse  to  2400  K  exploding 
into  an  ideal  diatomic  gas  at  300  K  and  one  atmosphere.  The  nucleation 
rates  that  result  for  droplet  initial  conditions  of  2400  K  and  1600  K  are 
also  presented. 

Figures  B  — 1  — 1  through  B-1 -4  show  a  time  series  of  density, 
pressure,  temperature,  and  velocity  respectively  for  a  one  micron  droplet 
initially  heated  to  2400  K.  The  spherically  symmetric  explosion  of  a 
hot,  dense  gas  has  been  well  studied  because  of  the  interest  in 
explosives.  Brode11^  has  given  a  very  good  description  of  the  structure 
and  behavior  of  the  flow,  to  which  we  refer  anyone  desiring  a  detailed 
explanation  of  the  hydrodynamics.  Briefly,  the  system  manifests  two 
shock  waves.  One  shock  wave  exists  in  the  water  region  and  faces  toward 


the  origin.  The  second  shock  exists  in  the  air  and  fanes  -uiv  fr<-r 
origin.  Ihe  contact  discontinuity,  which  is  the  remnant  ot  -  h  m 


discontinuity,  separates  the  water  from  the  air. 
discontinuity  is  characterized  by  a  discontinuous  change  i- 

temperature,  and  continuous  pressure  and  velocity  profiles.  r't 


are  discontinuous  in  all  of  the  flow  variables.  The  region  bet  wee-  *n-' 
water  shock  and  the  origin  is  a  region  of  strong  expansior. . 

Notice  the  very  high  temperatures  to  which  the  surrounding 


ideal  gas  is  heated  by  the  passage  of  the  shock.  The  lowering  of  tne 
breakdown  threshold  in  air  in  the  presence  of  fine  particles  is  well 
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knovm.  Steverxj  ing  c 1  has  suggested  that  one  of  the  possible  mechanisms 


for  this  phenomenon  is  shock  heating  of  the  air.  We  do  not  consider 
breakdown  here.  However,  our  results  support  the  possibility  of 
ionization  in  the  shock  heated  air.  Substantial  ionization  can  occur  at 


the  high  temperatures  found  in  the  simulation.  Avalanche  breakdown  could 
subsequently  occur  in  the  presence  of  intense  EM  R.  It  must  be  remarked, 
nowever,  that  at  high  fluences,  breakdown  by  ionization  processes  in  the 


droplet  phase  may  also  occur. 

The  region  between  the  origin  and  the  rear  facing  (water)  shock 
undergoes  extremely  strong  isentropic  expansion  before  colliding  with  the 
near  ric.ng  shook,  TV:e  expansion  causes  very  rapid  cooling  of  the  water 
vapor  in  the  region.  As  the  expansion  proceeds,  the  equation  of  state 
begins  to  predict  nonphysical  (negative)  pressures  in  this  region  of  tne 


flow. 


indicates  that  the  Euler  equations  are  no  longer  sufficient 


The  nonphysical  results  appear  because  the  water  cools  to  well 
below  the  point  at  which  it  should  condense.  The  Euler  equations  make  no 
provision  for  a  phase  change  in  the  fluid.  Thus,  as  the  expansion 
proceeds  the  water  vapor  becomes  more  and  more  supersaturated  until  the 
model  breaks  down. 

By  acting  strangely,  the  equation  of  state  is  saying  that  all 
is  not  well  with  the  simulation.  This  is  the  most  compelling  reason  to 
use  a  realistic  equation  of  state.  Under  most  conditions,  the 
differences  between  water  vapor  and  an  ideal  gas  are  small,  at  least  not 
large  enough  to  drastically  change  the  hydrodynamics.  However,  the 
difference  between  condensing  water  vapor  and  an  ideal  gas  is  very  great. 
Condensation  can  therefore  affect  the  hydrodynamics  (and  the 
electrodynamics)  cond  siderably.  Simulating  the  physical  properties 
realistically  often  indicates  when  and  where  additional  physics  becomes 
important.  This  behavior  might  otherwise  not  be  identified  and  cause 
incorrect  results. 

As  a  first  step  towards  incorporating  the  effects  of 
condensation  on  the  flow  we  have  computed  the  nucleation  rate  in  the 
water  vapor  as  a  function  of  time  and  space.  Figures  B-l  -5  and  B-1  -6  show 
nucleation  rates  for  initial  temperatures  of  2400  K  and  1600  K 
respectively.  The  logarithm  of  the  nucleation  rate  as  a  function  of 
position  is  shown  for  various  times. 

Comparison  of  Figure  B-1 -5  with  Figures  B-1  -1  through  B-1  -4 
shows  tfat  nucleation  first  develops  in  the  region  of  lowest  density  and 


temperature.  This  is  the  region  that  has  undergone  tne  greatest  amount 
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Figure  B-1-5.  Time  series  of  the  logarithm  of  nucleation  rate  versus 
radius  for  an  initial  droplet  temperature  of  2400  K. 

of  expansion  and  cooling.  As  time  proceeds,  the  nucleating  regions  grow 
toward  the  origin.  The  nucleation  rates  start  out  small  but  grow  to  very 
large  values  as  the  expansion  proceeds  and  the  vapor  becomes  more 
supersaturated.  The  nucleation  rates  continue  to  increase  because  there 
is  no  mechanism  to  relax  the  high  supersaturation.  Figure  B-1-6  shows 
that  ultimately  the  super  sa  turation  grows  so  large  that  it  drives  the 
exponential  term  in  the  nucleation  rate  to  one  and  the  rate  reaches  its 
maximum  possible  value  (according  to  classical  nucleation  theory).  This 
point  is  nonphysical  and  under  these  conditions  the  equation  of  state 
predicts  extremely  low  and  sometimes  negative  pressures.  The  nonphysical 
regions  can  be  avoided  only  by  incorporating  the  nucleation  and  growth 
equations.  These  would  allow  the  vapor  to  condense  and  reduce  the 
excessive  supersaturation.  Note  that  whereas  classical  nucleation  theory 
is  sufficient  to  describe  when  and  where  nucleation  should  occur,  there 
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Figure  8-1-6.  Time  series  of  the  logarithm  of  nucleation  rate  versus 
radius  for  an  initial  droplet  temperature  of  1600  K. 

is  some  doubt  as  to  its  ability  to  correctly  predict  nucleation  rates  in 
this  problem.  Since  the  flow  timescale  is  so  short,  the  onset  of  high 
supersa tur3 tion  occurs  very  quickly.  Establishing  a  "detailed  balance" 
between  evaporation  and  condensation  rates,  a  requirement  of  the 
classical  theory,  may  not  be  a  realistic  assumption  in  this  case. 

Strong  nucleation  occurs  after  only  a  brief  expansion  for  ail 
of  the  initial  temperatures  considered,  from  1600  K  through  2800  K.  As 
one  might  expect,  for  lower  initial  temperatures,  nucleation  becomes 
important  earlier  and  the  system  expands  less  before  the  simulation 
becomes  nonphysical .  It  seems  that  nucleation  and  particle  growth  should 
be  considered  under  almost  all  conditions  of  strong  convection. 

For  an  initial  condition  of  2400  K,  the  nucleation  rate  .'  ..us 
iorurti/  i  r  cm  its  maximum  value  in  the  field  to  zero  when  the  ex;  v\. 

.•  1  :d  es  with  tne  shock  wivt  ind  is  re  he  ate  1 .  In  this  reg  ion  ' 


particles  foraed  previously  will  probably  be  vaporized. 


For  the  ’.over 


initial  temperature  '  1  600  K)  ,  the  temperature  rise  at  the  shock  or.ly 
decreases  the  value  of  the  nucleation  rate,  which  still  remain"  quite 
high.  This  results  in  the  "arm"  in  the  profiles  of  Figure  3-1  .  The 

nucleation  rate  goes  to  zero  only  at  the  contact  discontinuity. 

Chang^k  has  reported  observing  what  appear  to  be  jets 
emanating  from  single  droplets  upon  intense  laser  irradiation.  The  jets 
originate  at  hot  spots  on  either  side  of  the  droplet,  along  the  axis  of 
the  beam,  and  have  a  cloudy  appearance.  The  droplet  seems  to  maintain 
its  basic  shape  and  does  not  appear  to  be  shattering.  This  observation 
indicates  that  the  jets  may  be  composed  of  a  hot  stream  of  vaporized 
material  whTch  subsequently  condenses  to  form  fine  particles.  Scattering 
of  light  by  the  particles  could  explain  the  cloudy  appearance  of  the 
jets.  Thus,  convection  dominated  hydrodynamics  coupled  with  nucleation 
and  particle  growth  provide  a  possible  explanation  for  this  phenomenon. 

There  are  several  areas  that  need  to  be  considered  to  assure 
that  the  hydrodynamic  simulation  is  realistic. 

As  can  be  seen  from  Figures  B-1  -1  through  B-1-4,  the  relevant 
timescale  for  large  charges  in  the  flow  variables  is  on  the  order  of  0.1 
nanosecond.  Ir.  order  that  no  appreciable  motion  occurs  during  absorption 
of  energy,  the  length  of  the  laser  pulse  must  be  on  the  order  of  0.01 
nanoseconds  for  a  i  micron  droplet.  From  the  equation  of  state,  the 
energy  required  to  raise  the  droplet  from  ambient  conditions  to  2400  K  is 
about  6  x  i  ' h-  f/kg.  T.n  i  s  ir  r  1  i»  •  that  the  droplet  must  absorb  tne 
energy  at  a  rate  of  about  6  x  *  O'  1  W/  kg  over  the  heating  ;°ri>d. 


Assuming  that  the  optical  properties  of  ilu  water  droplet  remain  constant 


during  the  heating  period,  the  intensity  required  to  accomplish  the 
heating  (see  Armstrong'^  or  pend  le  ton 1  ^3)  iS  on  the  order  of  1011  -  lO1^ 
W/cm’.  The  breakdown  potential  for  clean  air  is  on  the  order  of  10^ 
W/cm^.  This  indicates  that  in  practice,  heating  a  droplet  to  iii/h 

temperatures  before  substantial  motion  occurs  ana  in  the  absence  of 
nonlinear  optical  effects  in  either  the  air  or  the  droplet  might  be 
difficult.  Nonlinear  optical  effects  may  dominate  at  these  power 
densities . 

A  strongly  absorbing  droplet  might  allow  the  proper  heating 
rate  without  such  a  high  power  density,  but  this  sidesteps  the  real 
problem.  To  avoid  very  high  power  densities,  and  to  extend  the  range  of 
validity  of  the  simulation,  the  heating  process  must  be  modelled  more 
realistically.  Very  high  power  densities  can  be  avoided  if  the  heating 
process  is  allowed  to  occur  over  a  longer  period.  In  general,  motion 
will  occur  during  the  heating  process,  and  this  will  affect  subsequent 
absorption.  Even  for  a  droplet  that  stays  in  the  Rayleigh  size  limit 
while  expanding,  the  effect  of  varying  temperature  and  density  on  the 
absorption  must  be  accounted.  In  addition,  since  the  vaporization 
process  is  not  even  approximately  quasi-stationary,  the  kinetics  of  the 
phase  change  should  be  considered. 

3y  requiring  instantaneous  heating  of  the  droplet,  the  current 
simulation  predicts  the  most  dynamic  response  possible  for  a  given  input 
energy.  Thus,  it  traces  an  envelope  of  limiting  behavior.  Simulations 
with  a  more  realistic  treatment  of  heating  should  fall  within  the 
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Finally,  since  the  relevant  timescale  for  flow  is  so  short,  the 
partitioning  of  energy  amongst  the  translational,  rotational,  and 
vibrational  modes  of  the  molecule  should  be  considered.  On  one  hand,  the 
absorption  of  energy  from  the  laser  will  cause  a  departure  from 
equilibrium  amongst  the  various  modes.  The  molecules  will  move  toward  an 
equilibrium  state  through  collisional  relaxation.  If  the  time  for  this 
relaxation  is  of  the  same  order  as  the  convective  timescale,  then  the 
initial  nonequilibrium  state  will  affect  the  hydrodynamics. 

Nonequilibrium  partitioning  of  energy  can  also  arise  as  a 
result  of  the  expansion.  Che  of  the  implications  of  using  an  equilibrium 
equation  of  state  is  that  energy  contained  in  all  modes  is 
instantaneously  available  to  drive  fluid  motion.  However,  depending  on 
the  temperature  and  density,  vibrational  relaxation  times  can  be  longer 
than  the  convective  timescale  of  the  expansion.  When  this  occurs, 
vibrational  energy  will  be  available  only  at  a  rate  determined  by 
collisional  relaxation.  In  order  to  account  for  these  effects 
explicitly,  the  transport  equations  need  to  be  reformulated  in  terms  of 
the  individual  molecular  energy  modes  (see  Anderson10®). 

d.  Conclusions. 

We  are  modelling  the  hydrodynamic  response  of  a  small  water 
droplet  to  a  very  large  pulse  of  laser  energy.  Our  current,  simulation 
assumes  isotropic  absorption  and  a  droplet  which  is  inertialiy  locked  in 
place  during  the  absorption  process.  Upon  heating  to  very  high 
temperatures,  the  droplet  explodes.  We  have  solved  the  transport 

1  (I 


equations  that  govern  the  explosion,  incorporating  a  realistic  equation 
of  state  for  water.  The  requirement  of  instantaneous  heating  manifests 
the  msot  dynamic  response  possible  for  a  given  input  energy*  Thus,  the 
current  simulation  represents  a  limiting  case  in  the  hydrodynamic1, 
addition,  we  have  computed  nucleation  rates  for  the  water  v.-.  for 
expands. 


Ln  the  course  of  the  explosion  the  system  manifests  two  ct ,  ./ 
waves,  one  in  the  water  vapor  and  one  in  the  air.  There  is  also  a  reg  ion 
of  strong  expansion  and  cooling  between  the  water  shock  and  the  center  of 
the  droplet.  This  is  consistent  with  previous  studies  of  explosions  in 
curvilinear  coordinates. 

The  air  is  shock  heated  to  temperatures  high  enough  to  oh  un¬ 
ionization.  If  the  ionized  air  is  subjected  to  more  radiation,  avalanche 
breakdown  of  the  air  can  occur.  Our  simulation  substantiates  the  tr.eory 
that  shock  nesting  of  air  car,  be  at  least  partly  responsible  for 
enhanced  breakdown. 

For  all  initial  droplet  temperatures  considered  ( 1 600 
strong  nucleation  occurs  in  certain  region  of  the  flow  after  a  mi 
expansion.  If  condensation  of  the  vapor  is  not  allowed,  the  soluti~’ 
soon  becomes  nonphysical.  Without  using  a  realistic  equation  -f  state 


this  behavior  would  likely  go  unnoticed.  Thus,  in  the  ex 


vaporization  revim.  ,  condensation  appears  to  be  an  important  me,'!-  >r.\ 
which  should  be  considered.  The  presence  of  fine  pai  titles 


appreciably  affect  both  the  hydrodynamics  and  the  electrodynamics 


Finally,  a  realistic  heating  process  need  3  to  be  included. 
Heating  a  water  droplet  to  high  temperatures  before  appreciable  nation 
occurs  requires  very  high  po  wer  densities.  The  power  den  si  ties  c  -r.  be 
lowered  by  relaxing  the  constraint  of  very  little  motion  dur  Lr.v  ire 
heating  process.  However,  this  requires  accounting  for  t.he  variation  of 
absorption  with  temperature  and  density  as  the  motion  occurs ,  still 
assuming  that  the  "droplet"  remains  in  the  Rayleigh  size  limit  during 
absorpt  ior. . 

We  are  currently  incorporating  a  realistic  heating  rate  for  the 
droplet.  In  addition,  we  are  accounting  for  the  possibility  of 
condensation,  in  the  flow  by  including  the  particle  growth  equations. 
These  additions  will  add  a  new  dimension  of  realism  to  the  simulation, 
yielding  information  that  can  be  verified  experimentally  and  providing 
new  insight  into  the  problem. 

3 -? .  Rapid  evaporation  of  droplets  in  laser  beams . 

a.  Introduction, 

There  is  considerable  current  interest  in  the  problem  of  high 
energy  laser  beam  propagation  through  an  atmosphere  containing  an 
absorbing  aerosol.  When  intense  electromagnetic  radiation  interacts  with 
absorbing  particles,  a  variety  of  diverse  phenomena  have  been  noted, 
including  thermal  blooming  (self  steering/defocusiog  of  the  beam'  ,  and 
optical  breakdown.  Fach  of  these  effects  is  due  in  part  to  the 
absorption  and  consequent  comprising  the  aerosol.  Thus  we  first  d  meet 


Two  extremes  of  heating  may  be  identified  for  single  particle 
absorption.  For  highly  intense  beams,  and/or  strongly  absorbing 
droplets,  there  is  the  explosive  vaporization  regime,  in  which  droplet 
deformation  and  shattering  may  occur.  This  regime  is  dominated  by 
convective  effects,  in  which  pressure  gradients  are  responsible  for,  and 
govern,  mass  transfer  from  the  droplet.  At  the  other  extreme,  droplet 
interaction  with  less  intense  beams  falls  into  what  may  be  termed  the 
slow  evaporation  regime.  By  contrast,  this  regime  is  dominated  by 
diffusive  effects,  in  which  concentration  driving  forces  determine  mass 
transfer  rates  from  the  particle.  Carls  and  Brock1^  have  recently 
studied  the  case  of  droplet  explosion,  shock  wave  formation,  and 
expansion.  Here  we  study  the  case  of  slow  evaporation,  that  is, 
evaporation  in  which  the  droplet  temperature  never  exceeds  its  boiling 
point,  and  in  which  a  constant  pressure  field  is  assumed  to  prevail 
around  the  droplet. 

Under  this  restriction  of  constant  pressure,  we  treat 
rigorously  the  time  dependent  continuum  transport  equations  which 
describe  the  hydrodynamics  of  single  droplet  evaporation,  and  develop  a 
computer  code,  !  K!>r  >  1  -niKir  : .  !<•!■!  i.i  I  ct  :  \ ■  •<  ■  r.i  i  i .  a  ■  .  for  their 

numerical  solution.  A  slip  flow  analysis  is  introduced  to  account  for 
noncontinuum  effects.  We  compare  our  solution  with  a  quasi-static 
approximation  to  the  evaporation  process  which  has  been  widely  used  in 
the  past, 1  1 35,  1 36  arKj  so  examine  the  validity  of  that  approximation. 

Further,  from  a  unified  treatment  of  the  droplet  mass 


ml  energy 


balances,  we  examine  the  consequences  of  various  simplifying  assumptions 
that  have  often  been  implicit  in  the  past. 

b.  Theory^ 


The  physical  system  to  be  described  consists  of  an  isolated 
spherical  droplet  (or  particle)  in  an  effectively  infinite  host  gas 
illuminated  by  laser  radiation.  Hie  amount  of  heat  an  individual  droplet 
absorbs  depends  on  the  particle's  complex  refractive  index  and  its  Mie 
size  parameter  (the  ratio  of  the  particle's  circumference  to  the 
wavelength  of  the  incident  radiation).  We  have  utilized  Pendleton's1^ 
code,  PENDSOR,  to  determine  the  Mie  size  parameter  below  which  the 
droplet  heats  uniformly.  PENDSOR  computes  a  dimensionless  heating  source 
function  along  specified  chords  of  a  sphere.3  For  droplet  size 
parameters  below  1.0,  PENDSOR  shows  little  variation  in  source  function 
along  a  diameter.  Thus,  for  incident  radiation  at  10.6  um,  restriction 
of  the  analysis  to  particles  which  heat  uniformly  places  an  upper  limit 
on  particle  radius  of  -2.0  um.  This  restriction  makes  the  problem 
spherically  symmetric.  Extension  of  the  model  to  larger  droplets  which 
do  not  absorb  uniformly  would  require  a  solution  of  both  the  Maxwell 
equations1^  and  the  transport  equations  *39  over  the  droplet  interior. 


aPendleton  defines  a  dimensionless  source  function,  or 
normalized  power  density  as: 


s  .  :s/Elnc;2 


where  E  is  a  complex  representation  of  the  electrical  field  induced  by 
the  laser  within 
electrical  field . 1 


droplet. 


magm! 


of  the  incident 


For  an  incompressible  droplet  of  density  .  p,  and  constant  neat 
capacity  Cp  ~  cv,  the  general  energy  conservation  equation  at  my 
interior  point  has  the  form:1^ 

p  ■/  t(  CppT^+1  /2v^)  +  ip.  *[(CpTp  +  P  / .  p  +  1  /2v^)  v]  -  -  ■  ’  ;  q+  q  r 

where  Tb  is  the  droplet  temperature  at  any  point  (r,',,‘)  within  the 
droplet,  and  v^  the  velocity  at  that  point.  The  conductive  heat  flux  i: 
the  droplet  surface,  q,  is  given  by  a  -Kp-Tp,  where  Kp  is  the  thermal 
conductivity  for  the  drop,  qr  accounts  for  the  laser  heating  source.  We 
consider  only  a  lumped  parameter,  macroscopic  description  of  droplet 
behavior  in  terms  of  an  average  droplet  temperature  Tp.  This  is  achieved 
by  integrating  Equation  ( B  —2  —  1 )  term  by  term  over  the  droplet  volume, 
with  the  definition  Tp  =  ,'vTLdV/  -"vdV,  where  V  is  the  droplet  volume. 
Utilizing  the  spherical  symmetry  of  the  problem,  application  of  GAuss' 
divergence  theorem  yields. 

CpDd/dt  (4/3  R^.pTp)  -  4-R2KpPTL/V)iR  +  -R2Qa5sI0  ( E-  - 

where  terms  representing  the  internal  kinetic  energy  and  the  work  due  tc 
droplet  expansion  have  been  a  glee  ted.  The  droplet  is  irradiated  with  a 
uniform  beam  of  intensity  Iq,  and  we  account  for  the  laser  heating  source 
with  a  Mie  absorption  efficiency 138  Qabs.  For  water  droplets  below  in 
in  radius,  Qabs  i 3  a  linear  function  of  radius.13®  Hence  the  1  .ser 
heating  term  can  be  written  as  4/3"  R  3  J  q,  where  is  a  ?on:-’  •»  r.  t  , 
bulk  absorption  coefficient. 


Consistent  with  the  macroscopic  description  of  droplet 
behavior,  we  donot  calculate  the  internal  droplet  temperature  profile  to 
evaluate  the  surface  flux  term,  Kq(  'T-^/  r)  I  p.  Instead,  the  term  is 
replaced  through  the  use  of  an  interfacial  energy  balance.  Equating  the 
heat  flux  from  the  droplet  interior  at  the  surface  (conductive  only)  with 
that  leaving  the  surface  (conductive,  kinetic,  and  convective  energy 


fluxes,  and  latent  heat  loss)  to  the  surrounding  host  gas  yields: 


U  H2{-Kgf  Tg/  r)iR  +  (■  vCpgTg)  !  R  +  1/2.  avad} 
HvapdM/dt  (B-2-3) 


where  T\.  t,r)  is  the  host  gas  temperature  at  any  point  outside  the 
droplet,  K,  is  its  thermal  conductivity,  and  Cn~  its  heat  capacity.  is 
the  total  mass  concentration  of  the  gas,  and  .  ,  the  mass  concentration 

a 

of  the  vapor.  va  is  the  radial  outward  directed  velocity  of  the  vapor 
species,  ana  v  the  radial  mass  average  velocity  of  the  gas.  Substitution 
of  Equation  (3-2-3)  in  Equation  (B-2-2)  and  rearrangement  yields: 


DCpDU/3  -R3dTD/dt  +  H  R2- DCpDTDdR/dt  =  4  R2 {Kg(  -Tg/  r)  1  R 

-  ■  ■  vCpgTg)  ’  R  -  1/2  ava3)  ♦  HvapdM/dt  +  4/3  R3  IQ  (B-2^4) 


Comparison  of  Equation  (B-2-4)  with  previous  analyses  1 3^’ 1  3^  shows  two 
hitherto  neglected  terms:  one  representing  bulk  convection  of  energy 
away  trom  the  interface  by  the  moving  host  gas,  (.  vrpgTg!R),  and  a  term 
which  accounts  for  the  droplet  radius  shrinking  with  time. 


We  will  examine  the  significance  of  these  terms. 


We  derive  a  droplet  mass  balance  taking  into  account  the  fact 
that  the  surface  is  shrinking  with  a  velocity  dP/dt.  Mass  flux  of  the 
vapor  species  relative  to  the  moving  boundary  is  given  by: 

1  /( 4~  r2)  dM/dt  =  -  .  a(va  dR/dt)  (B-2-5) 

where  M(  t)  is  the  droplet  mass.  Noting  that  ava'  flux 

relative  to  a  stationary  co-ordinate  system,  is  made  up  of  diffusive  and 
convective  components,1^  we  write: 

'  ala  =  ~  Dab  ■  a  +  al  (B-2-6) 

where  a  is  the  mass  fraction  of  the  vapor  species  (  .  a=  /  .) .  Equation 
'.B-2-6)  requires  a  specification  of  v(v=  •ava  +  .  bvb,  where  vb  is  the 
absolute  veolcity  of  the  inert  component  of  the  host  gas).  The  inert 
component  is  assumed  to  be  insoluble  in  the  droplet,  and  hence  at  the 
interface,  v  b  =  dR/dt.  The  velocity  of  the  water  vapor  species  is 
specified  as  v3-  .  /U  R^+dR/dt,  where  .  is  the  volumetric  flowrate  of 
water  vapor  from  the  surface.  These  definitions,  with  the  assumption 
that  the  host  gas  behaves  ideally  and  with  M(  t)  =  (4/3)  P-(t)  D,  lead  to: 

v  =  dR/dt  (.  -  D)  /.  (3-2-7) 

substituting  (B-2-6)  and  (B-2-7)  in  (B-2-5),  and  solving  for  dR/dt: 

dR/dt  =  ■  CAB  a/  [  D(1-  a)]  '3-2-9) 

Tne  numerator  represents  normal  Fickian  diffusion  from  a  stationary 
surrj'e.  The  denominator  accounts  for  mass  transfer  from  the  surface  by 
convection,  i  .e  .  bulk  movement,  of  the  host  gas  at  the  sur  f  i  'e . 

Equations  ( B  -2  -4)  and  (B-2-8)  are  two  ord  ir.  -.ry  differential 
equations  which,  with  the  initial  conditions  zT-  ,  and  R(  tsO)  =  R(y 

describe  the  droplet  behavior  over  its  1  i  f«t.  ime.  ••••  w ,  tne  energy 


conservation  Equation  .  B-2-4)  requires  a  specification  of  the  host  gas 
energy  flux,  Kg(  -Tg/  -r)  !  p,  at  the  surface,  and  the  mass  conservation 
Equation  (3-2-8)  requires  a  mass  flux,  n These  fluxes  can  only 
be  obtained  by  a  simultaneous  solution  in  time  of  the  equations  of 
continuity,  energy,  and  motion,  which  govern  respectively,  the 
concentration,  temperature,  and  velocity  profiles  in  the  atmosphere 
surrounding  the  drop. 

We  limit  our  simulation  to  the  case  of  an  isobaric  gas  around 
the  droplet,  with  the  argument 1 ^  that  any  pressure  wave  in  the  host  gas 
induced  by  evaporation  at  the  surface  is  a  transient  which  propagates  at 
the  local  speed  of  sound.  Thus,  for  the  time  scale  in  which  we  are 
interested,  pressure  equilibration  is  complete.  In  contrast  to  earlier 
stud  ies, 1 35, 136,  1  ^0  the  mediura  surrounding  the  droplet  was  not  assumed  to 
have  a  constant  density.  While  the  assumption  of  an  incompressible  fluid 
is  appealing  in  that  only  the  isobaric  case  is  considered,  temperature 
variations  of  up  to  30%  will  cause  similar  vairations  in  density.  With 
these  restrictions,  the  transport  equations  in  the  host  gas  surrounding 
the  spherical  droplet  are  well  known.  ^9 


DAB(T)  =  Dabo  (T/T0)  ( B -?-!?:• 

Kg(T)  =  KgO  (T/T0)  (  B -1  it' 

where  the  subscript  ’O'  denotes  the  reference  state.  Shear  stress  shown 
is  for  the  Newtonian  fluid,  with  viscosity,  •  Although  the  host  gas  : 
assumed  ideal,  the  inclusion  of  a  more  rigorous  equation  of  :  i- 

straightforward.  The  system  (B-2-9)  -  (B-2-12)  is  highly  coupled  an  i 
nonlinear,  and  can  only  be  solved  numerically. 

For  purposes  of  numerical  calculation  we  consider  a  spherical 
water  droplet  in  an  initially  ( t  =0)  uniformly  saturated  atmosphere  of 
water  vapor  in  air,  stationary  relative  to  the  sphere.  At  t=0  tne 
droplet  is  exposed  to  a  beam  of  intensity  Iq.  The  initial  conditions  for 
the  system  are  thus  : 

T(  r  ,t  =0 )  =  Tdo;  V(r,t=0)  =  0;  ■  a(  r,t=0)=  •  a0  ’  B -2  - 1 R } 

The  interior  interfacial  boundary  conditions  in  temperature  ar.d 
concentration  are  given  by  solution  to  the  ordinary  differential 
equations  ( B  —2  —4 )  and  (B-2-8).  Departure  from  continuum  behavior  is 
accounted  for  through  the  introduction  of  concentration  and  temperature 
j urn p  cond i t ions , 1 ^  appro pr ia te  for  the  si  i  p  reg ime . 1  -  -  For  exterior 
interfacial  conditions  the  medium  is  assumed  to  be  undisturbed. 
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and  C-j  are  respectively  the  temperature  and  conne.ntrat  i "  r. 
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motion.  “jss  concentration  at  the  sphere  surface,  a3,  i-3  calculated 

using  Antoine's  expression  for  vapor  pressure,  with  the  Kelvin  surface 

,  1 1 1? 

curvature  correction  included. 

The  system  of  Equations  (B-2-9)  through  (B-2-lt)  is  highly 
coupled  and  nonlinear,  but,  with  the  restriction  of  isobaric  evaporation 
is  parabolic  in  nature.  The  boundary  conditions  to  Equations  (B-2-9) 
through  (3-2-11)  are  specified  by  the  time  integration  of  Equations  (E-2- 
A'1  and  (B-2-9).  The  system  (B-2-9)  through  (B-2-11)  presents  a  moving 
boundary  (Stefan-type)  problem,  since  the  sphere  interface  recedes 
with  time.  Although  the  actual  radial  domain  of  interest  is  semi¬ 
infinite,  the  computational  doraan  has  been  arbitrarily  limited  to  100 
droplet  radii.  Results  from  the  simulation  show  the  host  gas  undisturbed 
at  TO  droplet  radii  and  validate  this  assumption.  When  solving  the 
transport  equations  for  the  host  gas  we  have  chosen  to  ignore  the  moving 
boundary  problem  entirely,  since  complete  disappearance  of  the  droplet 
constitutes  only  a  1?  change  in  the  computational  domain.  For  each  of 
the  partial  differential  equations  (3-2-9)  through  (B-2-11),  3  simple 
implicit  finite  differencing  scheme  was  used.  All  coupling  terms  and 
non!  inear  it  ies  have  peer,  evaluated  explicitly  at  the  previous  time  step. 

.  SAim-r ;  Aes.ilt?  and  Discussion. 

The  dual  r.ians  have  been  limited  to  the  case  of  a  pure  water 
droplet  irnlii'el  with  monochromatic  laser  light  at  1  Oo?5  m. 
v*  n  s  i '  n  \  f  oh*'  -  y.  -i  to  include  droplets  comprising  a  non  ideal  aqueous 


solution  requires  a  stra  ight  for  ward  modification  of  the  vapor  pressure 
calculation.1^  Table  B-2-1  lists  the  physical  property  data  used. 

In  this  section  we  first  discuss  the  importance  of  terms 
previously  neglected  in  the  droplet  energy  balance.  We  then  examine  the 
effect  on  the  solution  of  the  transport  equations  of  both  a  temperature 
dependent  diffusion  coefficient,  and  of  the  temperature  and  concentration 
jump  boundary  conditions.  Next  we  examine  how  laser  energy  incident  on  a 
particle  is  distributed  in  various  forms  over  the  particle's  lifetime. 
Finally,  the  transient  solution  described  is  compared  with  a  quasi-steady 
approximation  to  the  evaporation  process.  The  limitations  of  this 
pseud  o- ste  ad  y  approach  are  noted,  and  we  show  how  its  accuracy  and 
applicability  may  be  improved  in  several  ways:  through  the  use  of  a 
rigorous  energy  balance,  a  temperature  dependent  diffusion  coefficient, 
and  a  proper  incorporation  of  the  slip  flow  boundary  conditions. 

We  start  by  assessing  the  relative  importance  of  four  different 
effects  on  droplet  heating.  These  include:  (a)  accounting  for  the 
shrinking  droplet  radius  with  the  term  4  R?  pCpnTDdR/dt  in  'Equation  B-2- 
4);  (b)  including  convective  energy  losses  from  the  droplet,  with  the  term 
(  vCp^Tgjlp  in  (Equation  B-2-4);  (o'1  accounting  for  the  temperature 
dependence  of  the  diffiusion  coefficient  with  Equation  (B-2-1  (d) 

extending  the  continuum  analysis  in  an  approximate  rnitter  to  the  slip 
flow  regime  with  the  jump  boundary  conditions,  Equation  'B-2-1  t).  None 
of  these  effects  have  been  accounted  in  previous  anal yses. 1  1  5 4  %  1 144 

The  individual  contribution  of  each  of  these  effects  t  •  d r'q-let  behavior 


Table  B-2-1.  Physical  Constants  Used  in  the  Simulation  (at  298  K) 


Symbol 


Value 


Units 


53. 

J  kg  ~ 1 K  ~ 

>8. 

J  kg^K-1 

1.5 

d imensionless 

2.3 

dimen  sionless 

2.6x10-5 

m^sec-1 

2.259x106 

J  kg-1  1 

0.025 

Jm-1sec-1K-1 

0.02265 

0*4 

3 

1 

UJ 

'0 

kg  m“5 

18 

K 

18 

K 

0.019 

dimensionless 

7.87x10u 

m_1 

1.786 

dimensionless 

1.5 

d  imensionless 

6.6x10“® 

m 

1 . 8xl0~5 


kgm-1sec-1 


was  determined  by  'switching  off1  other  terms  in  the  simulation. 


B-2-1  and  B -2-2  compare  these  effects  .in  the  long-term  heatm;  ce:«*v  1 
of  a  droplet  initially  1  .m  in  radius.  In  these  figures  eaen  of  *w.v- 
above  effects  is  compared  with  a  base  line  (curve  A),  in  wr.  i  c  h  -v. : 
conductive  and  latent  heat  losses  from  the  droplet  were  cor. >1  ’ 
i"  which  diffusion  coefficient  D^g  was  held  constant,  while  the  si  . :  f. 
corrections  to  the  boundary  conditions  were  ignored.  In  Figure  r  -2  - 1 
which  depicts  droplet  temperature  variation  with  time,  and  in  Figure  B-' 
2,  which  shows  howdroplet  radius  varies  with  time,  each  of  the  c  u  rv  e  > 
shown  represents  the  simultaneous  solution  of  Equations  (B-2 -4,  and  ,  B-2  - 
3)  through  (B-2-16).  Although  the  droplet  temperature-time  curves  i~ 
Figure  B  —2  - 1  all  vary  in  their  prediction  of  peak  droplet  temperature 
they  are  similar  in  that  none  show  a  steady-state  droplet  temperature 


Similar  behavior  has  been  reported  by  Sageev  and  Seinfeld, 


1  36  . 


an; 


> 


Caledonia  and  Teare.7^0  We  have  modified  the  simulation  so  as  1'  ;oi. 
the  droiet  size  constant,  and  observed  that  droplet  temperature  increases 
mono  ton  icaliy  with  time,  asymptotically  approaching  a  steady  state  -v 
which  point  heat  losses  by  evaporation,  conduction,  etc.,  balance  the 
laser  heating  source.  Thus  the  nonstationary  temperature  is  1  k  sc 
to  the  fact  that  the  droplet  is  shrinking.  Th  i  s  suggests,  an  "rro 
previous  analyses  which  have  assumed  that  droplets  rise  quick:. 
steady  temperature ,  at  which  they  remain  while  ev  ape  rat  i  nv  ,J'' 

Figure  3  -2  -i  makes  it  clear  that  failure  to  accouh  1 
rat  s  '  f  change  of  droplet  shrinkage  (  the  term  f  r  ,  :  .  >,  g  c, 

-  ■  will  result  In  a  significant  underpred  in  tier,  o  t.c-  .  • 


Figure  3-2-1. 


alOC’C;TIU£)  (S£C) 

Droplet  tern  pera t ure  behavior  in  response  to  long  time 

heating  for  a  1  pm  drop.  Iq=1  0^  W/m^. 

(A)  Conductive  and  latent  heat  losses  accounted  for 

( B)  Curve  (A)  with  the  temperature  dependence  of  DAB 
included . 

'  C)  Curve  (A)  with  convective  heat  loss  from  drop  includ¬ 
ed  . 

(D)  Curve  (A)  with  droplet  mass  loss  included. 

(E)  Curve  (A)  with  the  incorporation  of  slip  flow  boun¬ 
dary  conditions. 

(  F)  True  solution,  i  .e  .  Curve  (A)  with  all  above  effects 
included  . 
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t.  siz-^  behavior  in  response  to  long  time  neat  mg  - 
.ions  n  s  in  Figure  8  - 1 .  Note  that  curves  ■  B'  and 

•.distinguishable  from  curve  '  j\]  , 
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temperature,  and  a  consequent  overpr  ed  ic  tion  of  its  lifetime  (Figure  B-2- 
2).  In  the  context  of  fog  'hole-boring'  it  may  be  desirable  to  stay  in 
the  slow  evaporation  regime  while  burning  through  a  cloud  in  order  to 
avoid  the  shock  formation  and  refractive  index  changes  which  might 
accompany  more  rapid  evaporation  and  perhaps  droplet  explosion.  These 
effects  may  serve  to  further  scatter  and  attentua'e  any  ensuing  laser 
beam.  Thus,  an  accurate  prediction  of  droplet  peak  temperature  is 
essential  when  one  wants  to  determine  the  onset  of  droplet  boiling  and 
perhaps  explosion,  or  to  study  the  thermal  blooming  effects  caused  by 
highly  nonuniform  temperature  fields  in  the  atmospheric  medium.  Although 
not  shown  on  these  plots,  kinetic  energy  losses  and  accounting  for  the 
temperature  dependence  of  have  contributed  negligibly  to  the  droplet’s 
behavior.  Inclusion  of  convective  heat  losses  (curve  C)  from  the  droplet 
surface  has,  as  might  be  expected,  lowered  its  peak  temperature. 
Accounting  for  the  temperature  dependence  of  DAg  (curve  B)  similarly 
lowers  the  peak  droplet  temperature.  Curve  E  in  Figures  B-2-1  and  B-2-2 
shows  the  effect  of  extending  the  continuum  bound  iry  conditions  (obtained 
with  Ct=Cd=0)  to  the  slip  flow  regime.  liie  temperature  discontinuity  at 
the  droplet  surface  has  served  to  lower  the  effective  driving  force  for 
evaporation,  hence  causing  the  droplet  to  rise  to  a  higher  peak 
temperature  in  order  to  continue  evaporating  at  the  same  rate  (Figure  B- 
2-2).  Figures  B-2-1  and  B-2-2  make  it  clear  that  while  trie  effects  of 

convective  energy  loss  from  the  droplet,  temperature  dependent  D*b>  and 
slip  flow  boundary  conditions,  are  of  lesser  importance  than  accounting 


for  the  rate  of  irople 
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Next , 

to  further  examine 

ho  w 

energy  is 

dist'ibuted  in 

various 

forms  within  the 

droplet,  Equation 

(  B-2  - 

4'  is 

rewr 

itten  : 

eH  +  Ek'E  +  ECV  + 

eCD  + 

E  VAP 

=  Et 

*  eh  < 

B-2-17) 

in  which: 

EH  =  c 

pD'  D4/3  R3dVdt 

eke 

=  1 /2A  ava3 

Ecd  = 

-AKg(  Tg/  r)  ;  R 

ECV 

-  v Epg^'1  1  R 

Evap  = 

~HvapdM/dt 

el 

=  4 /3  R  3  Iq 

Er  = 

-A-  DCp DTDd  R/dt 

A 

=  4-R2 

Terms  on  the  left  hand  side  of  Equation  (B-2-17)  represent  modes  in  which 
energy  is  dissipated  in  the  drop,  i.e.:  droplet  heating,  kinetic  energy, 
convective,  conductive,  and  latent  heat  losses.  Terms  on  the  right  hand 
side  represent  the  laser  source  of  heat,  and  the  ’artificial  source’  due 
to  the  shrinking  radius.  Each  term  in  Equation  (B-2-17)  is  normalized  by 
the  right  hand  side,  yielding: 

+  ^KE  +  qCD  +  QCV  +  ^VAP  =  1  (B-2-18) 

in  which  Qx  =  EX/(ER+EL).  This  formulation,  in  terms  of  the  fractional 

energies,  provides  a  convenient  means  of  studying  energy  distribution 
in  an  absorbing  droplet.  We  have  deliberately  paralleled  Armstrong's 
analysis  here  in  order  to  allow  an  immediate  comparison  with  his 
results.1^  Figures  B-2-3  and  3-2-4  show  these  fractional  energies  for 
both  small  drop/ low  intensity,  and  large  drop/high  intensity  extremes. 
Initially  the  droplet  is  at  equilibrium  at  the  ambient  temperature,  and 
hence  cannot  lose  energy  by  conduction,  evaporation,  etc.;  thus  100$  of 


V- 

......... 


.  i. 


.c.  -iUV.  JVjV. 


the  source  energy  goes  into  droplet  heating.  As  the  drop  heats  it  li¬ 
able  to  lose  energy  through  the  conductive,  evaporative,  and  convective 
mechanisms  discussed,  whose  rates  increase  with  increasing  temperature. 
Figure  8-2  -3  shows  latent  heat  losses  asymptotically  approaching  of 

the  source,  conductive  losses  -20$,  and  convective  losses  -101.  For  this 
extreme  of  small  droplet/low  intensity,  similar  results  have  been 
reported  by  Armstrong,1^  shows  asymptotic  limits  of  72$  for  latent 

heat  losses,  and  28$  for  conduction  (convective  losses  were  not  accounted 
for  in  his  analysis).  For  the  opposite  extreme  of  a  large  drop  subjected 
to  an  intense  beam,  Armstrong1^  has  predicted  that  latent  heat  losses 
will  asymptotically  absorb  100$  of  the  source  energy,  and  conduction  0$. 
However,  due  to  the  nature  of  the  approximation  used  in  his  solution, 
Armstrong  was  unable  to  verify  this  result.  When  we  ignore  convective 
losses  from  the  droplet  (the  term  Efv  in  Equation  B  —2  —17)  in  our 
simulation,  although  not  shown  we  find  exactly  this  behavior.  At  this 
extreme,  the  droplet  is  evaporating  at  a  sufficiently  high  rate  that 
conductive  heat  losses,  while  still  present,  are  negligibly  small 
compared  to  latent  heat  losses.  Figure  B -2-4  shows  this  case  when  the 
term  E^y  is  included.  Latent  heat  losses  asymptotically  approach  86$, 
cond uc  tiv elosses  are  insignificant,  and  convective  losses,  following  the 
same  trend  as  latent  heat  losses,  approach  14$.  As  expected,  for  this 
case  of  faster  evaporation,  convective  losses  are  larger  than  those  shown 
in  Figure  3-2-3-  These  high  evaporation  losses  in  comparison  to 
conductive  losses  occur  because  vapor  pressure,  and  hence  ev  a  p  '■  r  a  ’  i  -  r 
rate,  is  an  exponential  function  of  temperature,  wnile  conduct iv---  0  3ser. 


depend  only  linearly  on  temperature.  Note  that  eonvectiv  ->  losses,  unlike 
conductive  losses,  do  not  disappear,  since  they  are  due  purely  to  the  net 
efflux  of  material  from  the  particle.  Figures  B-2-3  and  B -2 -4  should  not 
be  construed  as  implying  that  droplet  temperature  asym  ptotically 
approaches  a  steady  state  (since  the  fractional  energy  used  in 

droplet  heating,  approaches  zero).  They  give  this  appearance  only 
because  the  time  scales  shown  represent  such  a  small  portion  of  the  life 
of  the  drop.  However  these  figures  do  lend  credence  to  the  assumption  of 
a  quasi-steady  state  droplet  temperature  for  heating  times  very  short 
compared  with  the  droplet  lifetime.  Such  a  situation  develops  when  the 
drop  is  exposed  to  a  pulsed  laser  source.  135*136 

Finally,  we  compare  our  transient  model  of  droplet  evaporation 
with  a  quasi-steady  approximation  to  the  process  which  has  been  used  by 
several  previous  inv  est  ig  a  to  rs. 1  3*1,  1  35,  1  36  This  approximate  analysis 
makes  a  great  simplification  in  the  model  development  by  assuming  a 
pse  ud  o-stead  y  state  approximation  to  the  transport  equations  (B-2-9) 
through  (B-2-11)  in  which  the  explicit  time  derivatives  are  ignored. 
This  approach1^  argues  that  the  characteristic  relaxation  times  for 
changes  in  gas  phase  temperature  and  concentration  profiles  are  far 
smaller  than  for  changes  in  the  droplet  temperature.^1  Equations  (B-2-9) 
and  (3-2-10)  thus  reduce  to  two  second  order  ordinary  differential 


°This  is  easily  verified  for  the  air  water  system,  and  laser 
heating  times  of  the  order  of  1  sec.  It  is,  however,  easy  to  imagine  a 
situation  (e.g.,  for  a  short  laser  pulse  (<<  1  sec)  or  for  the  case  in 
which  thermal  d  i  f  f  us  iv  i  t  ie  s  of  particle  and  hose  gas  are  -  equal)  in 
'ft".  Lch  if  m3y  not  be  a  valid  assumption. 
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equations.  Williams  1  '. :  so  us  re.*  'his  approach  in  some  detail  and  shows 

that,  with  the  further  assumption  of  an  incompressible  atmosphere,  with 
constant  binary  diffusion  coefficient  DAg  anda  thermal  conductivity  Kg*  a 
straightforward  analytical  integration  gives  radial  profiles  for  mass 
fraction  .  a  and  host  gas  temperature  Tg»  and  obviates  numerical  solution 
of  the  system  (3-2-1  tr.rough  C  B  — .?  —  1  6 ) .  We  have  generalized  Williams 
approach1 3^  somewhat  by  integrating  the  steady-state  form  of  equations 
(B-2-9)  and  (B-2-’9),  with  the  temperature  and  concentration  jump 
boundary  conditions  of  Equation  (3-2-16).  Di  f  feren  tia  t  i  ng  these  quasi¬ 
steady  profiles  yields  fluxes  of  mass  and  energy: 

■  dAB  '  a7  r'  R  =  dM/dt  (  a  -  a)  /(<4  R2)  /  {£  XP  [dM/d  t/(  4  "  R .  DAB)  ]-F}  (B-2-21) 

Kg  ,Tg/  -r!  R  =  dM/dt(Tg0-Tg,  )Cpg/(4  R2) /{EXP[CpgdM/dt/(4  -RKg)]-G)  (B-2-22) 

The  subscripts  'O'  and  *•  '  represent  quantities  at  the  sphere  surface  and 
at  infinity  respectively.  The  quantities  F  and  G  can  be  viewed  as 
correction  factors  for  the  slip  flow  regime: 

F  =  1  -  Cd  gdM/ dt  /  (n-R^g) 

G  =  1  -  CtgC  gdM/dt  /  (4rR2Kg) 

Clearly,  for  the  continuum  case,  F=G=1,  and  equations  (B-2-21)  and  (F-2- 
22)  reduce  to  those  originally  presented  by  Williams. ^4  previous 
investigators’35<  1  3 have  substituted  the  continuum  forms  of  the  flux 
expressions  (3-2-21  anc  3.-2 -22)  into  the  droplet  energy  and  mass  balances 
Equations  'B-2-U)  and  B-3-S).  Substitution  of  expression  (B-2-22)  in 
(3-2-3)  and  rearrangement  yields  the  following  quasi-steady  expression 
for  the  mass  balance. 


+ 


aO  "  a  ■'  1 '  1  “  aO  '  ^ 


y:y 

V.y. 


dM.  dt 


’  A2  ln  [(*•’ 


(  1  -F) 


( B-2  -2  3) 


Note  that  if  the  continuum  restriction  is  made  then  F  =  1,  and  Equation  P- 
2-2  3 1  is  identical  to  that  used  by  both  Sageev  and  -..'.fell1-'  an 
Armstrong.1^  In  this  pseudo-steady  approach  the  ordinary  differential 
euqations  (8-2-4)  with  (B-2-21),  and  (B-2-2  3)  provide  a  complete  descrip¬ 
tion  of  droplet  behavior.  With  the  restriction  to  continuum  i  -  h  <  • :  <  -  . 
and  some  further  simplifying  assumptions,  Armstrong1^  nas  achieved  an 
analytical  solution  of  the  system  (B-2-4),  (B-2-21),  and  (B-2-23),  which 
is  valid  for  short  heating  times  during  which  the  droplet  may  be  assumed 
to  move  quickly  to  a  steady  pool  temperature.  Sageev  and  Seinfeld1^ 
have  extended  this  to  droplets  comprised  of  aqueous  solutions. 
Computationally,  solution  the  pair  (B-2-21)  and  ( B -2  -22)  is 
straightforward.  It  should  be  noted  that  the  correction  factor  F  is 
defined  in  terms  of  the  evaporation  rate  dM/dt.  However,  this  poses  no 
problem  In  the  numerical  solution  of  the  pair  (6-2-4)  and  (B-2-23). 
Equation  (B-2-23)  is  solved  iteratively  at  each  time  step  for  the 
evaporation  rate  dM/dt. 

In  Figures  B-2-5  and  B-2-6  the  solutions  to  equations  (B-2-4), 
(B-2-21),  and  (B-2-23)  are  compared  with  our  transient  solution 
(Equations  (B-2-4)  and  (8-2-8)  through  (B-2-16)),  for  a  1  .  m  dro,  let 
irradiated  with  light  of  intensity  1C>9  W/m^. 

Previous  inv es t ig a  tor s  1  34 ,  1  35,  140  have  used  an  energy  balance 
similar  to  Equation  (B-2-4),  but  neglected  the  terms  E  q^,  and  Fp*  —  c.‘- . 

in  order  to  allow  a  direct  comparison  with  their  work,  in  curve  (  R  f 
have  held  F  =3  =  1 ,  and  ignored  the  term  s  representing  convective  -■  ■ 
loss,  E ry,  and  droplet  shrinkage,  Fp,  and  furthermore,  held  dili  ; 
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Figure  3-2-5.  Droplet  temperature  comparison  of  pseudo-steady  state 
solution  and  full  transient  solution  for  long-time  droplet 
heating.  Conditions  as  in  Figure  B-2-1 
A  -  Rigorous  transient  solution 

3  -  pseudo  steady  solution  of  previous  investigators 
C  -  pseudo  steady  solution  using  rigorous  energy  balance, 
temperature  dependent  D^g  and  slip  flow  boundary 
eond itions . 
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Droplet  si^e  comparison  of  pseudo-steady  3tate  solution 
and  full  transient  solution  for  long-time  droplet  heating, 
.'lame  as  Figure  B-2-5. 


coefficient  constant.  This  form  of  the  p  -ste  niy  so'.  u‘  :<n  provides  a 
poor  prediction  of  droplet  peak  temperature,  and  we  immediately  suspect 
this  to  be  caused  by  the  difference  in  energy  balances  used,  more  than  by 
the  nature  of  the  quasi-steady  assumption.  Inclusion  of  slip  boundary 
conditions  the  rigorous  energy  balance,  and  a  temperature 

dependent  diffusion  coefficent,  provides  an  improved  prediction  of 
droplet  behavior.  Although  each  effect  has  not  been  singled  out 
explicitly  as  in  Figure  B-2-1,  their  relative  magnitudes  are  similar  to 
those  already  shown.  While  the  temperature  profiles  (Figure  B-2-5)  are 
still  in  slight  disagreement,  the  radial  profiles  agree  closely  (Figure 
3-2-6).  When  studying  atmospheric  extinction  of  electromagnetic 
radiation,  the  size  of  the  absorbing  particle  is  of  greater  importance 
than  its  temperature  (in  the  absence  of  thermal  blooming  effects). 
Figure  B-2-6  shows  that  the  quasi-steady  assumption  may  be  made  to 
provide  an  adequate  description  of  droplet  shrinkage  only  if  the  rigorous 
energy  balance  Equation  (B-2-41  is  use  d .  Greater  accuracy  can  be 
achieved  by  further  use  of  a  temperature  dependent  diffusion  coefficient, 
and  the  slip  flow  correction  factors  F  and  G. 

Figure  3-2-7  shows  droplet  response  to  a  pulsed  laser  source. 
Here  various  sized  droplets  have  been  irradiated  with  a  r  microsecond 
pulse  of  10°  W/m--.  It  can  be  seen  that  the  smaller  drops  attain  a  quasi¬ 
steady  temperature  during  passage  of  the  pulse.  Also  shown  for 
comparison  are  results  obtained  using  the  quasi-steady  analysis  and 
neglecting  the  effects  discussed  above.  As  shown  by  Figure  B-2-1, 
neglect  of  these  effects  results  in  an  erroneous  prediction  of  droplet 
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Figure  Troplet  te  m  p-ra  t.  ur  e  response  to  a  pulsed  laser  source. 

Pulse  length  ;  5  seconds,  1q  *  10^  W/m-,  initial  radius 
indicated  is  in  m . 

A  -  full  transient  solution 

B  -  pseudo-steady  state  solution  of  previous  investiga¬ 
tors  . 

peak  temperature,  even  for  the  case  of  heating  times  very  short  in 
can  pari  son  to  the  droplet's  lifetime. 

d  .  Conclusions . 

We  have  made  a  rigorous  theoretical  development  of  the 
equations  which  describe  the  continuum  evaporation  of  a  single  droplet 
suspended  in  an  inert  host  gas,  and  extended  thi  s  a  pprox  im  ate  1  y  to  the 
noncontinuum  regime  through  a  slip  flow  analysis. 1  33  A  computer  code 
I  PIPE’1  for  numerical  solution  of  the  system  has  been  developed  and 
utilized  to  compare  this  rigorous  development  with  a  simplified  analysis 
which  uses  a  pse  di  -  -  s*  *>-;•!  y  approx  imation  to  the  hydrodynamic  equations 
governing  droplet  evaporation.  This  pseudo- steady  approximation  allows  a 


straightforward  analytical  spatial  integration  of  the  transport 
equations,  and  has  been  commonly  used  in  the  past.  Its  solution  is  also 
attractive  in  that  it  affords  a  reduction  of  several  orders  of  magnitude 
in  the  computational  time  required,  to  characterize  a  droplet's  behavior 
over  its  lifetime.  However,  the  pse  ud  o-  stead  y  solution  can  be  made  to 
agree  approximately  with  the  full  transient  solution  only  if  a  rigorous 
form  of  the  droplet  energy  balance  is  used,  as  done  here.  The  most 
important  feature  of  this  rigorous  balance  is  that  it  accounts  for  the 
rate  of  charge  of  droplet  radius.  Significant  differences  in  prediction 
of  peak  droplet  temperature  result  from  neglect  of  the  droplet  shrinkage 
term.  A  further,  although  less  important  improvement  in  the  prediction 
using  the  pseudo  steady  approximation  can  be  achieved  if  a  temperature 
dependent  diffusion  coefficient  is  used,  and  if  the  appropriate  jump 
boundary  conditions  in  temperature  and  concentration  are  included. 

For  the  case  of  a  water  droplet  irradiated  with  a  1 0.6  m  laser 
1  ig  ht  of  pu  1  se  leng  th  g  reate  r  than  the  ord  er  of  a  microsecond,  we  hav  e 
satisfied  ourselves  that  the  pseudo-steady  solution  developed  here  can  be 
made  to  give  an  adequate  qualitative  description  of  droplet  behavior.  We 
are  using  the  pseud  o-  ste  ad  y  solution  to  model  the  laser-forced 
evaporation  of  a  polyd  isper sed  aerosol  distribution.  This  involves  both 
an  approximate,  moment  method,114^  and  a  high  accuracy  finite  element  code 
developed  by  Tsang  and  Brock. 1  ^  In  this  context,  the  computational 
speed  of  the  pse  ud  o- s  te  ad  y  solution  is  especially  important.  We  will 
elaborate  or  our  results  in  a  later  publication.  With  this  pol  yd  ispersed 


model  'cmplete,  a  strsight  forwarl  extension  allows  one  to  tackle 
laser  hole-boring  problem,  modeling  the  propagation  of 
electromagnetic  beam  through  a  po  1  yd  i  spe  r  sed  aerosol  cloud,  a  prob 
which  up  to  now  appears  to  have  been  solved  for  the  limiting  case  o 
mono-dispersed  aerosol. 
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C.  AEROSOL  FILTRATION 


We  have  only  recently  started  work  on  aerosol  filtration.  Two 
types  of  filters  have  been  studied  theoretically.  The  fibrous  filter  and 
the  membrane  filter.  Here  only  the  work  on  membrane  filtration  is 
discussed . 

1 .  Membrane  Filters. 

a.  In trod uction. 

For  aerosol  filtration,  fibrous  filters  have  played  a  dominant 
role  and  so  have  received  most  attention.  By  contrast,  membrane  filters, 
although  widely  used,  have  not  been  studied  extensively.  The  structure  of 
membrane  filters  exerts  a  somewhat  limiting  influence  on  their 
applicability,  but  they  are  still  used  widely  and  they  are  useful 
especially  in  the  analysis  of  aerosols11*?  due  to  their  own  characteristic 
advantages  such  as  a  nearly  absolute  efficiency,  light  and  non- 
hygropscopic  character,  capture  of  particles  mostly  on  the  surface  or  in 
the  thin  surface  layer,  slow  clogging  of  its  pores,  wide  thermal  range, 
resistance  against  acids,  bases  and  certain  solvents,11*?  to  name  a  few. 
In  comparison  with  fibrous  filters,  membrane  filters  have  been  studied 
theoretically  very  little  since  they  have  complex  structure  and  filtra¬ 
tion  mechanisms  which  don’t  permit  simplifying  assumptions.  Therefore, 
exact  theoretical  analysis  requires  numerical  analysis. 

In  recent  times,  the  finite  element  method  has  been  widely  used 
to  solve  fluid  dynamics  problems.  It  can  handle  problems  with  non- 


regular  geometries  more  easily  than  finite  difference  methods.  In  this 
study,  a  finite  element  method  based  on  Galerkin  discretization  was  used 
to  investigate  the  membrane  filtration  problem. 

b .  Mathematical  Formulation  and  a  Finite  Element  Scheme. 

The  mathematical  description  of  steady  fluid  motion  for  the 
membrane  filter  problem  is  assumed  to  be  given  by  the  following 
equations : 

(continuity  equation) 

V  *  V  =  0  (C-1-1) 

(momentum  equation) 

p  V  *  V  V  =  p_f+divo  =  -Vp  +  p_f  +  V  •  x  (C-1  -2) 

Here  the  fluid  is  assumed  to  be  incompressible  and  Newtonian.  In  these 
equations,  V  represents  the  fluid  velocity,  p  the  density,  f  the  body 
force  vector  per  unit  mass,  x  deviatoric  stress  tensor,  a  total  stress 
tensor,  p  pressure.  For  the  constitutive  equation  a  Newtonian  fluid  was 
assumed.  Then 

x  =  2n  D  (C-1 -3) 

where  D  is  rate  of  strain  tensor.  In  the  analysis,  we  also  need  the 
diffusion  equation  for  aerosol  concentration, 

v  •  V  Cft  =  Da  v2ca  +  Rft  ( C-1  -4) 

where  CA  is  a  concentration  of  aerosol,  DA  is  the  diffusivity  of  aerosol 
in  the  fluid,  PA  is  the  rate  of  formation  of  aerosols  per  unit  volume  of 
the  fluid.  With  suitable  boundary  conditions,  these  four  equations  form 
the  basis  for  analysis  by  the  finite  element  method  used  in  this  study. 


Using  the  variational  statement  for  these  equations  implicitly  included 


with  finite  element  interpolation  for  the  independent  variables,  v,  p 
and  C  yields  the  standard  finite  element  equations.  Since  the  finite 
element  scheme  used  in  this  study  has  been  fully  described  elsewhere, ^ 
the  derivation  is  not  given  here;  we  only  mention  briefly  the  main 
features  of  this  scheme.  Basically  the  finite  element  method  code  used 
is  designed  for  steady  state,  incompressible,  two  dimensional  (plane  or 
axisymmetric  without  torsion)  problems.  It  is  based  on  the  Galerkin 
discretization  procedure,  solving  simultaneously  equations  (C-1-1)  to  (C- 
1-4)  in  their  full  nonlinear  forms.  To  handle  the  nonlinear  terms, 
iteration  is  done  until  convergence  occurs.1**®  The  computer  program  has 
been  amply  tested  for  correct  simulations. 


c.  Diffusion  of  Aerosols  in  Laminar  Flow  in  a  Cylindrical  Tube 


Before  solving  the  membrane  filtration  problem,  the  problem  of 
mass  transfer  of  aerosols  with  axial  diffusion  in  laminar  flow  through  a 
cylindrical  tube  was  solved  first  since  its  analytical  solution  as  a 
series  of  Bessel  functions  is  already  known.1  **9  This  provides  a  test  of 
our  numerical  procedures. 

When  small  particles  are  suspended  in  a  fluid  flowing  through  a 
channel,  the  random  Brownian  diffusion  of  the  particles  may  bring  them 
into  contact  with  the  channel  walls  where  they  are  absorbed  and  normally 
lose  any  electric  charge  they  carry.  The  steady  state  mass  diffusion 
equation  governing  the  concentration  of  particles  within  a  fluid 
moving  in  a  circular  pipe  is  given  by  the  diffusion  equation  (C-1-4). 


m 


Cylindrical  coordinates  (r,z)  with  the  axis  of  the  tube  at  r=0  are  used 


for  this  problem. 

Many  authors  have  studied  this  problem  under  varying 
assumptions,  e.g..  Tan  and  Hsu,  1^9»150  Ingham, 1  5 1»  1 52  Qormley  and 
Kennedy, ^53,  Davies.^1*,  155  Some  of  them  solved  this  problem  under  the 
assumption  of  a  fully  developed  or  a  flat  axial  velocity  profile  at  tube 
entry  and  some  neglected  axial  diffusion  (This  can  be  a  valid  assumption 
provided  the  Peclet  number  is  larger  than  100.  However  for  Pe  <  100,  the 
effect  of  axial  diffusion  is  still  significant  near  the  channel 
entrance 1 50) .  The  entrance  effect  in  a  short  tube  was  considered  by 
Ingham *52  With  no  axial  diffusion. 

In  this  study  we  will  solve  full  diffusion  equation,  making 
none  of  these  previous  assumptions.  With  RA  =  0,  equation  (C-1-4)  can  be 
written  as 

vz9CA  vr9CA  (3(rCA)  (92CA) 

-  +  -  +  Da  -  +  - —  (C-1-5) 

9z  3r  r9r  9r  9z^ 

with  the  following  boundary  conditions  imposed  on  this  equation 

CA  (0’r)  =  C0 

CA  (  ,r)  =  0 

acA 

-  (z,0)  =  0  (C-1-6) 

9r 

Ca(z,R)  =  0 

The  first  boundary  condition  is  the  assumption  of  uniform  concentration 
of  particles  at  the  tube  inlet  and  the  second  represents  decay  of 
concentration  to  zero  for  very  long  tubes,  the  third  one  is  the 


axisymmetric  condition  at  the  center  of  the  tube  and  the  fourth  is  the 
condition  of  perfect  absorption  at  the  tube  wall.  The  initial  velocity 
profile  was  assumed  that  for  Poiseuille  flow  vz/vq  =  2  (l-r^/R2),  where  R 
is  tube  radius.  The  schematic  diagram  of  this  problem  is  shown  in  Figure 
C-1  -1 . 


The  solution  of  this  problem  can  be  found  following  Tan  and 
Hsu1^  and  defining  a  parameter  F^  as  the  function  of  aerosols  arriving 
at  a  distance  z. 


vzCA2  7Tr  dr 


IT  R^CqVo 


(C-1 -7) 


Figure  C-1-1.  Geometry  of  diffusion  inside  tube. 


(where  y  =D^z/ vqR2 =2z/RPe) 

The  following  solution  was  derived  by  Tan  and  Hsu11^ 


F  =  _jj  I  Dn  { Rn * ( 1 )  +  6n2  /  1  ORndn}  exp(-Bn2  )  (C-l-8) 

M  n=1  0 

where  Rj,  is  the  eigen  function  of 

V  +  R n'/y  +  Bn2  C(1-n2>  +  (Bn2/Pe2) ]Rn=0  (C-1-9) 

with  Rn'(0)=0  and  Rn(1)=0,  Bn  is  the  eigenvalue  of  equation  (C -1 -5) , 

£  =  y/2,  ?i=  r/R  and  Pe=2Rvg/D£.  The  values  of  were  tabulated  by  Tan 

and  Hsu^O  fOP  pe=1,10, 100. 

To  solve  this  problem  by  a  finite  element  method,  a  total  of 
120  elements  (6x20)  was  used.  Calculated  values  of  are  tabulated  for 
Pe=1,10,100  in  Tables  C-1-1,  C-1  -2,  and  C— 1  —3  and  Tan  and  Hsu's  result  is 
also  presented  for  comparison.  Generally  the  agreement  is  very  good. 
The  differences  between  the  FEM  solution  and  those  of  Tan  and  Hsu's  are 
less  than  0.5*.  The  integration  of  equation  (C-1 -7)  was  done  using 


Table  C-1-1.  Values  of  F  for  Diffusion  of  Aerosols  in 

y 

a  Tube  When  Pe  =  1 


F^  by  FEM 

F^  by  Tan  and  Hsu 

0 

0. 99564  9 

1. 

0.005 

0.984864 

0.9885 

0.01 

0.974130 

0.9773 

0.02 

0.952  884 

0  .  95  5  3 

0.05 

0.891686 

0.8928 

0. 1 

0.798612 

0.7991 

0.2 

0.642655 

0.6431 

0.3 

0.519083 

0.5196 

0.4 

0.420359 

0.  421 C 

0.5 

0.341043 

0.3416 

Table  C-1-2.  Values  of  Fy  for  Diffusion  of  Aerosols  in 

a  Htbe  When  Pe  =  10 


F  ,  by  ?m 

Fy  by  Tan  and  Hsu 

0 

0.995649 

1. 

0.005 

0.990372 

0.994 

0.01 

0.995085 

0.994 

0.02 

0.974554 

0.977 

0.05 

0.944011 

0.9542 

0.1 

0.897130 

0.8975 

0.2 

0.81616 

0.8166 

0.3 

0.748123 

0.7488 

0.4 

0.689215 

0.6901 

0.5 

0.63723 

0.6382 

1 . 

0.442648 

0.444 

2. 

0.222972 

0.224 

3. 

0. 113101 

0.  114 

4. 

0.057516 

0.058 

5. 

0  .  0292  5  6 

0.0295 

Table  C— 1  —3 -  Values  of  Fy  for  Diffusion  of  Aerosols  in 
a  Tube  When  Pe  =  100 


Fy  by  FEM 

Fy  by  Tan  and  Hsu 

0 

0.995649 

1. 

0.01 

0.992  5  8 

0.9955 

0.02 

0.989477 

0.992 

0.05 

0.980677 

0.9824 

0.1 

0.967846 

0.9689 

0.2 

0.946274 

0.947 

0.3 

0.928518 

0.929 

0.4 

0.912596 

0.9132 

0.5 

0.898638 

0.8989 

1 . 

0.840628 

0.8406 

2. 

0.754646 

0.7545 

3. 

0.688676 

0.6872 

4. 

0.632621 

0.6305 

5. 

0.584127 

0.581 

10. 

0.396452 

0.  3968 

Simpson’s  rule  in  the  FEM  calculation,  due  to  the  coarse  grid  in  a  row  a 
small  error  can  arise.  When  y=0,  F  should  be  1,  but  the  FEM  solution 

y 

shows  small  deviations  from  this.  If  we  use  more  elements  in  a  row, 
better  results  are  obtained  but  at  increased  computation  cost.  The  small 
difference  can  be  reduced,  so  we  have  shown  that  the  FEM  solution  is 
accurate  enough  to  be  accepted.  From  these  tables  we  can  see  also  that 
decreasing  Pe  increases  F^  since  a  decrease  of  Peclet  number  means 
physically  an  enhancement  of  axial  diffusion  over  radial  diffusion, 
reducing  the  loss  of  particles  to  the  tube  walls  and  increasing  F 
values . 

d .  The  Computational  Results  of  Membrane  Filtration  and  Discussion. 

With  assurance  to  the  accuracy  of  the  present  code,  membrane 
filtration  was  investigated  next.  Since  there  are  so  many  pores  on  the 
surface  of  membrane  filters,  it  is  impossible  to  include  all  the  pores  in 
the  simulation.  Instead,  we  do  the  calculations  using  a  single  pore 
model,  a  procedure  analogous  to  that  used  in  the  fibrous  filter  problem. 

Consider  the  membrane  filtration  process  as  occurring  in  a 
single  pore  membrane  as  shown  in  Figure  C-1-2.  Dimensionless  variables 
are  used : 

r*  r  r/R,  z*  =  z/R 

C*=Ca/Cq  v«=v/v0 

u*  =  u/v0  DA»  =  DA/2RV0 

where  R  i3  the  reference  length,  vq  is  reference  velocity. 


d=0 

Figure  C-1-2.  Geometry  for  membrane  filter. 

r-|Z:  radial  coordinates;  u*,  v*  radial  and  axial  veloci¬ 

ties; 

C:  concentration;  Rm:  outer  radius  of  pore. 

The  following  boundary  conditions  are  used: 
v*  =  1,  u*  =0,  C*  =  1  at  far-  upstream  ( z*  =  zmin) 

u*  =  v*  =  0,  C*  =  0  at  pore  wall  and  membrane  surface 

u*  =  0  and  C*/  r*  =  0  at  pore  center  line.  (C-1-10) 

The  far  upstream  position  was  taken  as  3  times  the  pore  radius. 
The  pore  radius  was  set  as  follows.  If  A  is  the  area  of  the  whole 
filter,  L  the  thickness,  Ap  the  area  of  one  pore  and  Np  the  number  of 
pores  per  unit  surface  area,  then  the  volume  of  the  pores  Vp  is 

Vp  =  ApNpLA  (C-1-11) 

Since  the  volume  of  the  filter  Vf  is  Vf=AL,  the  porosity  =Vp/Vf,  then 
equation  (C-1-11)  yields 

c  =  ApNp  (C-1-12) 

The  porosity  can  be  expressed  as  a  function  of  pore  diameter  Dp  and  the 
axial  distance  between  pores  d  from  simple  geometrical  considerations11*? 

t  =  7’Dp2/(2/"3(j2)  ( C  — 1  —13) 
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The  porosity  reaches  the  maximum  value  for  d  =  Dp»  i.e., 
porosity  of  membrane  filters  is  accordingly  much  lower  than  the  porosity 
of  fibrous  filters.  In  the  limit  of  a  single  pore  membrane,  the  porosity 
is  e=1/Rm2»  so  from  equations  (  C— 1 —12)  and  ( C— 1 —1 3 ) 

Rm  =  {^2/(2/Id2)  }0.5  (C-1-m) 

Hence  minimum  R^  reaches  almost  1.05  when  the  membrane  has  the  maximum 
porosity . 

In  this  study,  four  different  values  of  Rm  were  taken:  1.05, 
1.1,  1.2,  1.3  respectively,  and  accordingly  the  porosity  was  90.7$, 

82.64$,  69.44$,  59.17$.  In  this  study  two  different  mechanisms  of 

membrane  filtration  were  investigated:  diffusional  deposition  inside  the 
pore  and  inertial  deposition  on  the  filter  surface.  For  diffusion  inside 
the  pore,  the  efficiency  was  defined  as 

Ed  =  1  "  Cf«/q»  ( C— 1  —15) 

where  Cj*  is  the  concentration  at  pore  inlet,  q*  is  the  concentration  at 
the  pore  outlet.  For  inertial  deposition,  particle  trajectories  were 
calculated  based  on  the  equation  of  particle  motion  following  Stoke' s  law 

du  * 

Stk  — —  +  Uf*  -  u  * 
dt  P 

dVpi 

Stk  -  +  Vf*  _  yD*  (C-1-16) 

dt  v 

where  Stk  is  Stokes  number  defined  as 

Stk  =  mv0/3TT  n0dpR  (C-1  -17) 

In  the  above  equations,  Uf*  and  Vf*  are  the  radial  and  axial  fluid 
velocities.  Up  and  Vp  are  radial  and  axial  particle  velocities,  dp  is  the 
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particle  diameter,  nQ  is  the  fluid  viscosity,  and  m  is  the  mass  per  unit 
particle.  Since  particle  movement  does  affect  fluid  momentum,  equation 
(C-1-16)  is  solved  using  the  Runge-Kutta-Fehlberg  method  with  velocity 
fields  generated  from  the  FEM  solution.  Far  upstream,  particles  were 
assumed  to  have  the  same  velocities  as  that  of  the  fluid. 


Inertial  deposition  efficiency  wa3  defined  as  follows 

Si  =  (R22  -  R12)/Rm2 


CC-1  -18) 


where  R2  is  the  starting  radial  position  of  a  particle  far  up  stream 
which  hits  the  membrane  skirt  (r*  =  Rm)  and  R ^  is  the  starting  radial 
position  of  a  particle  far  upstream  viiich  hits  the  pore  skirt  (r*=l) . 

Figure  C-1-3  shows  a  fluid  stream  line  and  Figure  C  —1  — M  shows 
typical  particle  trajectories.  Figure  C-1 -5  shows  the  axial  velocity 
profiles  at  some  positions.  We  can  see  that  at  the  pore  inlet  the  axial 
velocity  is  neither  fully  developed  nor  in  plug  flow.  However  the  flow 
develops  very  fast,  so  that  except  for  very  short  pores  (very  thin 
membranes)  the  general  assumption  of  a  fully  developed  velocity  profile 
seem s  to  be  rea sonable.  Table  C -1 sho ws  the  variation  of  d i f fusion 
efficiency  and  inertial  efficiency  with  different  porosity  membranes.  As 
the  porosity  decreases,  the  inertial  deposition  increases  due  to 
increased  membrane  surface  area  and  the  initial  concentration  at  pore 
inlet  decreases  due  to  reduced  porosity.  The  diffusion  efficiency 
increases  slightly  owing  to  decreased  inlet  axial  velocity  (slow 
convection),  but  due  to  decreased  initial  concentration  the  total 
deposition  by  diffusion  decreases.  In  this  computation  when  exact 
collision  at  pore  inlet  and  outer  membrane  skirt  did  not  occur. 
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Figure  C-1-5.  Iso-axial  velocity  profiles,  ^=1.3. 


interpolation  was  done  to  get  the  inertial  deposition  efficiency.  This 
can  introduce  a  very  slight  error  in  the  calculation  of  inertial 
deposition,  but  it  i3  at  least  qualitatively  correct.  Another  fact  to  be 
mentioned  is  that  the  particle  which  moves  into  the  pore  does  not  deposit 
on  the  pore  wall,  due  to  the  laminar  flow  of  the  fluid  inside  the  pore. 
This  is  confirmed  by  particle  trajectory  calculations;  the  particle  which 
enters  the  pore  just  passes  out  the  pore. 

The  computational  result  in  Table  C  —  1  —4  is  for  Pe  =  1.  Peclet 
number  was  varied  to  determine  the  concentration  change.  The  result  is 
summarized  in  Table  C-1-5.  Increasing  Pe  increases  convection.  However 
due  to  increased  axial  diffusion,  the  diffusion  efficiency  was  reduced. 
The  total  diffused  concentration  at  the  wall  is  also  presented  in  Ihble 
C-1-5.  Increasing  Pe  raises  the  diffused  concentration  but  when  Pe  goes 
over  100  it  drops.  This  means  that  axial  diffusion  is  very  important  for 
diffusion  in  the  pore  when  the  fluid  moves  slowly  or  diffusivity  of  the 
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particles  is  very  large.  Figures  C-l  -6  through  C-1  -8  show  iso-concentra¬ 
tion  profiles  when  Pe  =  1, 10,100  respectively.  These  figures  show  that,  as 


already  mentioned,  high  Pe  causes  high  concentrations  at  the  pore  inlet 
(z*=0);  the  higher  Pe  is,  the  steeper  the  concentration  gradient  is  near 
the  pore  wall. 

Up  to  this  point,  we  have  just  looked  at  the  problem  with  the 
initial  velocity  v*  =  1.  Next  the  initial  velocity  far  upstream  was  varied 
as  3,5,10,  so  that  the  Peclet  number  was  also  changed  3»5, 10.  The  result 
is  shown  in  Table  C-1 -6.  As  previously  seen,  increasing  Pe  decreases  the 
diffusion  efficiency.  When  the  initial  axial  velocity  is  10,  inertial 
deposition  rises  very  sharply,  while  the  concentration  at  pore  inlet 
changes  remarkably  due  to  the  large  convection. 

In  the  calculation  of  particle  trajectories,  the  particle  size 
was  assumed  to  be  uniform.  For  different  sizes  of  particles,  the 
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diffusivity  and  Stokes  number  change.  The  Stokes  number  is  proportional 
to  the  square  of  particle  diameter  while  diffusivity  decreases  inversely 
with  particle  d iameter. 1 ^6  To  examine  the  selective  characteristics , 


t 


a 


which  can  be  defined  as  an  efficiency  dependence  of  the  filter  on  the 


particle  size,  particle  size  was  changed  and  accordingly  the  Stokes 


number  and  diffusivity  were  changed.  The  computational  results  are 


presented  in  Table  C -1 -7.  As  expected,  the  diffusion  efficiency 


decreases  due  to  increased  Pe  and  inertial  efficiency  increases.  Total 


efficiency  increases  slightly  with  increasing  particle  size.  From  this 


limited  calculation,  it  is  not  quite  clear  if  the  effect  of  decreasing  of 
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Figure  C-1-6.  Iso-concentration  profiles  when  Pe=1,  Rm= 1  - 3 • 
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Figure  C-1-7.  Iso-concentration  profiles  when  Pg=io,  ^=1.3 
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Figure  C-1-8.  Iso-concentration  profiles  when  Pe=100,  Rm=  1*3* 


Table  C— 1 —6 -  The  Effect  of  Initial  Velocity  on  Diffusion  and 
Inertial  Efficiency  when  R_=1.2,  Stk=0.2 


Table  C-1-7.  The  Effect  of  Particle  Size  on  Diffusion  and 
Inertial  Efficiency,  ^=1.2 


V* 

Diffusion 

Efficiency 

Inertial 

Efficiency 

F  at 
Pore  Inlet 

1. 

0.8949 

7.257E-4 

0.006096 

2. 

0.8925 

0.084998 

0.007324 

4. 

0.8880 

0.124378 

0.009582 

7. 

0.8809 

0.170056 

0.012407 

10. 

0.8737 

0.179799 

0.014742 

20. 

0.8483 

0. 184490 

0.020805 

diffusion  efficiency  is  greater  than  that  of  increasing  inertial 
efficiency.  And  it  is  not  clear  that  there  should  exist  a  minimum  point 
of  total  efficiency  at  some  particle  radius,  as  Pich  insisted J1*?  This 
point  needs  more  careful  investigation. 

e.  Conclusions. 

Using  a  finite  element  method,  membrane  filtration  efficiency 
has  been  studied.  First  diffusional  deposition  in  a  cylinder  was 
investigated  and  the  results  were  compared  to  analytical  results.  The 
agreement  was  very  good.  Using  the  tested  FEM  code,  computations  were 
carried  out  next  for  membrane  filters.  The  following  conclusions  wer 
made  for  membrane  filtration. 

(1)  Using  the  same  size  particle,  increasing  membrane  area  (or 
decreasing  porosity)  increases  inertial  deposition  efficiency  and 
slightly  increases  the  diffusion  efficiency. 


(2)  Increasing  flow  velocity  increases  Inertial  deposition 
efficiency  but  decreases  diffusion  efficiency. 

(3)  Increasing  particle  size  increases  inertial  deposition 
efficiency  but  decreases  diffusion  efficiency.  The  total  efficiency  did 
not  show  a  minimum  over  computed  range. 

(4)  Increasing  Pe  results  in  high  concentration  at  the  pore 
inlet  but  increasing  convective  flow  decreases  diffusion  efficiency  in 
the  pore. 

This  work  is  being  continued  and  compared  with  available 
experimental  data.  In  addition,  the  transmission  problem  is  being 
investigated  for  fibrous  filters  for  both  spherical  and  nonspherical 


particles 
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